other ODE integrators in CAMB
Posted: December 06 2018
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
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