Improvement to regular grid interpolation
Dear scipy developers, One could use *scipy.interpolate.RegularGridInterpolator* at different times with different values but the same source grid and target points. Currently this would recompute the indices each time the function is called. I would like to add the possibility to provide the target points at initialisation and the values when calling the object. Best regards, Edouard Goudenhoofdt
Edouard, I think this thread touched on that a bit: https://mail.python.org/pipermail/scipy-dev/2019-May/023540.html So Simon Clift might have already some functional code for that which could serve as a starting point. I personally would be certainly interested in the functionality. Dieter On 18/09/2019 11:50, Edouard Goudenhoofdt wrote:
Dear scipy developers,
One could use /scipy.interpolate.RegularGridInterpolator/ at different times with different values but the same source grid and target points. Currently this would recompute the indices each time the function is called. I would like to add the possibility to provide the target points at initialisation and the values when calling the object.
Best regards,
Edouard Goudenhoofdt //
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
This is clearly a duplicate of your previous post (How can I search the mailing list?). I think this is a basic task for many geoscientists. I did not make the effort to understand Simon Clift smarter approach. But I already implemented a quick solution using a switch for the two use cases. May I submit a pull request? On Wed, Sep 18, 2019 at 12:47 PM Dieter Werthmüller <Dieter@werthmuller.org> wrote:
Edouard,
I think this thread touched on that a bit: https://mail.python.org/pipermail/scipy-dev/2019-May/023540.html
So Simon Clift might have already some functional code for that which could serve as a starting point.
I personally would be certainly interested in the functionality.
Dieter
On 18/09/2019 11:50, Edouard Goudenhoofdt wrote:
Dear scipy developers,
One could use /scipy.interpolate.RegularGridInterpolator/ at different times with different values but the same source grid and target points. Currently this would recompute the indices each time the function is called. I would like to add the possibility to provide the target points at initialisation and the values when calling the object.
Best regards,
Edouard Goudenhoofdt //
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
On Wed, Sep 18, 2019 at 2:33 PM Edouard Goudenhoofdt <egouden@gmail.com> wrote:
This is clearly a duplicate of your previous post (How can I search the mailing list?).
Yeah, good question .... that's kind of problematic at the moment. Eventually we'll switch to Mailman 3 which has a sane search interface. Sorry, we don't have a better answer than using Google et al. right now. I think this is a basic task for many geoscientists.
I did not make the effort to understand Simon Clift smarter approach. But I already implemented a quick solution using a switch for the two use cases. May I submit a pull request?
It's now come up a couple of times, so clearly there's interest. Yes, opening a PR would be very helpful. Cheers, Ralf
On Wed, Sep 18, 2019 at 12:47 PM Dieter Werthmüller < Dieter@werthmuller.org> wrote:
Edouard,
I think this thread touched on that a bit: https://mail.python.org/pipermail/scipy-dev/2019-May/023540.html
So Simon Clift might have already some functional code for that which could serve as a starting point.
I personally would be certainly interested in the functionality.
Dieter
On 18/09/2019 11:50, Edouard Goudenhoofdt wrote:
Dear scipy developers,
One could use /scipy.interpolate.RegularGridInterpolator/ at different times with different values but the same source grid and target points. Currently this would recompute the indices each time the function is called. I would like to add the possibility to provide the target points at initialisation and the values when calling the object.
Best regards,
Edouard Goudenhoofdt //
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
Hi Edouard, I'll keep an eye out for your pull request, or perhaps you could send me a note on your use case to my gmail address. My approach just consists of saving the interpolation weights in a sparse matrix then doing a sparse multiply for subsequent interpolations. Weights can be either multi-linear or using barycentric co-ordinates on regular grids; it's an arbitrary choice unless you have a simplicial grid layout you prefer. Cheers -- Simon On Wed, Sep 18, 2019 at 8:50 AM Ralf Gommers <ralf.gommers@gmail.com> wrote:
On Wed, Sep 18, 2019 at 2:33 PM Edouard Goudenhoofdt <egouden@gmail.com> wrote:
This is clearly a duplicate of your previous post (How can I search the mailing list?).
Yeah, good question .... that's kind of problematic at the moment. Eventually we'll switch to Mailman 3 which has a sane search interface. Sorry, we don't have a better answer than using Google et al. right now.
I think this is a basic task for many geoscientists.
I did not make the effort to understand Simon Clift smarter approach. But I already implemented a quick solution using a switch for the two use cases. May I submit a pull request?
It's now come up a couple of times, so clearly there's interest. Yes, opening a PR would be very helpful.
Cheers, Ralf
On Wed, Sep 18, 2019 at 12:47 PM Dieter Werthmüller < Dieter@werthmuller.org> wrote:
Edouard,
I think this thread touched on that a bit: https://mail.python.org/pipermail/scipy-dev/2019-May/023540.html
So Simon Clift might have already some functional code for that which could serve as a starting point.
I personally would be certainly interested in the functionality.
Dieter
On 18/09/2019 11:50, Edouard Goudenhoofdt wrote:
Dear scipy developers,
One could use /scipy.interpolate.RegularGridInterpolator/ at different times with different values but the same source grid and target points. Currently this would recompute the indices each time the function is called. I would like to add the possibility to provide the target points at initialisation and the values when calling the object.
Best regards,
Edouard Goudenhoofdt //
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
-- 1129 Ibbetson Lane Mississauga, Ontario L5C 1K9 Canada
Dear Simon, Thank you for your comment. For the sake of time I just kept the current implementation while saving weights. I guess it could later be used as a reference to benchmark your approach. Best regards, Edouard On Wed, Sep 18, 2019 at 3:45 PM Simon S. Clift <ssclift@gmail.com> wrote:
Hi Edouard,
I'll keep an eye out for your pull request, or perhaps you could send me a note on your use case to my gmail address. My approach just consists of saving the interpolation weights in a sparse matrix then doing a sparse multiply for subsequent interpolations. Weights can be either multi-linear or using barycentric co-ordinates on regular grids; it's an arbitrary choice unless you have a simplicial grid layout you prefer.
Cheers -- Simon
On Wed, Sep 18, 2019 at 8:50 AM Ralf Gommers <ralf.gommers@gmail.com> wrote:
On Wed, Sep 18, 2019 at 2:33 PM Edouard Goudenhoofdt <egouden@gmail.com> wrote:
This is clearly a duplicate of your previous post (How can I search the mailing list?).
Yeah, good question .... that's kind of problematic at the moment. Eventually we'll switch to Mailman 3 which has a sane search interface. Sorry, we don't have a better answer than using Google et al. right now.
I think this is a basic task for many geoscientists.
I did not make the effort to understand Simon Clift smarter approach. But I already implemented a quick solution using a switch for the two use cases. May I submit a pull request?
It's now come up a couple of times, so clearly there's interest. Yes, opening a PR would be very helpful.
Cheers, Ralf
On Wed, Sep 18, 2019 at 12:47 PM Dieter Werthmüller < Dieter@werthmuller.org> wrote:
Edouard,
I think this thread touched on that a bit: https://mail.python.org/pipermail/scipy-dev/2019-May/023540.html
So Simon Clift might have already some functional code for that which could serve as a starting point.
I personally would be certainly interested in the functionality.
Dieter
Dear scipy developers,
One could use /scipy.interpolate.RegularGridInterpolator/ at different times with different values but the same source grid and target
On 18/09/2019 11:50, Edouard Goudenhoofdt wrote: points.
Currently this would recompute the indices each time the function is called. I would like to add the possibility to provide the target points at initialisation and the values when calling the object.
Best regards,
Edouard Goudenhoofdt //
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
-- 1129 Ibbetson Lane Mississauga, Ontario L5C 1K9 Canada _______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
The pull request has been made. On Mon, Sep 23, 2019 at 11:29 AM Edouard Goudenhoofdt <egouden@gmail.com> wrote:
Dear Simon,
Thank you for your comment. For the sake of time I just kept the current implementation while saving weights. I guess it could later be used as a reference to benchmark your approach.
Best regards,
Edouard
On Wed, Sep 18, 2019 at 3:45 PM Simon S. Clift <ssclift@gmail.com> wrote:
Hi Edouard,
I'll keep an eye out for your pull request, or perhaps you could send me a note on your use case to my gmail address. My approach just consists of saving the interpolation weights in a sparse matrix then doing a sparse multiply for subsequent interpolations. Weights can be either multi-linear or using barycentric co-ordinates on regular grids; it's an arbitrary choice unless you have a simplicial grid layout you prefer.
Cheers -- Simon
On Wed, Sep 18, 2019 at 8:50 AM Ralf Gommers <ralf.gommers@gmail.com> wrote:
On Wed, Sep 18, 2019 at 2:33 PM Edouard Goudenhoofdt <egouden@gmail.com> wrote:
This is clearly a duplicate of your previous post (How can I search the mailing list?).
Yeah, good question .... that's kind of problematic at the moment. Eventually we'll switch to Mailman 3 which has a sane search interface. Sorry, we don't have a better answer than using Google et al. right now.
I think this is a basic task for many geoscientists.
I did not make the effort to understand Simon Clift smarter approach. But I already implemented a quick solution using a switch for the two use cases. May I submit a pull request?
It's now come up a couple of times, so clearly there's interest. Yes, opening a PR would be very helpful.
Cheers, Ralf
On Wed, Sep 18, 2019 at 12:47 PM Dieter Werthmüller < Dieter@werthmuller.org> wrote:
Edouard,
I think this thread touched on that a bit: https://mail.python.org/pipermail/scipy-dev/2019-May/023540.html
So Simon Clift might have already some functional code for that which could serve as a starting point.
I personally would be certainly interested in the functionality.
Dieter
Dear scipy developers,
One could use /scipy.interpolate.RegularGridInterpolator/ at different times with different values but the same source grid and target
On 18/09/2019 11:50, Edouard Goudenhoofdt wrote: points.
Currently this would recompute the indices each time the function is called. I would like to add the possibility to provide the target points at initialisation and the values when calling the object.
Best regards,
Edouard Goudenhoofdt //
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
_______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
-- 1129 Ibbetson Lane Mississauga, Ontario L5C 1K9 Canada _______________________________________________ SciPy-Dev mailing list SciPy-Dev@python.org https://mail.python.org/mailman/listinfo/scipy-dev
participants (4)
-
Dieter Werthmüller -
Edouard Goudenhoofdt -
Ralf Gommers -
Simon S. Clift