CosmoCoffee Forum Index CosmoCoffee

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

axionCAMB and cosmomc: how to transfer the CAMB params to CMB params
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
chen lu



Joined: 23 Jul 2016
Posts: 4
Affiliation: ITP

PostPosted: April 06 2017  Reply with quote

Hi,
I tried to link cosmomc with axionCAMB, a modified code of the old camb in cosmomc which adds something about ultra-light axions.

I have added the two new initial params(max and omaxh2) to cosmomc successfully. But it stopped when the file named *.out showed as follows :
...
max_eta_k = 14000.00
transfer kmax = 5.000000
adding parameters for: lowl_SMW_70_dx11d_2014_10_03_v5cAp
adding parameters for: plik_dx11dr2_HM_v18_TT
Fast divided into 1 blocks
23 parameters ( 9 slow ( 0 semi-slow), 14 fast ( 0 semi-fast))

Maybe it occurred because some params about axions are not included in cosmomc. Then I want to transfer them from axionCAMB to CMB because they have been defined and calculated in axionCAMB. So can you guys tell me where I should modify the code in cosmomc/source to do it?

Tks!
Back to top
View user's profile  
Luke Hart



Joined: 13 Jul 2015
Posts: 24
Affiliation: University of Manchester

PostPosted: April 07 2017  Reply with quote

Have a look at Calculator_CAMB.f90 in the source folder. I did all my mods in the Calc_CMBtoCAMB bit I think but it might be different for axion stuff I'm not sure? But you just have to find out where the parameters you want are in Camb


(usually denoted by P%ParamType%newparam_CAMB = CMB%new_param)

I've recently done this so if you get stuck, give me another shout! Hope this helps :)
Back to top
View user's profile  
chen lu



Joined: 23 Jul 2016
Posts: 4
Affiliation: ITP

PostPosted: April 11 2017  Reply with quote

Thank you for your help.
I have added the new initial params (P%new_params = CMB%new_params) in the subroutine CAMBCalc_CMBToCAMB. But it did not work. Maybe there are other new derived params which are calculated in axionCAMB and I have added the derived params' definitions in cosmomc, but I don't know how to transfer them from axionCAMB to cosmomc. Do you have any ideas about this?
Thank you very much.
Back to top
View user's profile  
Luke Hart



Joined: 13 Jul 2015
Posts: 24
Affiliation: University of Manchester

PostPosted: April 11 2017  Reply with quote

Have you made sure that your parameter has also been added to the .paramnames file in the cosmomc directory? I had a problem with that when I set it up.

Chen if you have a look at EFTCAMB and compare to normal CAMB, there is also a version of CosmoMC called EFTCosmoMC you can compare and see how they added in their effective field parameters?

Luke
Back to top
View user's profile  
Luke Hart



Joined: 13 Jul 2015
Posts: 24
Affiliation: University of Manchester

PostPosted: April 11 2017  Reply with quote

It definitely looks like (given the number of slow parameters) that CosmoMC is reading in your parameters. You might want to try and see if theres another auxillary file (.o or .e files) that will give you a better idea of the error? If not, then go with the EFTCAMB method I mentioned before.
Back to top
View user's profile  
chen lu



Joined: 23 Jul 2016
Posts: 4
Affiliation: ITP

PostPosted: April 13 2017  Reply with quote

Thank you very much!

I'm sure I have added the two new initial parameters successfully because I can check them in chains/test.ranges.

Maybe there's something else wrong and I will check things more than the params.

Tks again for your reply.
Back to top
View user's profile  
Display posts from previous:   
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software All times are GMT + 5 Hours
Page 1 of 1

 
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.