[SciPy-user] SciPy-user Digest, Vol 65, Issue 50

Christopher MacMinn cmac at mit.edu
Thu Feb 19 16:03:21 EST 2009


>
> On Wed, Jan 21, 2009 at 1:22 PM, Rob Clewley <rob.clewley at gmail.com>
> wrote:
>
>> 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
> >
> > Ryan Gutenkunst already wrapped it while working on the SloppyCell
> package. See
> >
> > http://osdir.com/ml/python.scientific.devel/2005-07/msg00028.html
> >
> > with a link there to the code. I've never tried it myself or even
> > looked at it, FYI :)
> > -Rob
> >
>
> PS There's some mention of Ryan's lsodar.pyf in the trunk of scipy SVN, as
> per
>
>
> projects.scipy.org/scipy/scipy/browser/trunk/scipy/integrate/setup.py?rev=4763
>
> but I don't know if it's still there. If it is, is the associated pyd
> now shipped with Scipy? I haven't installed a new version for months.
> -Rob



Sorry for letting this drop for... err... awhile.  I don't speak
"developer", but I take this to mean that this functionality is available in
a few other python packages (PyDSTool, CVODE via PySundials), and that there
may or may not be pieces of lsodar already in scipy.

As a MATLAB to python/numpy/scipy convert who really still has one foot on
each log, as it were, having this 'root finding' functionality readily
available in python/scipy/numpy would be a big plus.  Consider this one vote
for finishing the wrapping of lsodar.  Of course, I should not say such
things without volunteering to help... I would be happy to contribute, but I
don't even know where to start.

Best, Chris
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20090219/815b769b/attachment.html>


More information about the SciPy-User mailing list