[SciPy-user] scipy.interpolate.Rbf() MemoryError for more than 4000 (x, y) pairs ?

[josef.pktd at gmail.com]

On Fri, Apr 3, 2009 at 12:39 PM, Pauli Virtanen <pav at iki.fi> wrote:
> Fri, 03 Apr 2009 12:22:26 -0400, josef.pktd kirjoitti:
> [clip]
>> Besides the memory error, I don't think rbf can do this, rbf uses all
>> points to do the interpolation, the kernel is number of observation
>> squared. In some examples, that I tried, if the number of points is to
>> large, then you don't get a smooth solution, the interpolation has too
>> much variation (like fitting a very high order polynomial).
>
> Btw, the different Rbf functions (thin-plate, linear, etc.) seem to
> behave slightly differently in how stably they can fit many points. I'd
> imagine why this happens is answered somewhere in the literature...
>

I had only tried and looked at the code for gaussian and multiquadratic.
linear and thin-plate work in my example, but according to the
comments 'thin-plate' is the first that has a memory error.

Josef
-------------- next part --------------
A non-text attachment was scrubbed...
Name: rbf100_thinplate.png
Type: image/png
Size: 43707 bytes
Desc: not available
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20090403/f29f1fb9/attachment.png>