Run CosmoMC without analytical Hubble parameter

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
zahra safari
Posts: 1
Joined: October 27 2018
Affiliation: Univercity of Kurdistan

Run CosmoMC without analytical Hubble parameter

Post by zahra safari » December 30 2018

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

Antony Lewis
Posts: 1457
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: Run CosmoMC without analytical Hubble parameter

Post by Antony Lewis » January 05 2019

You need to at least also update the derivs() function (for perturbation evolution) in equations.f90,

Post Reply