Computation of N0 and N1 biases in the context of CMB weak lensing (flat-sky approximation). Derivatives of N1 biases with respect to the lensing potential power-spectrum available.
License: GNU License (see the LICENSE file for details) covers all files in the lensingbiases repository unless stated otherwise.
The Python code has the following dependencies:
- numpy, pylab, argparse
- f2py
In most machines, those packages are automatically installed. Make sure you update your PYTHONPATH to use the code. Just add in your bashrc:
BIASPATH=/path/to/the/package
export PYTHONPATH=$PYTHONPATH:$BIASPATHThe package is written in python, but internal operations are done in Fortran. Therefore you need to compile it before using the routines. We provide a setup file. Just run:
python setup.py build_ext --inplace --fcompiler=gfortranDepending on system, you may need to specify a different fortran compiler. For the purist, we also provide a Makefile for a direct compilation. For extended multipole ranges, computation can be long. So make sure that you are using several processors by adding in your bashrc:
export OMP_NUM_THREADS=nUser has to provide the range of multipole, beam width, noise level, etc. In order to ease the start, we provide both a test script and a launcher:
- LensingBiases.py (see the main)
- launch.sh
For usage on clusters, we also provide a batch (SLURM).