CosmoCoffee

Dear Cosmo Coffee forum,

I am implementing an alternative DE model in CAMB. However, the dynamical background equations seem to be stiff and the default integrator dverk available in $CAMB/subroutines.f90 is not able to solve these equations. On the other hand, these same equations can be solved numerically in python & mathematica. I have heard that there is a special library (ODEPACK) written in fortran 77 with which is possible to solve stiff equations. Is it possible to plug this library into CAMB, so that we can use the subroutines therein to solve stiff ODES?


Regards

There's not built-in support for different integrators, but it should be relatively straightforward to replace dverk with another one.