Interpolating multidimensional normal spline