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
other ODE integrators in CAMB
-
- Posts: 29
- Joined: June 10 2012
- Affiliation: University of Mexico
-
- Posts: 1944
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: other ODE integrators in CAMB
There's not built-in support for different integrators, but it should be relatively straightforward to replace dverk with another one.