[SciPy-user] ode45 replacment

A. M. Archibald peridot.faceted at gmail.com
Sun Nov 5 03:48:42 CST 2006


On 04/11/06, Steve Schmerler <elcorto at gmx.net> wrote:
> A. M. Archibald wrote:
> > On 03/11/06, Ryan Krauss <ryanlists at gmail.com> wrote:
> >> What is the best scipy replacement for Matlab's ode45?  Is it
> >> scipy.integrate.odeint?  What is the difference between odeint and
> >> ode?
> >
> > Either one will probably suit your needs. For practical purposes, the
> > difference is one of interface - odeint has a more procedural
> > interface, while ode has a more object-oriented one. Underneath, both
> > use well-established fortran ODE integrators (odeint uses LSODA from
> > ODEPACK, and ode uses VODE). odeint may also have problems on a 64-bit
> > system, IIRC.
> >
>
> I'm planning to set up a 64bit machine and use odeint.
> In case you remember any thread on the list or can point me in the right
> direction, I'd be happy to hear about it. A quick search of the list
> revealed nothing specific...

I had some difficulties with a FORTRAN wrapper that assumed machine
integers were 32 bits, but the current version of odeint appears to be
fine. If that problem does occur, it is easily fixed with a bit of
grepping through the python interface file.

A. M. Archibald


More information about the SciPy-user mailing list