[SciPy-Dev] Splines

[Charles R Harris]

On Thu, Feb 21, 2013 at 3:33 AM, Pablo Winant <pablo.winant at gmail.com>wrote:

> Hi,
>
> It's good to see people working in that direction. I will be more than
> happy to participate if I can.
>
> As a user : I really wish I had a cubic spline interpolation with
> natural bounding conditions (which implies linear extrapolation). This
> is what is implemented in Matlab's griddedInterpolator and is missing
> from the existing options in scipy (right?). Being able to evaluate the
> derivatives is also a big advantage.
>
> As a developper : there is the Einspline library that has a
> straightforward implementation in C. In particular, the representation
> of spline objects (C structs) and the low-level API are very
> straightforward and may be a good of inspiration. The library is
> currently GPL, but the author told me it could be made BSD if some m4
> macros are removed.
> For what it's worth, I have made a simple Cython wrapper around some of
> its functions
> (https://github.com/albop/dolo/tree/master/dolo/numeric/interpolation)
> and was considering repackaging it. I had some other plans (like writing
> code for higher dimensions using some kind of templating, and updating
> SIMD evaluation). If there is a more elegant alternative, I would be
> happy to jump on the bandwagon.
> I agree completely about the remarks on the spline format: it should be
> left as simple as possible. Having it isomorphic to C structs would be a
> good thing too as it permits easy experiments with extensions (for
> instance PyCuda kernels)
>

C structs tend to be non-portable, although we could work around that. Do
you use PyCuda? Is it mostly float32?

<snip>

Chuck
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-dev/attachments/20130221/7041bdd2/attachment.html>