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

Evgeni Burovski evgeny.burovskiy at gmail.com
Mon Nov 23 03:28:45 EST 2015


Hi David,

This sounds like a useful feature.
I guess a pull request would be welcome. (even though I'm not a user, and I
won't be able to review it).

Evgeni
23.11.2015 10:47 пользователь "David Mikolas" <david.mikolas1 at gmail.com>
написал:

> 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
>>
>
>
> _______________________________________________
> 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/b00a68b0/attachment.html>


More information about the SciPy-Dev mailing list