An analytic approximation to the cardinal functions of Gaussian radial basis functions on an infinite lattice

Gaussian radial basis functions (RBFs) have been very useful in computer graphics and for numerical solutions of partial differential equations where these RBFs are defined, on a grid with uniform spacing h, as translates of the ''master'' function @f(x;@a,h)=exp(-[@a^2/h^2]x^2) where @a is a user-choosable constant. Unfortunately, computing the coefficients of @f(x-jh;@a,h) requires solving a linear system with a dense matrix. It would be much more efficient to rearrange the basis functions into the equivalent ''Lagrangian'' or ''cardinal'' basis because the interpolation matrix in the new basis is the identity matrix; the cardinal basis C"j(x;@a,h) is defined by the set of linear combinations of the Gaussians such that C"j(kh)=1 when k=j and C"j(kh)=0 for all integers k j. We show that the cardinal functions for the uniform grid are C"j(x;h,@a)=C(x/h-j;@a) where C(X;@a)~(@a^2/@p)sin(@pX)/sinh(@a^2X). The relative error is only about 4exp([email protected]^2/@a^2) as demonstrated by the explicit second order approximation. It has long been known that the error in a series of Gaussian RBFs does not converge to zero for fixed @a as h->0, but only to an ''error saturation'' proportional to exp([email protected]^2/@a^2). Because the error in our approximation to the master cardinal function C(X;@a) is the square of the error saturation, there is no penalty for using our new approximations to obtain matrix-free interpolating RBF approximations to an arbitrary function f(x). The master cardinal function on a uniform grid in d dimensions is just the direct product of the one-dimensional cardinal functions. Thus in two dimensions C(X,Y;@a)~(@a^4/@p^2)sin(@pX)sin(@pY)/[sinh(@a^2X)sinh(@a^2Y)]. We show that the matrix-free interpolation can be extended to non-uniform grids by a smooth change of coordinates.