A fast and well-conditioned spectral method for singular integral equations

详细信息    查看全文

• 作者：Richard Mikael Slevinsky^a ; ^{Richard.Slevinsky@umanitoba.ca} ; Sheehan Olver^b ; ^{Sheehan.Olver@sydney.edu.au}
• 关键词：Spectral method ; Ultraspherical polynomials ; Singular integral equations
• 刊名：Journal of Computational Physics
• 出版年：2017
• 出版时间：1 March 2017
• 年：2017
• 卷：332
• 期：Complete
• 页码：290-315
• 全文大小：2331 K
• 卷排序：332

文摘

We develop a spectral method for solving univariate singular integral equations over unions of intervals by utilizing Chebyshev and ultraspherical polynomials to reformulate the equations as almost-banded infinite-dimensional systems. This is accomplished by utilizing low rank approximations for sparse representations of the bivariate kernels. The resulting system can be solved in O(m2n)O(m2n) operations using an adaptive QR factorization, where m is the bandwidth and n   is the optimal number of unknowns needed to resolve the true solution. The complexity is reduced to O(mn)O(mn) operations by pre-caching the QR factorization when the same operator is used for multiple right-hand sides. Stability is proved by showing that the resulting linear operator can be diagonally preconditioned to be a compact perturbation of the identity. Applications considered include the Faraday cage, and acoustic scattering for the Helmholtz and gravity Helmholtz equations, including spectrally accurate numerical evaluation of the far- and near-field solution. The Julia software package SingularIntegralEquations.jl implements our method with a convenient, user-friendly interface.