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


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