CosmoCoffee Forum Index CosmoCoffee

 
 FAQFAQ   SearchSearch  MemberlistSmartFeed   MemberlistMemberlist    RegisterRegister 
   ProfileProfile   Log inLog in 
Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

[CosmoMC] Failure to compute Planck likelihood in test case
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
David Parkinson



Joined: 19 Jan 2005
Posts: 22
Affiliation: KASI

PostPosted: July 24 2017  Reply with quote

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
Quote:

[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:
Acib_217
cib_index
xi_sz_cib
Asz
ps_A100_100
ps_A143_143
ps_A143_217
ps_A217_217
ksz_norm
gal545_A100
gal545_A143
gal545_A143_217
gal545_A217
calib_100T
calib_217T
Aplanck
Using clik with likelihood file
./data/clik/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5cAp.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_v5cAp.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:
Aplanck
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:
Aplanck
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_v5cAp
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_v5cAp: 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_v5cAp
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 )


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.
Back to top
View user's profile [ Hidden ] Visit poster's website
David Parkinson



Joined: 19 Jan 2005
Posts: 22
Affiliation: KASI

PostPosted: August 02 2017  Reply with quote

I did manage to solve this. Cython was not being properly install by the waf tool that came with the Planck likelihood, as the resource location they had specified for cython in the wscript was out of date. They suggested that it should be downloaded from http://cython.org/release/Cython-0.14.1.tar.gz, but this link is dead, and cython is now available from https://pypi.python.org/pypi/Cython/.

Someone please ask the Planck team to update wscript accordingly.
Back to top
View user's profile [ Hidden ] Visit poster's website
Display posts from previous:   
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software All times are GMT + 5 Hours
Page 1 of 1

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum


Powered by phpBB © 2001, 2005 phpBB Group. Sponsored by WordWeb online dictionary and dictionary software.