Why is ay(1) going to infinity or NaN in derivs in CAMB

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
srijita sinha
Posts: 11
Joined: April 12 2019
Affiliation: Indian Institute of Science Education and Research Kolkata

Why is ay(1) going to infinity or NaN in derivs in CAMB

Post by srijita sinha » April 20 2020

HI,

I modified the equations_ppf.f90 for interacting DM-DE model. CAMB is working if I set the interaction = 0 in the modified code. But when interaction != 0 then a = ay(1) is becoming infinity or NaN, making others NaN too and the program stops. For non-interacting or LCDM a is in the range.
Can anyone tell me what is going wrong?

Many thanks in advance.

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

Re: Why is ay(1) going to infinity or NaN in derivs in CAMB

Post by Antony Lewis » April 20 2020

I expect you have a bug, compile and run with debug options to stop rather than generate NaNs.
(equations_ppf.f90 does not exist any more from camb version 1+)

Post Reply