A new class of oscillatory radial basis functions

Radial basis functions (RBFs) form a primary tool for multivariate interpolation, and they are also receiving increased attention for solving PDEs on irregular domains. Traditionally, only nonoscillatory radial functions have been considered. We find here that a certain class of oscillatory radial functions (including Gaussians as a special case) leads to nonsingular interpolants with intriguing features especially as they are scaled to become increasingly flat. This flat limit is important in that it generalizes traditional spectral methods to completely general node layouts. Interpolants based on the new radial functions appear immune to many or possibly all cases of divergence that in this limit can arise with other standard types of radial functions (such as multiquadrics and inverse multiquadratics).