https://github.com/julienpeloton/lensingbiases
Computation of N0 and N1 biases in the context of CMB weak lensing
https://github.com/julienpeloton/lensingbiases
biases lensing
Last synced: over 1 year ago
JSON representation
Computation of N0 and N1 biases in the context of CMB weak lensing
- Host: GitHub
- URL: https://github.com/julienpeloton/lensingbiases
- Owner: JulienPeloton
- License: gpl-3.0
- Created: 2016-12-02T17:17:21.000Z (over 9 years ago)
- Default Branch: master
- Last Pushed: 2020-05-20T14:01:48.000Z (about 6 years ago)
- Last Synced: 2025-01-16T04:10:44.392Z (over 1 year ago)
- Topics: biases, lensing
- Language: Fortran
- Homepage:
- Size: 1.39 MB
- Stars: 7
- Watchers: 2
- Forks: 4
- Open Issues: 1
-
Metadata Files:
- Readme: README.md
- License: LICENSE
Awesome Lists containing this project
README
Lensingbiases
==
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.
### Before starting
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:
```bash
BIASPATH=/path/to/the/package
export PYTHONPATH=$PYTHONPATH:$BIASPATH
```
### Compilation of Fortran
The 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:
```bash
python setup.py build_ext --inplace --fcompiler=gfortran
```
Depending 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:
```bash
export OMP_NUM_THREADS=n
```
### Usage
User 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).