[SciPy-user] odeint limitations?

Pearu Peterson pearu at cens.ioc.ee
Thu Mar 7 12:28:55 CST 2002


Hi,

I have used regularly vode for tasks with 256 upto 12288 coupled
nonlinear ordinary diff. equations. Scipy has interface to that routine
implemented in integrate.ode. The question of how long it would take does
not depend on the interface, it depends
1) on the size of the problem
2) on the properties of the system (stiff, non-stiff)
3) on whether you have choosen the right solver for the problem and tuned
   its parameters to optimal (again, depends on the system)
4) on how efficiently you can calculate the right-hand-side of the
   equations
5) on the hardware you have
etc.

And when the problem is of the size of hundred or thousands equations, it
is definitely _not_ recommended to calculate the jacobian in Python or
even in C or Fortran, _unless_, the jacobian is banded and  really simple.


More information about the SciPy-user mailing list