Page 1 of 1

Error in derived parameter numbers in CosmoMC-May2020

Posted: December 28 2020
by Saisandri Saini
Hello,
I am doing a quick test run with ModeChord where I am calculating the inflationary models with natural potential. I have successfully complied for cosmomc using make in source directory. I am getting the following run time error


Number of MPI processes: 1
file_root:natural
Random seeds: 22276, 3913 rand_inst: 1
read jla dataset data/Pantheon/full_long.dataset
reading WL data set: DES_1YR_final
Doing non-linear Pk: T
Doing CMB lensing: T
Doing non-linear lensing: T
TT lmax = 2500
EE lmax = 2500
ET lmax = 2500
BB lmax = 2500
PP lmax = 2500
lmax_computed_cl = 2500
Computing tensors: F
max_eta_k = 14000.0000
transfer kmax = 10.1999998
adding parameters for: JLA
adding parameters for: BK15_dust
adding parameters for: DES_1YR_final
Fast divided into 2 blocks
Block breaks at: 14
33 parameters ( 6 slow ( 0 semi-slow), 27 fast ( 0 semi-fast))
skipped unused params: calPlanck acib217 xi asz143 aps100 aps143 aps143217 aps217 aksz kgal100 kgal143 kgal143217 kgal217 cal0 cal2
starting Monte-Carlo
num_derived = 72 ; ix, Theory%numderived = 71 23
TP_CalcDerivedParams error in derived parameter numbers
MpiStop: 0
application called MPI_Abort(MPI_COMM_WORLD, 128) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=128
:
system msg for write_line failure : Bad file descriptor

I will be grateful if this could be resolved.

Regards,
Saisandri

Re: Error in derived parameter numbers in CosmoMC-May2020

Posted: December 29 2020
by Antony Lewis
Looks like an issue with new derived parameters? If this works OK with standard cosmomc it's probably an issue with ModeChord.

Re: Error in derived parameter numbers in CosmoMC-May2020

Posted: January 01 2021
by Saisandri Saini
Hi Antony,

Thank you for reply.
I have added new derived parameters in params_CMB.paramnames. Do I have to change other file too in paramnames for new derived parameters?

Re: Error in derived parameter numbers in CosmoMC-May2020

Posted: January 01 2021
by Antony Lewis
You need to define them in code as well of course, see e.g.
viewtopic.php?t=2128