Get likelihood using clik_example_(py,C,f90)

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
Matteo Cataneo
Posts: 1
Joined: March 18 2015
Affiliation: Dark Cosmology Centre

Get likelihood using clik_example_(py,C,f90)

Post by Matteo Cataneo » March 18 2015

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

Post Reply