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

[Juan Luis Cano]

For the rest of the subscribers, the pull request is in:

https://github.com/scipy/scipy/pull/9290

On Wed, Sep 19, 2018 at 1:44 PM David Hagen <david at drhagen.com> wrote:

> 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
>>
> _______________________________________________
> SciPy-Dev mailing list
> SciPy-Dev at python.org
> https://mail.python.org/mailman/listinfo/scipy-dev
>


-- 
Juan Luis Cano
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20180920/435c0d56/attachment.html>