Nishanth Sasankan
Joined: 07 Sep 2016 Posts: 26 Affiliation: university of notre dame

Posted: November 21 2017 


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, w_{de} 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 a_{ppf}= a
w_{ppf}=w_{de}
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 w_{de}(a) and grho_de(a) was changed to interpolate using the w_{de} 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 h^{2} = 0.022600
Om_c h^{2} = 0.112000
Om_nu h^{2} = 0.000640
Om_Lambda = 0.724000
Om_K = 0.000000
Om_m (1Om_KOm_L) = 0.276000
100 theta (CosmoMC) = 0.929762
N_{eff} (total) = 3.046000
1 nu, g= 1.0153 m_{nu}*c^{2}/k_B/T_nu0= 353.71 (m_{nu}= 0.060 eV)
Reion opt depth = 0.0900
Age of universe/GYr = 17.719
zstar = 1088.72
r_{s}(zstar)/Mpc = 146.38
100*theta = 0.930036
DA(zstar)/Gpc = 15.73971
zdrag = 1059.70
r_{s}(zdrag)/Mpc = 149.01
k_{D}(zstar) Mpc = 0.1392
100*theta_D = 0.143348
z_{EQ} (if v_{nu}=1) = 3216.47
k_{EQ} Mpc (if v_{nu}=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 
