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@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@scipy.org
https://mail.scipy.org/mailman/listinfo/scipy-dev