[Numpy-discussion] A step toward merging odeint and ode
Sun Aug 12 05:41:50 CDT 2012
I made a pull request  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!
More information about the NumPy-Discussion