interpol/intergrid: fast N-d map_coordinates + non-uniform grid