[SciPy-Dev] dense output option in integrate.ode

[David Mikolas]

One way that dense output for dopri5 and dop853 can be very useful is if
the integration is expensive/long and you don't want to repeat it, but you
want to obtain results at new time points at a later date, or even iterate
on it - for example, find the time of closest approac. This is done by
saving the interpolation coefficients. I put a simple example here, though
there is no saving to disk yet.

    http://pastebin.com/e6qNjbL9   dendop_test_v00.py

I wonder if this could be developed into an option to return an
interpolator object or function. If I can help let me know.

On Sat, Nov 14, 2015 at 12:19 AM, Evgeni Burovski <
evgeny.burovskiy at gmail.com> wrote:

> Hi David, Hi Jim,
>
> > I am new to SciPy-dev. Will the dense output option in
> scipy.integrate.ode
> > become available in 0.17.0? This is a feature already available in the
> > original FORTAN, but wasn't implemented in the wrapper.
>
> If the feature is sent as a pull request against the scipy master
> branch, the PR is reviewed by the maintainers of the integrate package
> and merged into master before the release split, then yes, it would be
> available in 0.17.0.
>
> So far I do not see any progress towards it.
>
>
> > I wrote these dense output extensions that you listed:
> >
> >
> https://github.com/jddmartin/scipy/tree/dense_output_from_dopri5_and_dop853
> > https://github.com/jddmartin/dense_output_example_usage
> >
> > but I didn't issue any pull request to scipy.  I sent this message to the
> > scipy developers list:
> >   http://article.gmane.org/gmane.comp.python.scientific.devel/19635/
> > explaining the changes.  I was hoping for some feedback before issuing a
> > pull request.
>
> Ah, I see that the email likely fell through the cracks back in April.
> Sorry about that.
> You might want to ping that email thread once more or send a pull
> request on github (or both).
>
> <snip>
>
> Cheers,
>
> Evgeni
> _______________________________________________
> SciPy-Dev mailing list
> SciPy-Dev at scipy.org
> https://mail.scipy.org/mailman/listinfo/scipy-dev
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20151123/4a717522/attachment.html>