Page 1 of 1

Run CosmoMC without analytical Hubble parameter

Posted: December 30 2018
by zahra safari
hi Antony

I want to constraint a dark energy model with CosmoMC. Model has two coupled differential equation for evaluate the Hubble parameter. I modified the function "dtauda" in equations_ppf.f90, and add Runge Kutta subroutine to subtoutins.f90 for numerical solving this equations. After compile and run code, I see that number of iterations very very largest than my code step size. Step size is 10 and iteration is 911520!!! Please help me to debug the code.
Thank you.
zahra

Re: Run CosmoMC without analytical Hubble parameter

Posted: January 05 2019
by Antony Lewis
You need to at least also update the derivs() function (for perturbation evolution) in equations.f90,