CosmoCoffee Forum Index CosmoCoffee

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

CosmoMC: does not generate .converge file
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Ira *Wolfson



Joined: 24 Jan 2013
Posts: 54
Affiliation: BGU

PostPosted: May 03 2017  Reply with quote

Hi,

For some reason my CosmoMC runs do not generate the .converge file with convergence statistics.

CosmoMC runs fine, and checking the chains with GetDist mid-run generates graphs that approximate the known results nicely.

However, I use the .converge file to asses how much time is left for the run, and plan workload ahead....


My .ini file:

Code:

#general settings
#Bicep-Keck-Planck, varying cosmological parameters
DEFAULT(batch2/BKPlanck.ini)

#Planck 2015, default just include native likelihoods (others require clik)
#DEFAULT(batch2/plik_dx11dr2_HM_v18_TT.ini)
#DEFAULT(batch2/lowTEB.ini)
#DEFAULT(batch2/lowl.ini)
DEFAULT(batch2/lensing.ini)

#Other Likelihoods
#DEFAULT(batch2/BAO.ini)
#DEFAULT(batch2/WiggleZ_MPK.ini)
#DEFAULT(batch2/MPK.ini)
#DEFAULT(batch2/WL.ini)

#general settings
DEFAULT(batch2/common.ini)

#e.g. to vary r in addition to standard 6:
#(for r>0 also need compute_tensors=T)
compute_tensors = T
param[r] = 0.03 0 2 0.04 0.04

redo_likelihoods = T
redo_theory =T
redo_cls = T
redo_pk = T
redo_skip = 0
redo_outroot =
redo_thin = 1
redo_add = F
redo_from_text = F

#high for new runs
MPI_Max_R_ProposeUpdate = 30

propose_matrix= planck_covmats/base_TT_lowTEB_plik.covmat

#Folder where files (chains, checkpoints, etc.) are stored
root_dir = chains/BKPNS/

#Root name for files produced
file_root=BKP_ns
#action= 0 runs chains, 1 importance samples, 2 minimizes
#use action=4 just to quickly test likelihoods
action = 0

num_threads = 8

#if you want to get theory cl for test point
#test_output_root = output_cl_root

start_at_bestfit =F
feedback=2
use_fast_slow = T

checkpoint = T

#sampling_method=7 is a new fast-slow scheme good for Planck
sampling_method = 7
dragging_steps  = 3
propose_scale = 2

#Set >0 to make data files for importance sampling
indep_sample=10

#these are just small speedups for testing
get_sigma8=T

#Uncomment this if you don't want one 0.06eV neutrino by default
#num_massive_neutrinos=
#to vary parameters set param[name]= center, min, max, start width, propose width
#param[mnu] = 0 0 0 0 0


#for PICO install from https://pypi.python.org/pypi/pypico and download data file
#cosmology_calculator=PICO
#pico_datafile = pico3_tailmonty_v34.dat
#pico_verbose=F

rand_seed =



Back to top
View user's profile [ Hidden ] Visit poster's website
Antony Lewis



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

PostPosted: May 03 2017  Reply with quote

CosmoMC generates a .converge_stat if run with MPI (after a while and burn in complete). .converge files are produced by getdist.
Back to top
View user's profile [ Hidden ] Visit poster's website
Ira *Wolfson



Joined: 24 Jan 2013
Posts: 54
Affiliation: BGU

PostPosted: May 03 2017  Reply with quote

I have
Code:
burn_in = 0
, running with MPI...

Still no .converge_stat file generated.

I am running Cosmomc from the master branch downloaded on April 30th.
Back to top
View user's profile [ Hidden ] Visit poster's website
Jason Dossett



Joined: 19 Mar 2010
Posts: 97
Affiliation: The University of Texas at Dallas

PostPosted: May 03 2017  Reply with quote

With burn_in = 0, when running in MPI CosmoMC calculates its own burn in period (every parameter you are varying has to have changed at least 51 times).

Check the outputs from stdout, you should see something like

Quote:

Chain 1 MPI done ''burn'', Samples = xxxx, like = xxxx
Time: xxxx output lines= xxxx


If you don't see that, your job has not completed burn in.

Another way to check and see if your chains have finished burn in is to check if you have any .chk files if checkpointing is enabled.
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.