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