cosmomc: adding an extra parameter

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Antony Lewis
Posts: 1354
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: cosmomc: adding an extra parameter

Post by Antony Lewis » August 18 2016


Suresh Kumar India
Posts: 14
Joined: May 29 2015
Affiliation: BITS-Pilani, India

cosmomc: adding an extra parameter

Post by Suresh Kumar India » August 25 2016

Thanks. But still it is not working. Following the read me, I did the following to add a new parameter namely "delta_coupled".

(1) In CosmologyTypes.f90, below the code line

Code: Select allp>

 real(mcp) w, wa 
I added

Code: Select allp>

real(mcp) delta_coupled
(2) In CosmologyParameterizations.f90, I modified the code line

Code: Select allp>

call this%SetTheoryParameterNumbers(16,last_power_index)
by

Code: Select allp>

call this%SetTheoryParameterNumbers(17,last_power_index)
(3) In Calculator_CAMB.f90, below the code line

Code: Select allp>

 real(dl) neff_massive_standard 
I added

Code: Select allp>

real(dl) delta_coupled
In the same file, below the code line

Code: Select allp>

 wa_ppf = CMB%wa 
I added

Code: Select allp>

delta_coupled = CMB%delta_coupled
(4) In params_CMB.paramnames, below the code line

Code: Select allp>

wa            wa       #wa variation
I added

Code: Select allp>

delta_coupled delta_coupled
(5) In params_CMB_defaults.ini, I added

Code: Select allp>

param[delta_coupled] = 0 −1 1 0.01 0.01
After the above changes, when I run the code, it stops with the following error:

Code: Select allp>

Fast divided into            1  blocks
 25 parameters (11 slow ( 0 semi-slow), 14 fast ( 0 semi-fast))
 starting Monte-Carlo
 Compile with CosmoRec to use fdm
 MpiStop:            0
Please reply in detail so that I can run the code correctly with the new parameter.

Antony Lewis
Posts: 1354
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: cosmomc: adding an extra parameter

Post by Antony Lewis » August 25 2016

I think you need to get/set the value of your new CMB%delta_coupled parameter in CosmologyParameterizations.f90, updating indexing of other params if necessary.

Suresh Kumar India
Posts: 14
Joined: May 29 2015
Affiliation: BITS-Pilani, India

cosmomc: adding an extra parameter

Post by Suresh Kumar India » August 25 2016

It is working now. Thanks.

Suresh Kumar India
Posts: 14
Joined: May 29 2015
Affiliation: BITS-Pilani, India

cosmomc: adding an extra parameter

Post by Suresh Kumar India » August 26 2016

I see that the parameter cs2_lam is used in equations.f90, but it is not listed in params_CMB_defaults.ini and params_CMB.paramnames just like wlam. Is it hard coded or listed somewhere else?

Antony Lewis
Posts: 1354
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: cosmomc: adding an extra parameter

Post by Antony Lewis » August 27 2016

I don't think it's interfaced to CAMB.

R Purnalingam
Posts: 19
Joined: March 03 2017
Affiliation: IIT KANPUR
Contact:

cosmomc: adding an extra parameter

Post by R Purnalingam » April 05 2017

I am also getting the same error .Please explain how to resolve it?

Post Reply