Singular quadrature rules and fast convolutions for Fourier spectral methods

arXiv1211.6074MaRDI QIDQ6237489FDOQ6237489

Authors: Jae-Seok Huh, George Fann 

Publication date: 26 November 2012

Abstract: We present a generic scheme to construct corrected trapezoidal rules with spectral accuracy for integral operators with weakly singular kernels in arbitrary dimensions. We assume that the kernel factorization of the form, $K=alphaphi+widetildeK$ with smooth $alpha$ and $widetildeK$, is available so that the operations on the smooth factors can be performed accurately on the basis of standard Fourier spectral methods. To achieve high precision results, our approach utilizes the exact evaluation of the Fourier coefficients of the radial singularity $phi$, which can be obtained in arbitrary dimensions by the singularity isolation/truncation described in this article. We provide a complete set of formulas for singularities of the type: $log(r)$ and $r^{-u}$. Convergence analysis shows that the constructed quadrature rules exhibit almost identical rate of convergence to the trapezoidal rule applied for non-singular integrands. Especially, for smooth data, the corrected trapezoidal rules converge super-algebraically.

This page was built for publication: Singular quadrature rules and fast convolutions for Fourier spectral methods

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6237489)