\input zb-basic
\input zb-ioport
\iteman{io-port 06071015}
\itemau{Fasshauer, Gregory E.; McCourt, Michael J.}
\itemti{Stable evaluation of Gaussian radial basis function interpolants.}
\itemso{SIAM J. Sci. Comput. 34, No. 2, A737-A762 (2012).}
\itemab
Summary: We provide a new way to compute and evaluate Gaussian radial basis function interpolants in a stable way with a special focus on small values of the shape parameter, i.e., for ``flat'' kernels. This work is motivated by the fundamental ideas proposed earlier by Bengt Fornberg and his coworkers [{\it B. Fornberg}, {\it E. Larsson} and {\it N. Flyer}, SIAM J. Sci. Comput. 33, No. 2, 869--892 (2011; Zbl 1227.65018); {\it B. Fornberg} and {\it C. Piret}, SIAM J. Sci. Comput. 30, No. 1, 60--80 (2007; Zbl 1159.65307)]. However, following Mercer's theorem, an $L_2(\mathbb{R}^d, \rho)$-orthonormal expansion of the Gaussian kernel allows us to come up with an algorithm that is simpler than the one proposed by [Fornberg, Larsson, and Flyer, loc. cit.] and that is applicable in arbitrary space dimensions d. In addition to obtaining an accurate approximation of the radial basis function interpolant (using many terms in the series expansion of the kernel), we also propose and investigate a highly accurate least-squares approximation based on early truncation of the kernel expansion.
\itemrv{~}
\itemcc{}
\itemut{radial basis functions; Gaussian kernel; stable evaluation; Mercer's theorem; eigenfunction expansion; QR decomposition}
\itemli{doi:10.1137/110824784}
\end