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.
Why is ay(1) going to infinity or NaN in derivs in CAMB
-
- Posts: 32
- Joined: April 12 2019
- Affiliation: Indian Institute of Science Education and Research Kolkata
-
- Posts: 1941
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: Why is ay(1) going to infinity or NaN in derivs in CAMB
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+)
(equations_ppf.f90 does not exist any more from camb version 1+)