minimally coupled dark energy solver introduced in camb for later use in cosmomc

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
Nishanth Sasankan
Posts: 27
Joined: September 07 2016
Affiliation: university of notre dame

minimally coupled dark energy solver introduced in camb for

Post by Nishanth Sasankan » 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 (1-Om_K-Om_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

Post Reply