[SciPy-user] ODE solvers and Scipy

Robert Clewley rclewley at cam.cornell.edu
Sat Aug 12 20:02:01 CDT 2006


On Sat, 12 Aug 2006, A. M. Archibald wrote:

> More seriously, if there are parts of the solution space that cause
> problems (the solution approaches a singularity, or the edge of a
> coordinate patch), there is no really effective way to bring the
> solution up to the edge and then stop. This is because lsoda doesn't
> do that - that feature is in lsodar.f (
> http://www.netlib.org/alliant/ode/prog/lsodar.f ), which has not
> (yet?) been wrapped in python.

I believe that lsodar *has* been wrapped as part of the sourceforge 
"SloppyCell" project for systems biology modelling. That package 
provides event detection for stopping integration under circumstances 
such as you mention, as does PyDSTool.

Rob


More information about the SciPy-user mailing list