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.0072019202949190E003 1.6916759104756984E005 1.0000000000000001E005
Warning: Rombint failed to converge;
integral, error, tol: 7.8744507688045242E004 1.2436078607158674E004 1.0000000000000001E005
Warning: Rombint failed to converge;
integral, error, tol: 1.2034253426512022E004 4.5098700359258448E004 1.0000000000000001E005
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
minimally coupled dark energy solver introduced in camb for later use in cosmomc

 Posts: 27
 Joined: September 07 2016
 Affiliation: university of notre dame