[SciPy-user] integrate.odeint and event handling

Warren Weckesser warren.weckesser@enthought....
Wed Jan 21 12:18:07 CST 2009


Christopher W. MacMinn wrote:
> Hello -
>
> I was wondering if integrate.odeint offers any event handling  
> capabilities.
>
>
>   
<snip>
> Thanks!
>
> Best, Chris
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>
>   

Hi Chris,

As Rob Clewley pointed out, odeint does not provide event detection.  I
don't think SciPy's ode class does, either.  Rob's PyDSTool is one
alternative (and it provides a lot of other nice tools to go along with
the ODE solver); another is PySUNDIALS, as mentioned in the thread to
which Rob provided a link.

odeint is a wrapper for the LSODA solver in the Fortran ODEPACK
library.  This library also includes LSODAR, which is LSODA with
root-finding (aka event detection).  Does anyone want to take a stab at
wrapping LSODAR?  The wrapping of LSODA with odeint provides a good
starting point, and an ODE solver with root-finding would be a great
addition to SciPy.

Warren


-- 
Warren Weckesser
Enthought, Inc.
515 Congress Avenue, Suite 2100
Austin, TX  78701
512-536-1057 x249



More information about the SciPy-user mailing list