Re: ENH: add cumulative_simpson integration to scipy.integrate

On Mon, Mar 20, 2023 at 2:05 PM Prav <[email protected]> wrote:

> Hello everyone!
>
> The cumulative integration and differentiation of numerical samples is a
> common operation in the engineering domains. Scipy currently only provides
> the integrate.cumulative_trapezoid function for cumulative integration.
> However, this is insufficiently accurate for engineering purposes.
> The integrate.simpson method is available but only as a complete integral,
> not a cumulative integral.
>
> I propose to add a cumulative integration function using Simpson's 1/3
> rule (https://en.wikipedia.org/wiki/Simpson%27s_rule).
> CAE software such as ANSYS/LS-DYNA provide simpson's integration as a
> default processing option. Adding this to Scipy.integrate will help expand
> the utility of the library in the engineering domains.
>
> I have added a PR for this here: https://github.com/scipy/scipy/pull/18151
>
> I would really appreciate your reviews and opinions.
>

Hi Praveer, thanks for proposing this feature and sharing your PR and
notebook with examples. A more accurate cumulative integration method seems
like a very reasonable addition to `scipy.integrate` to me. And your PR
looks pretty polished already. So +1 from me.

The replies/review here may be a bit slow because we're really short on
folks with algorithmic expertise in this area - our resident experts are in
low-activity mode on SciPy right now. Please do ping again if it takes too
long.

Cheers,
Ralf

_______________________________________________
SciPy-Dev mailing list -- [email protected]
To unsubscribe send an email to [email protected]
https://mail.python.org/mailman3/lists/scipy-dev.python.org/
Member address: [email protected]