Planck 2018 data + SNe in CosmoMC-Planck2018 branch

[Erick Almaraz]

Dear all,

I have downloaded the CosmoMC package (planck2018 branch) and successfully installed in my lap (using ifort)

Code: Select all

ealmaraz@ealmaraz:~$ ifort --version
ifort (IFORT) 15.0.7 20160518
Copyright (C) 1985-2016 Intel Corporation.  All rights reserved.

Without modifying the code, the program runs well with Planck 2018 data, eg,:

Code: Select all

DEFAULT(batch3/lowl.ini)
DEFAULT(batch3/plik_rd12_HM_v22_TT.ini)

Also, the code runs well when I use jointly Planck 2018 data and BAO information:

Code: Select all

DEFAULT(batch3/lowl.ini)
DEFAULT(batch3/plik_rd12_HM_v22_TT.ini)
DEFAULT(batch3/BAO.ini)

However, when I tried to incorporate SNe information (JLA or Pantheon) the program simply crashes. Here is the output with feedback = 2:

Code: Select all

ealmaraz@ealmaraz:~/software/cosmomc/lcdm/v2.2$ mpirun -np 1 ./cosmomc test.ini 
 Number of MPI processes:           1
 file_root:test
 Random seeds: 29972,  5528 rand_inst:   1
 Using clik with likelihood file 
 ./data/clik_14.0/low_l/commander/commander_dx12_v3_2_29.clik
 Initialising clik...
----
clik version 3be036bbb4f9
  gibbs_gauss b13c8fda-1837-41b5-ae2d-78d6b723fcf1
Checking likelihood './data/clik_14.0/low_l/commander/commander_dx12_v3_2_29.clik' on test data. got -11.6257 expected -11.6257 (diff -1.07424e-09)
----
 Likelihood uses the following Cls:
   TT from l=0 to l=          29
 Clik will run with the following nuisance parameters:
 A_planck
 Using clik with likelihood file 
 ./data/clik_14.0/hi_l/plik/plik_rd12_HM_v22_TT.clik
 Initialising clik...
----
clik version 3be036bbb4f9
  smica
Checking likelihood './data/clik_14.0/hi_l/plik/plik_rd12_HM_v22_TT.clik' on test data. got -380.341 expected -380.341 (diff -1.08197e-08)
----
 Likelihood uses the following Cls:
   TT from l=0 to l=        2508
 Clik will run with the following nuisance parameters:
 A_cib_217
 cib_index
 xi_sz_cib
 A_sz
 ps_A_100_100
 ps_A_143_143
 ps_A_143_217
 ps_A_217_217
 ksz_norm
 gal545_A_100
 gal545_A_143
 gal545_A_143_217
 gal545_A_217
 A_sbpx_100_100_TT
 A_sbpx_143_143_TT
 A_sbpx_143_217_TT
 A_sbpx_217_217_TT
 calib_100T
 calib_217T
 A_planck
 JLA dataset name: PANTHEON
 JLA data file: data/Pantheon/lcparam_full_long_zhel.txt
 Number of SN read: 1048
 read jla dataset data/Pantheon/full_long.dataset
 Doing non-linear Pk: F
 Doing CMB lensing: T
 Doing non-linear lensing: F
 TT lmax =  2508
 EE lmax =  2500
 ET lmax =  2500
 lmax_computed_cl  =  2508
 Computing tensors: F
 max_eta_k         =    6645.000    
 transfer kmax     =    1.000000    
 adding parameters for: JLA
 JLA data param indices:          25          26
 adding parameters for: commander_dx12_v3_2_29
 commander_dx12_v3_2_29 data param indices:          27
 adding parameters for: plik_rd12_HM_v22_TT
 plik_rd12_HM_v22_TT data param indices:          28          29          30
          31          32          33          34          35          36
          37          38          39          40          41          42
          43          44          45          46          27
 Fast divided into            2  blocks
 Block breaks at:            8
 21 parameters ( 6 slow ( 2 semi-slow), 15 fast ( 0 semi-fast))
 starting Monte-Carlo
           1 Likelihood:   1.0000000E+30 Current Like:  1.0000000E+30
 Initialising BBN data PArthENoPE_880.2_standard.dat, col:           4
 Done. Interpolation table is           48  by           11
 Calling CAMB
 CAMB done
-------------------------------------------------------
Primary job  terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

  Process name: [[63218,1],0]
  Exit code:    174
--------------------------------------------------------------------------

Any idea about what is causing the problem when incorporating SNe data?

Regards

[Antony Lewis]

Not sure, try building and running cosmomc_debug

[Erick Almaraz]

Dear Antony,

It seems that the problem is in my lap, since I installed CosmoMC in a different machine and now everything goes well.

Regards