[CosmoMC] Failure to compute Planck likelihood in test case
Posted: July 24 2017
I can compile both the Planck likelihood and cosmomc source sucessfully, but when running the provided CosmoMC test file, it does not return the predicted value for the likelihood
I have tested the Planck likelihood using the code provided as part of the plc, and these do return the predicted values. I have also made the patch suggested here.[davidparkinson@polaris CosmoMC]$ mpirun -np 1 ./cosmomc_debug test.ini
Number of MPI processes: 1
file_root:test
Random seeds: 28259, 25208 rand_inst: 1
Using clik with likelihood file
./data/clik/hi_l/plik/plik_dx11dr2_HM_v18_TT.clik
----
clik version 6dc2a8cf3965
smica
Checking likelihood './data/clik/hi_l/plik/plik_dx11dr2_HM_v18_TT.clik' on test data. got -380.979 expected -380.979 (diff -8.68545e-09)
----
TT from l=0 to l= 2508
Clik will run with the following nuisance parameters:
A_cib_217
cib_index
xi_sz_cib
A_sz
ps_A_100_100
ps_A_143_143
ps_A_143_217
ps_A_217_217
ksz_norm
gal545_A_100
gal545_A_143
gal545_A_143_217
gal545_A_217
calib_100T
calib_217T
A_planck
Using clik with likelihood file
./data/clik/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5c_Ap.clik
BFLike Ntemp = 2876
BFLike Nq = 1407
BFLike Nu = 1407
BFLike Nside = 16
BFLike Nwrite = 32393560
cls file appears to have 5+ columns
assuming it is a CAMB file with l, TT, EE, BB, TE
info = 0
----
clik version 6dc2a8cf3965
bflike_smw
Checking likelihood './data/clik/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5c_Ap.clik' on test data. got -1e+30 expected -5247.87 (diff 1e+30)
----
TT from l=0 to l= 29
EE from l=0 to l= 29
BB from l=0 to l= 29
TE from l=0 to l= 29
Clik will run with the following nuisance parameters:
A_planck
Using clik with likelihood file
./data/clik/low_l/commander/commander_rc2_v1.1_l2_29_B.clik
----
clik version 6dc2a8cf3965
gibbs_gauss 1478fb2d-28fa-49ac-a8ae-677dbdc3600a
Checking likelihood './data/clik/low_l/commander/commander_rc2_v1.1_l2_29_B.clik' on test data. got -7.32304 expected -7.32304 (diff -2.52096e-10)
----
TT from l=0 to l= 29
Clik will run with the following nuisance parameters:
A_planck
Doing non-linear Pk: F
Doing CMB lensing: T
Doing non-linear lensing: T
TT lmax = 2508
EE lmax = 2500
ET lmax = 2500
BB lmax = 2500
PP lmax = 2500
lmax_computed_cl = 2508
Computing tensors: F
max_eta_k = 14000.00
transfer kmax = 5.000000
adding parameters for: lowl_SMW_70_dx11d_2014_10_03_v5c_Ap
adding parameters for: smica_g30_ftl_full_pp
adding parameters for: commander_rc2_v1.1_l2_29_B
adding parameters for: plik_dx11dr2_HM_v18_TT
Fast divided into 1 blocks
21 parameters ( 7 slow ( 0 semi-slow), 14 fast ( 0 semi-fast))
Time for theory: 5.66975
Time for lowl_SMW_70_dx11d_2014_10_03_v5c_Ap: 1.270103454589844E-002
loglike chi-sq
0.000 0.000 CMB: lensing = smica_g30_ftl_full_pp
0.000 0.000 CMB: plik = plik_dx11dr2_HM_v18_TT
0.000 0.000 CMB: lowTEB = lowl_SMW_70_dx11d_2014_10_03_v5c_Ap
0.000 0.000 CMB: lowl = commander_rc2_v1.1_l2_29_B
Test likelihoods done, parameter point rejected (logZero or outside prior)
Total time: 45 ( 0.01246 hours )