Local Accuracy for Radial Basis Function Interpolation on Finite Uniform Grids

We consider interpolation on a finite uniform grid by means of one of the radial basis functions (RBF) φ(r)=r^γ for γ>0, γ2 or φ(r)=r^γ ln r for γ2₊. For each positive integer N, let h=N⁻¹ and let {x_i: i =1, 2, …, (N+1)^d} be the set of vertices of the uniform grid of mesh-size h on the unit d-dimensional cube [0, 1]^d. Given f: [0, 1]^d→, let s_h be its unique RBF interpolant at the grid vertices: s_h(x_i)=f(x_i), i=1, 2, …, (N+1)^d. For h→0, we show that the uniform norm of the error f−s_h on a compact subset K of the interior of [0, 1]^d enjoys the same rate of convergence to zero as the error of RBF interpolation on the infinite uniform grid h^d, provided that f is a data function whose partial derivatives in the interior of [0, 1]^d up to a certain order can be extended to Lipschitz functions on [0, 1]^d.