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

[Nikolay Mayorov]

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> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20180918/25cf9a73/attachment-0001.html>