[SciPy-user] ODE solvers and Scipy
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.
More information about the SciPy-user