Aug. 13, 2006
1:02 a.m.
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