CosmoCoffee Forum Index CosmoCoffee

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

Error in CosmoMC regarding BAODR12 dataset
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Alvaro Hu



Joined: 19 Aug 2016
Posts: 2
Affiliation: University of Dallas

PostPosted: August 19 2016  Reply with quote

When I run jobs using the BAO.ini in batch 2 with BAODR11 and test the likelihoods, I receive this:

loglike chi-sq
6.079 12.157 CMB: lensing = smica_g30_ftl_full_pp
581.392 1162.783 CMB: plik = plik_dx11dr2_HM_v18_TT
5249.288 10498.575 CMB: lowTEB = lowl_SMW_70_dx11d_2014_10_03_v5cAp
7.639 15.278 CMB: lowl = commander_rc2_v1.1_l2_29_B
5.439 10.878 Hubble: H073p0
65.299 130.597 SN: Union2.1
0.044 0.087 BAO: 6DF
0.463 0.927 BAO: MGS
1.428 2.856 BAO: DR11CMASS
0.618 1.236 BAO: DR11LOWZ

But when I run the BAO.ini in batch 3, with BAODR12, the more recent dataset, I receive this:

loglike chi-sq
6.079 12.157 CMB: lensing = smica_g30_ftl_full_pp
0.000 0.000 CMB: plik = plik_dx11dr2_HM_v18_TT
5249.288 10498.575 CMB: lowTEB = lowl_SMW_70_dx11d_2014_10_03_v5cAp
0.000 0.000 CMB: lowl = commander_rc2_v1.1_l2_29_B
0.000 0.000 Hubble: H073p0
0.000 0.000 SN: Union2.1
0.000 0.000 BAO: 6DF
0.463 0.927 BAO: MGS
0.110 0.220 BAO: WIGGLEZ = WiggleZ
0.000 0.000 BAO: DR12CMASS
0.459 0.917 BAO: DR12LOWZ

where a lot of the likelihoods are 0.000. The only thing I changed was going from BAODR11 to BAODR12. Also, when I run the chains for BAODR12, I receive an error saying that it could not start after 1000 tries, but it works for BAODR11.
Back to top
View user's profile  
Antony Lewis



Joined: 23 Sep 2004
Posts: 1243
Affiliation: University of Sussex

PostPosted: August 25 2016  Reply with quote

Looks like something more general messed up, e.g. try a clean build (and check plik tests return sensible likelihoods when run separately).
Back to top
View user's profile [ Hidden ] Visit poster's website
María Pía Piccirilli



Joined: 17 Sep 2016
Posts: 5
Affiliation: Facultad de Ciencias Astronómicas y Geofísicas - Universidad de La Plata

PostPosted: September 17 2016  Reply with quote

Hi everyone!
I'm having the same trouble with BAO datasets.
Plank data alone or together with BAODR11 runs perfectly, but when I try BAODR12 in batch3 I get the error message: "MCMC.f90: Couldn't start after 1000 tries - check starting ranges".
I've tried a clean build with the same results as before.
Here I copy the output of the likelihood test, having exactly the same problems as cited in the first post by Alvaro.
Any clue?
Thanks in advance!

plik TT + lowTEB:
loglike chi-sq
581.156 1162.311 CMB: plik = plik_dx11dr2_HM_v18_TT
5249.287 10498.574 CMB: lowTEB = lowl_SMW_70_dx11d_2014_10_03_v5cAp
Test likelihoods done, total logLike, chi-eq = 5830.568 11661.135

plikTT + lowTEB + BAO DR11:
loglike chi-sq
581.156 1162.311 CMB: plik = plik_dx11dr2_HM_v18_TT
5249.287 10498.574 CMB: lowTEB = lowl_SMW_70_dx11d_2014_10_03_v5cAp
0.044 0.088 BAO: 6DF
0.463 0.927 BAO: MGS
0.110 0.221 BAO: WIGGLEZ = WiggleZ
1.431 2.861 BAO: DR11CMASS
0.620 1.240 BAO: DR11LOWZ
Test likelihoods done, total logLike, chi-eq = 5833.236 11666.472

plikTT + lowTEB + BAO DR12:
loglike chi-sq
0.000 0.000 CMB: plik = plik_dx11dr2_HM_v18_TT
5249.287 10498.574 CMB: lowTEB = lowl_SMW_70_dx11d_2014_10_03_v5cAp
0.000 0.000 BAO: 6DF
0.000 0.000 BAO: MGS
0.110 0.221 BAO: WIGGLEZ = WiggleZ
0.000 0.000 BAO: DR12CMASS
0.000 0.000 BAO: DR12LOWZ
Test likelihoods done, parameter point rejected (logZero or outside prior)
Back to top
View user's profile  
Alvaro Hu



Joined: 19 Aug 2016
Posts: 2
Affiliation: University of Dallas

PostPosted: September 18 2016  Reply with quote

What I ended up doing was replacing the BAODR12.ini in the job .ini file to BAODR12_RSD.ini, which allowed it to run. Hopefully that works for you as well
Back to top
View user's profile  
Antony Lewis



Joined: 23 Sep 2004
Posts: 1243
Affiliation: University of Sussex

PostPosted: September 19 2016  Reply with quote

Check you have the latest version from github, cf.

https://github.com/cmbant/CosmoMC/issues/6
Back to top
View user's profile [ Hidden ] Visit poster's website
María Pía Piccirilli



Joined: 17 Sep 2016
Posts: 5
Affiliation: Facultad de Ciencias Astronómicas y Geofísicas - Universidad de La Plata

PostPosted: September 21 2016  Reply with quote

Thank you very much Alvaro and Antony for your reply!

I've tried both suggestions. Working with RSD files leads me to something strange with the use of memory, as for running just a likelihood test requires more than 60gb which leads to the job beeing killed by the cluster queue manager.

Replacing data/sdss_DR12CMASS_bao.dataset and adding data/sdss_DR12CMASS_bao_dilation_qpm.dat and data/sdss_DR12CMASS_bao_dilation_md-patchy.dat downloaded from https://github.com/cmbant/CosmoMC/pull/5/commits/6e7b4d66ada21cb7cd529d34f11f873b3c4a0e49 didn't work at first. I tried getting the new version of bao.f90 and recompiling, but I get:
bao.f90(27): error #6366: The shapes of the array expressions do not conform. [MEASUREMENT_TYPES]
character(LEN=Ini_Enumeration_Len), parameter :: measurement_types(8)

What am I missing?
Thanks!
Back to top
View user's profile  
Antony Lewis



Joined: 23 Sep 2004
Posts: 1243
Affiliation: University of Sussex

PostPosted: September 21 2016  Reply with quote

I would just do a new full clone of the latest git repository rather than trying to get individual files.
Back to top
View user's profile [ Hidden ] Visit poster's website
María Pía Piccirilli



Joined: 17 Sep 2016
Posts: 5
Affiliation: Facultad de Ciencias Astronómicas y Geofísicas - Universidad de La Plata

PostPosted: September 23 2016  Reply with quote

Thank you very much, Antony!
A full clone worked perfectly!
Back to top
View user's profile  
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.