interpolation in N-D box grids, ~ 1M/sec in 4d, scipy.ndimage.map_coordinates + np.interp