An accurate boundary value problem solver applied to scattering from cylinders with corners

In this paper we consider the classic problems of scattering of waves from perfectly conducting cylinders with piecewise smooth boundaries. The scattering problems are formulated as integral equations and solved using a Nystr\"om scheme where the corners of the cylinders are efficiently handled by a method referred to as Recursively Compressed Inverse Preconditioning (RCIP). This method has been very successful in treating static problems in non-smooth domains and the present paper shows that it works equally well for the Helmholtz equation. In the numerical examples we specialize to scattering of E- and H-waves from a cylinder with one corner. Even at a size kd=1000, where k is the wavenumber and d the diameter, the scheme produces at least 13 digits of accuracy in the electric and magnetic fields everywhere outside the cylinder.Comment: 19 pages, 3 figure

Accuracy and Stability of Computing High-Order Derivatives of Analytic Functions by Cauchy Integrals

High-order derivatives of analytic functions are expressible as Cauchy integrals over circular contours, which can very effectively be approximated, e.g., by trapezoidal sums. Whereas analytically each radius r up to the radius of convergence is equal, numerical stability strongly depends on r. We give a comprehensive study of this effect; in particular we show that there is a unique radius that minimizes the loss of accuracy caused by round-off errors. For large classes of functions, though not for all, this radius actually gives about full accuracy; a remarkable fact that we explain by the theory of Hardy spaces, by the Wiman-Valiron and Levin-Pfluger theory of entire functions, and by the saddle-point method of asymptotic analysis. Many examples and non-trivial applications are discussed in detail.Comment: Version 4 has some references and a discussion of other quadrature rules added; 57 pages, 7 figures, 6 tables; to appear in Found. Comput. Mat

On the spectral distribution of kernel matrices related to\ud radial basis functions

This paper focuses on the spectral distribution of kernel matrices related to radial basis functions. The asymptotic behaviour of eigenvalues of kernel matrices related to radial basis functions with different smoothness are studied. These results are obtained by estimated the coefficients of an orthogonal expansion of the underlying kernel function. Beside many other results, we prove that there are exactly (k+d−1/d-1) eigenvalues in the same order for analytic separable kernel functions like the Gaussian in Rd. This gives theoretical support for how to choose the diagonal scaling matrix in the RBF-QR method (Fornberg et al, SIAM J. Sci. Comput. (33), 2011) which can stably compute Gaussian radial basis function interpolants