[SciPy-user] odeint limitations?

[Travis Oliphant]

> > Hi,
> > I am a new user of SciPy. It looks like a very convenient method to access
> > powerful numerical codes written in other languages from Python.
> >
> > At this point, I have an application in which I need to solve many (hundreds
> > or thousands) of coupled ode's. Is the 'odeint' functionality in SciPy up to
> > this task? I believe the code uses lsod(e|a|ar|es), which should, in theory,
> > be adequate. Are there any 'issues' associated with the wrapping of these
> > routines in Python? Anyone have experience using odeint for solving large
> > numbers of ode's?
> >


I've not tried it for such large scale problems (I've only gone to 4 or 5
coupled ode's).

As you say, the underlying code should handle it.

The only question is how long will it take.  Using Numeric arrays in your
function to compute the derivative is definitely recommended, to save
space.

Why don't you give it a try?

-Travis Oliphant