Hi there,
I tried to use the example codes in plc-1.0/bin directory to calculate the lensing likelihood given an input file as explained in the clik.pdf documentation, that is
./clik_example_C /path/to/lensing_likelihood_v4_ref.clik_lensing clfile
So far these are the combinations I tried for the clfile: an ASCII file with 2049 rows (from 0 to lmax=2048) with monopole and dipole set to zero, and with the following column(s)
1) C_l^{ΦΦ}
2) [l(l+1)]^2*C_l^{ΦΦ}/2π
3) l C_l^{ΦΦ}
4) l [l(l+1)]^2*C_l^{ΦΦ}/2π
In the first two cases I get complaints saying the 4098 entries are expected, in the last two no matter the input cosmology I always get an unreasonable value
$ ./clik_example_C ../lensing_likelihood_v4_ref.clik_lensing test.dat
Checking lensing likelihood '../lensing_likelihood_v4_ref.clik_lensing' on test data. got -5.81039
Log likelihood for file test.dat : -1.44273e+60
I also ran the script clik_get_selfcheck to gather some more info about the input clfile format but apparently it requires h5py and I have problem installing it.
Has anyone ever used these example scripts and got them to work?
Thanks in advance,
Matteo
Get likelihood using clik_example_(py,C,f90)
-
- Posts: 1
- Joined: March 18 2015
- Affiliation: Dark Cosmology Centre