CosmoCoffee Forum Index CosmoCoffee

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

minimally coupled dark energy solver introduced in camb for later use in cosmomc
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Nishanth Sasankan



Joined: 07 Sep 2016
Posts: 15
Affiliation: university of notre dame

PostPosted: November 21 2017  Reply with quote

I wanted to introduce a dynamic dark energy into CAMB, I am aware of the ppf module for camb, but i wanted to make a (w,a) array inside camb for future use with cosmomc, so I made a few modifications to camb.

Firstly I made a quintessence [V(phi)] solver which makes a table of a, wde and Grho(phi).
Then I set w not equal to −1 in params.ini (set it to w=−0.9)
Secondly i changed the subroutine DarkEnergy_ReadParams(Ini) to assign appf= a
wppf=wde
from the array generated by my solver
also nwmax and nw_ppf was changed to the length of the array made by the quintessence solver.

Thirdly subroutine wde(a) and grho_de(a) was changed to interpolate using the wde and Grho(phi) made by the solver.

after all this change I compiled camb and ran with params.ini

Now I have the following problem

Warning: Rombint failed to converge;
integral, error, tol: −3.0072019202949190E−003 1.6916759104756984E−005 1.0000000000000001E−005
Warning: Rombint failed to converge;
integral, error, tol: 7.8744507688045242E−004 1.2436078607158674E−004 1.0000000000000001E−005
Warning: Rombint failed to converge;
integral, error, tol: −1.2034253426512022E−004 4.5098700359258448E−004 1.0000000000000001E−005

after the first 2 warning lines, the code takes so much time before the 3rd warning comes up and then the following results show up

Reion redshift = 10.533
Om_b h2 = 0.022600
Om_c h2 = 0.112000
Om_nu h2 = 0.000640
Om_Lambda = 0.724000
Om_K = 0.000000
Om_m (1-Om_K-Om_L) = 0.276000
100 theta (CosmoMC) = 0.929762
Neff (total) = 3.046000
1 nu, g= 1.0153 mnu*c2/k_B/T_nu0= 353.71 (mnu= 0.060 eV)
Reion opt depth = 0.0900
Age of universe/GYr = 17.719
zstar = 1088.72
rs(zstar)/Mpc = 146.38
100*theta = 0.930036
DA(zstar)/Gpc = 15.73971
zdrag = 1059.70
rs(zdrag)/Mpc = 149.01
kD(zstar) Mpc = 0.1392
100*theta_D = 0.143348
zEQ (if vnu=1) = 3216.47
kEQ Mpc (if vnu=1) = 0.009817
100*theta_EQ = 0.758217
100*theta_rs_EQ = 0.417776
tau_recomb/Mpc = 284.76 tau_now/Mpc = 16024.4

Then it takes a very large amount of time to end.
Is there a mistake with the way I modified the code? or is this normal

Thanks
Nishanth
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.