[Numpy-discussion] A step toward merging odeint and ode

Fabrice Silva silva@lma.cnrs-mrs...
Sun Aug 12 05:41:50 CDT 2012

I made a pull request [1] to integrate the LSODA solver that is used in
odeint into the modular scipy.integrate.ode generic class. In a similar
way as for vode, it just wraps the already present lsoda.f file
(see .pyf file) and exposes it within an IntegratorBase subclass
adjusting the coefficients before calling lsoda.

Note that lsoda provide automatic switching between stiff and non-stiff
methods, feature that is not present in the available vode integrator.

Final note: tests are ok!


[1] https://github.com/scipy/scipy/pull/273 
Fabrice Silva

More information about the NumPy-Discussion mailing list