[SciPy-User] fast N-d interpolators Intergrid and Barypol
denis
denis-bz-py at t-online.de
Tue Mar 12 13:14:08 EDT 2013
Folks,
two small fast interpolators:
[Intergrid](http://denis-bz.github.com/docs/intergrid.html):
interpolate in an N-d box grid, uniform or non-uniform.
This is just a wrapper for scipy.ndimage.map_coordinates and numpy.interp.
[Barypol](http://denis-bz.github.com/docs/barypol.html):
interpolate in a uniform N-d box grid, using d + 1 corners of a simplex
(triangle, tetrahedron ...) around each query point.
>From Munos and Moore, "Variable Resolution Discretization in Optimal Control",
1999, 24p; see the pictures and example on pp. 4-5.
It's implemented in a C++ header barypol.h, with a Cython wrapper
(both of which more knowledgeable people could certainly improve).
In 4d, 5d, 6d, Intergrid does around 3M, 2M, .8M interpolations / second.
Barypol is ~ 5 times faster, but not as smooth.
(Both will of course become cache-bound for large grids.)
See interpol/test/*.log
and http://github.com/denis-bz/interpol .
Comments are welcome, testcases most welcome.
cheers
-- denis
More information about the SciPy-User
mailing list