Charles R Harris
Thu Feb 21 07:50:25 CST 2013
On Thu, Feb 21, 2013 at 3:33 AM, Pablo Winant <email@example.com>wrote:
> 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
> 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?
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-Dev