[SciPy-Dev] Improvement to regular grid interpolation

[Edouard Goudenhoofdt]

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 at 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 at gmail.com>
> wrote:
>
>>
>>
>> On Wed, Sep 18, 2019 at 2:33 PM Edouard Goudenhoofdt <egouden at 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 at 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 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
>>>
>> _______________________________________________
>> SciPy-Dev mailing list
>> SciPy-Dev at python.org
>> https://mail.python.org/mailman/listinfo/scipy-dev
>>
>
>
> --
> 1129 Ibbetson Lane
> Mississauga, Ontario
> L5C 1K9       Canada
> _______________________________________________
> SciPy-Dev mailing list
> SciPy-Dev at python.org
> https://mail.python.org/mailman/listinfo/scipy-dev
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20190923/89927205/attachment-0001.html>