fix a bug in scipy where interpolation with map_coordinates and order=1
yields a segfault if the array is float32 instead of float64
also, try to use a subset of the array for coordinate mapping if
possible