Repository navigation
Implementation of scipy.integrate.simps #2629
Description
Activity
Is this something that you're interested in implementing @Joshuaalbert ?
- Yes, I would like to. Is there a developers guide I should read first?…On 2017-08-26 08:47, Matthew Rocklin wrote: Is this something that you're interested in implementing @Joshuaalbert [1] ? -- You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub [2], or mute the thread [3]. * Links: ------ [1] https://github.com/joshuaalbert [2] #2629 (comment) [3] https://github.com/notifications/unsubscribe-auth/AOHv-LN7biOizmP-tkrGUsMWfmGNLr38ks5scBP4gaJpZM4PDTAL
I apologize for the long delay here. http://dask.pydata.org/en/latest/develop.html
This issue has been open for a while, contributions are welcome from anyone who'd like to tackle it.
@GenevieveBuckley/ @mrocklin may I work on this please?
Also maybe we should create adask-scipysort of repo?@mrinalsardar @mrocklin @GenevieveBuckley What is the status of this issue? Can I pick it up?
Based on the content of the issue, it doesn't looks like there has been any progress. Feel free to work on this if you'd like @wkrasnicki, one doesn't need to ask permission to pick up an issue
Hi everyone, first-time contributor here! I’m tackling this issue right now and will open a draft PR soon. Any feedback is appreciated!
Reacted by Wiktor Krokosz- added 6 commits that reference this issue
on Jun 14, 2025 I'd love to implement this!
Note:
scipy.integrate.simpshas been renamed toscipy.integrate.simpsonin newer scipy versions. I'll implementdask.array.simpsonthat wraps the scipy implementation to work on chunked arrays via the "overlapping" (overlap/ghost) pattern.Approach:
- Use
dask.array.lib.stride_tricksand the overlap pattern to handle chunk boundaries correctly - Apply Simpson's rule across each chunk with appropriate boundary handling
- Support the
axisanddx/xparameters
Would this be best added to
dask/array/routines.py? Happy to start a draft PR!- Use
Are there plans to implement
scipy.integrate.simps?It is a pretty useful method that would easily gain from chunking (when integration is done down one axis of a multi-N array).