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

Fabrice Silva silva at lma.cnrs-mrs.fr
Wed Aug 15 04:42:18 EDT 2012

Le mardi 14 août 2012 à 21:21 +0200, Ralf Gommers a écrit :
> On Sun, Aug 12, 2012, Fabrice Silva <silva at lma.cnrs-mrs.fr> wrote:
>         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.
> Does that mean that odeint can be made a wrapper around lsoda and that
> the odepack static extension can be completely removed?

Hi Ralf,
The pull request allows to run the integration using the object-oriented
interface ode, with the same solver than the odeint interface uses, i.e.
lsoda, extending the integrators available for the object-oriented

As I understand the scipy.integrate architecture, we are by now
* the odepack library, which has all the fortran sources required by
lsoda and vode at least.
* the _odepack extension, which defines the _odepack module needed by

This latter would be removable, and odeint a wrapper around the lsoda
pyf'ed function. I suppose you are talking about the _odepack extension,
am I wrong?

Fabrice Silva

More information about the NumPy-Discussion mailing list