ENH: Efficient interpolation on regular grids in arbitrary dimensions