Multivariate polyharmonic spline interpolation