[SciPy-Dev] Dormand-Prince 8(5,3) in new SciPy IVP integrators

David Hagen david at drhagen.com
Wed Sep 19 07:20:04 EDT 2018


I agree that this would be a good solver to add. Please mention @drhagen on
the PR also. For testing, the solver should be added to all the testing
sections that make sense. But we can discuss the details of that on the PR.

On Tue, Sep 18, 2018 at 2:37 AM Nikolay Mayorov <nikolay.mayorov at zoho.com>
wrote:

> Hi, Juan!
>
> I know that several people were interested to see this method implemented.
> I know it has a difficult logic for step size control so I couldn't
> implement it quickly. Now that you have a working version (say it might
> contain some small errors, but still the core logic is here)  conforming
> the new API I think you should definitely put the pull request. Perhaps
> adding a new simple example from celestial mechanics (say) to demonstrate
> how it works will be useful and interesting too (say compare lower order RK
> methods with it).
>
> Please create the PR and mention @nmayorov here.
>
> Best,
> Nikolay
>
>
>
> ---- On Mon, 17 Sep 2018 15:00:34 +0500 *Juan Luis Cano
> <juanlu001 at gmail.com <juanlu001 at gmail.com>>* wrote ----
>
> Hi all,
>
> We wanted to use the new initial value problem solvers present in SciPy
> but the DOP853 method was not present there, so we implemented it ourselves:
>
>
> https://github.com/poliastro/poliastro/blob/v0.11b1/src/poliastro/integrators.py
>
> We have done some soft validation (i.e. it has worked well so far for our
> specific use case) and we think it would be great to have it in SciPy,
> given that it was present in the "Old API" and that it's a widely used
> method in Celestial Mechanics.
>
> If I understood correctly, in-depth quantitative validation was not
> written for the current solvers:
>
> https://github.com/scipy/scipy/pull/6326#issuecomment-250830082
>
> Would you be interested in a PR? For testing, would it be enough to add
> the new method here?
>
>
> https://github.com/scipy/scipy/blob/4d16aea1267898e8532f962f0104d485578ab670/scipy/integrate/tests/test_ivp.py#L147
>
> Regards,
>
> --
> Juan Luis Cano
> _______________________________________________
> SciPy-Dev mailing list
> SciPy-Dev at python.org
> https://mail.python.org/mailman/listinfo/scipy-dev
>
>
>
> _______________________________________________
> SciPy-Dev mailing list
> SciPy-Dev at python.org
> https://mail.python.org/mailman/listinfo/scipy-dev
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20180919/41454a67/attachment.html>


More information about the SciPy-Dev mailing list