[SciPy-user] odeint limitations?
oliphant at ee.byu.edu
Thu Mar 7 09:02:44 CST 2002
> > 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
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
Why don't you give it a try?
More information about the SciPy-user