[SciPy-user] Polynomial interpolation

Anne Archibald peridot.faceted@gmail....
Mon Apr 28 15:11:23 CDT 2008


2008/4/28 Robert Kern <robert.kern@gmail.com>:
> On Mon, Apr 28, 2008 at 12:17 PM, Gael Varoquaux
>
> <gael.varoquaux@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