CosmoCoffee

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.

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+)