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

From MaRDI portal
Publication:680105

DOI10.1016/J.JCP.2016.12.009zbMATH Open1380.65446arXiv1507.00596OpenAlexW2963567562MaRDI QIDQ680105FDOQ680105


Authors: Richard M. Slevinsky, Sheehan Olver Edit this on Wikidata


Publication date: 22 January 2018

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: 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 calO(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 calO(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.


Full work available at URL: https://arxiv.org/abs/1507.00596




Recommendations




Cites Work


Cited In (23)

Uses Software





This page was built for publication: A fast and well-conditioned spectral method for singular integral equations

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