CosmoCoffee Forum Index CosmoCoffee

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

issue with CosmoMC chains when compiled with gfortran
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Elizabeth Gould



Joined: 30 Jan 2015
Posts: 10
Affiliation: University of Waterloo; Perimeter Institute for Theoretical Physics

PostPosted: May 04 2017  Reply with quote

When I run CosmoMC with check-pointing, stop running it, and start again later from the checkpoint, the code seems to have issues and stops. In addition, some of the chain files appear to stop around 20 lines in length, and GetDist fails to read them.

I currently have the Planck Likelihood Code installed with gfortran and CosmoMC compiled with gfortran. I had not had this issue when I used ifort. Due to an issue with the PLC, it had to be reinstalled, however, and it was reinstalled with gfortran.

Here are two samples of the last things written to the screen by CosmoMC:

Quote:
Fast divided into 1 blocks
35 parameters ( 9 slow ( 0 semi-slow), 26 fast ( 0 semi-fast))
1 Reading checkpoint from chains/CDM01_1.chk
Checking likelihood './data/clik/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5cAp.clik' on test data. got −5247.87 expected −5247.87 (diff 3.89829e−07)
—-
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
starting Monte-Carlo
Checking likelihood './data/clik/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5cAp.clik' on test data. got −5247.87 expected −5247.87 (diff 3.89829e−07)
—-
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
Chain:11 drag accpt: 0.379746825 fast/slow 121.696205 slow: 79
Chain:3 drag accpt: 0.333333343 fast/slow 147.111115 slow: 90
Chain:7 drag accpt: 0.326086968 fast/slow 146.782608 slow: 92
Chain:5 drag accpt: 0.288461536 fast/slow 138.500000 slow: 104
Chain:1 drag accpt: 0.300000012 fast/slow 149.020004 slow: 100
Chain:0 drag accpt: 0.441176474 fast/slow 133.507462 slow: 67
Chain:9 drag accpt: 0.229007632 fast/slow 153.236649 slow: 131
TObjectList: Unknown type to save



Quote:
Fast divided into 1 blocks
36 parameters (10 slow ( 0 semi-slow), 26 fast ( 0 semi-fast))
1 Reading checkpoint from chains/CDM03_1.chk
info = 0
info = 0
starting Monte-Carlo
—-
clik version 723c1a4b0580 MAKEFILE
bflike_smw
—-
clik version 723c1a4b0580 MAKEFILE
bflike_smw
Checking likelihood './data/clik/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5cAp.clik' on test data. got −5247.87 expected −5247.87 (diff 3.83458e−07)
—-
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
Checking likelihood './data/clik/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5cAp.clik' on test data. got −5247.87 expected −5247.87 (diff 3.83458e−07)
—-
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
Chain:9 drag accpt: 0.454545468 fast/slow 147.363632 slow: 66
Chain:0 drag accpt: 0.348837197 fast/slow 126.522385 slow: 67
Chain:10 drag accpt: 0.434782594 fast/slow 133.014496 slow: 69
TObjectList: Unknown type to save
Back to top
View user's profile  
Antony Lewis



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

PostPosted: May 08 2017  Reply with quote

Which gfortran version is this? (not sure I've ever tested checkpointing with gfortran myself).
Back to top
View user's profile [ Hidden ] Visit poster's website
Elizabeth Gould



Joined: 30 Jan 2015
Posts: 10
Affiliation: University of Waterloo; Perimeter Institute for Theoretical Physics

PostPosted: May 08 2017  Reply with quote

This is the version listed:
Quote:
gfortran –version
GNU Fortran (Ubuntu 4.9.4−2ubuntu1~14.04.1) 4.9.4
Back to top
View user's profile  
Antony Lewis



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

PostPosted: May 08 2017  Reply with quote

You need at least version 6 (I'm amazed it even compile with v4.9..)
Back to top
View user's profile [ Hidden ] Visit poster's website
Elizabeth Gould



Joined: 30 Jan 2015
Posts: 10
Affiliation: University of Waterloo; Perimeter Institute for Theoretical Physics

PostPosted: May 18 2017  Reply with quote

Thanks.

It appears that I did use gfortran 6.2. However, since both are installed it may have compiled partly with gfortran 4.9. I have recompiled it ensuring that the later version of gfortran is used for everything. I am currently rerunning it to see if that fixed the issue.
Back to top
View user's profile  
Elizabeth Gould



Joined: 30 Jan 2015
Posts: 10
Affiliation: University of Waterloo; Perimeter Institute for Theoretical Physics

PostPosted: June 08 2017  Reply with quote

I have tried rerunning CosmoMC ensuring I am using gfortran 6.2, and the issue has seemed to persist. For one run when it tries to restart from the checkpoint files, I get this:

Quote:

....

12 Reading checkpoint from chains/CDM22_12.chk
info = 0
—-
clik version 723c1a4b0580 MAKEFILE
bflike_smw
Checking likelihood './data/clik/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5cAp.clik' on test data. got −5247.87 expected −5247.87 (diff 5.67619e−07)
—-
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
11 Reading checkpoint from chains/CDM22_11.chk
info = 0
—-
clik version 723c1a4b0580 MAKEFILE
bflike_smw
Checking likelihood './data/clik/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5cAp.clik' on test data. got −5247.87 expected −5247.87 (diff 5.67619e−07)
—-
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
5 Reading checkpoint from chains/CDM22_5.chk
Chain:2 drag accpt: 0.566037714 fast/slow 140.622223 slow: 45
Chain:7 drag accpt: 0.555555582 fast/slow 138.679993 slow: 50
TObjectList: Unknown type to save


From another I get:

Quote:

....

reading BAO data set: DR11LOWZ
Doing non-linear Pk: F
Doing CMB lensing: T
Doing non-linear lensing: T
TT lmax = 2508
EE lmax = 2508
ET lmax = 2508
BB lmax = 2500
PP lmax = 2500
lmax_computed_cl = 2508
Computing tensors: F
max_eta_k = 14000.0000
transfer kmax = 5.00000000
adding parameters for: lowl_SMW_70_dx11d_2014_10_03_v5cAp
adding parameters for: smica_g30_ftl_full_pp
adding parameters for: BKPlanck_detset_comb_dust
adding parameters for: DR11CMASS
adding parameters for: DR11LOWZ
adding parameters for: MGS
adding parameters for: 6DF
adding parameters for: plik_dx11dr2_HM_v18_TTTEEE
Fast divided into 1 blocks
36 parameters (10 slow ( 0 semi-slow), 26 fast ( 0 semi-fast))
5 Reading checkpoint from chains/CDM23_5.chk
1 Reading checkpoint from chains/CDM23_1.chk
starting Monte-Carlo
Chain:7 drag accpt: 0.379746825 fast/slow 147.620255 slow: 79
TObjectList: Unknown type to save
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.