[SciPy-user] Polynomial interpolation
Anne Archibald
peridot.faceted at gmail.com
Mon Apr 28 16:11:23 EDT 2008
2008/4/28 Robert Kern <robert.kern at gmail.com>:
> On Mon, Apr 28, 2008 at 12:17 PM, Gael Varoquaux
>
> <gael.varoquaux at normalesup.org> wrote:
>
> > I am not suggesting to make a bad design choice, or to go
> > out of our way, I just want a two liner helper function.
>
> I'm afraid that *is* widely recognized as a bad design choice.
Well, I think the way to go is to mimic splrep/splev:
def krogh_rep(xi, yi):
return KroghInterpolator(xi, yi)
def krogh_ev(kr, x, der=0):
if der==0:
return kr(x)
else:
return kr.derivative(x, der=der)
The only danger is that if it ever raises an exception and the users
see how simple it is they will feel like idiots.
Anne
More information about the SciPy-User
mailing list