[SciPy-user] Polynomial interpolation

Eike Welk eike.welk@gmx....
Mon Apr 28 17:04:16 CDT 2008


On Monday 28 April 2008 23:24, Robert Kern wrote:
> On Mon, Apr 28, 2008 at 4:21 PM, Eike Welk <eike.welk@gmx.net> 
wrote:
> > On Monday 28 April 2008 22:36, Anne Archibald wrote:
> >  > 2008/4/28 Robert Kern <robert.kern@gmail.com>:
> >  > >  Well those certainly aren't useful. The only functions I
> >  > > would consider adding would be "one-shot" functions, e.g.:
> >  > >
> >  > >   def krogh(xi, yi, x):
> >  > >       return KroghInterpolator(xi,yi)(x)
> >  >
> >  > The problem here is that construction of the splines is an
> >  > order degree**2 process, so I want an interface that
> >  > encourages users to construct them once and for all. I think
> >  > such an approach also discourages people from just
> >  >
> >  > y_interp = krogh(all_my_data_x, all_my_data_y, x_interp)
> >  >
> >  > with hundreds of points, the results of which will be
> >  > meaningless and horrible.
> >
> >  You could store already constructed interpolation objects in a
> >  dictionary. (I didn't test it.):
>
> Arrays are unhashable and cannot be used as dictionary keys.
Ah, yes! 

One could however store a copy of the last x-data and y-data together 
with the interpolator. When x-data and y-data are the same at the 
next call, the interpolator could be reused.  

Kind regards,
Eike.


More information about the SciPy-user mailing list