CosmoCoffee

Posted: **May 03 2017**

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: Select all

#general settings
#Bicep-Keck-Planck, varying cosmological parameters
DEFAULT&#40;batch2/BKPlanck.ini&#41;

#Planck 2015, default just include native likelihoods &#40;others require clik&#41;
#DEFAULT&#40;batch2/plik_dx11dr2_HM_v18_TT.ini&#41;
#DEFAULT&#40;batch2/lowTEB.ini&#41;
#DEFAULT&#40;batch2/lowl.ini&#41;
DEFAULT&#40;batch2/lensing.ini&#41;

#Other Likelihoods
#DEFAULT&#40;batch2/BAO.ini&#41;
#DEFAULT&#40;batch2/WiggleZ_MPK.ini&#41;
#DEFAULT&#40;batch2/MPK.ini&#41;
#DEFAULT&#40;batch2/WL.ini&#41;

#general settings
DEFAULT&#40;batch2/common.ini&#41;

#e.g. to vary r in addition to standard 6&#58;
#&#40;for r>0 also need compute_tensors=T&#41;
compute_tensors = T
param&#91;r&#93; = 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 &#40;chains, checkpoints, etc.&#41; 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&#91;name&#93;= center, min, max, start width, propose width
#param&#91;mnu&#93; = 0 0 0 0 0


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

rand_seed =

Posted: **May 03 2017**

CosmoMC generates a .converge_stat if run with MPI (after a while and burn in complete). .converge files are produced by getdist.

Posted: **May 03 2017**

I have

Code: Select all

burn_in = 0

, running with MPI...

Still no .converge_stat file generated.

I am running Cosmomc from the master branch downloaded on April 30th.

Posted: **May 03 2017**

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

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.

CosmoCoffee

CosmoMC: does not generate .converge file

CosmoMC: does not generate .converge file

Re: CosmoMC: does not generate .converge file

CosmoMC: does not generate .converge file

CosmoMC: does not generate .converge file