CosmoCoffee Forum Index CosmoCoffee

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

cosmomc: adding an extra parameter
Goto page Previous  1, 2  
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Antony Lewis



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

PostPosted: August 18 2016  Reply with quote

It is in the readme:

http://cosmologist.info/cosmomc/readme.html#FAQ
Back to top
View user's profile [ Hidden ] Visit poster's website
Suresh Kumar India



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

PostPosted: August 25 2016  Reply with quote

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:
 real(mcp) w, wa

I added
Code:
real(mcp) delta_coupled


(2) In CosmologyParameterizations.f90, I modified the code line
Code:
call this%SetTheoryParameterNumbers(16,last_power_index)

by
Code:
call this%SetTheoryParameterNumbers(17,last_power_index)


(3) In Calculator_CAMB.f90, below the code line
Code:
 real(dl) neff_massive_standard

I added
Code:
real(dl) delta_coupled


In the same file, below the code line
Code:
 wa_ppf = CMB%wa

I added
Code:
delta_coupled = CMB%delta_coupled


(4) In params_CMB.paramnames, below the code line
Code:
wa            w_a       #w_a variation

I added
Code:
delta_coupled delta_coupled


(5) In params_CMB_defaults.ini, I added
Code:
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:
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.
Back to top
View user's profile  
Antony Lewis



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

PostPosted: August 25 2016  Reply with quote

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.
Back to top
View user's profile [ Hidden ] Visit poster's website
Suresh Kumar India



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

PostPosted: August 25 2016  Reply with quote

It is working now. Thanks.
Back to top
View user's profile  
Suresh Kumar India



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

PostPosted: August 26 2016  Reply with quote

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?
Back to top
View user's profile  
Antony Lewis



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

PostPosted: August 27 2016  Reply with quote

I don't think it's interfaced to CAMB.
Back to top
View user's profile [ Hidden ] Visit poster's website
R Purnalingam



Joined: 03 Mar 2017
Posts: 19
Affiliation: IIT KANPUR

PostPosted: April 05 2017  Reply with quote

I am also getting the same error .Please explain how to resolve it?
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
Goto page Previous  1, 2
Page 2 of 2

 
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.