Fast algorithms for Quadrature by Expansion I: Globally valid expansions
arXiv:1602.05301 · doi:10.1016/j.jcp.2017.04.062
Abstract
The use of integral equation methods for the efficient numerical solution of PDE boundary value problems requires two main tools: quadrature rules for the evaluation of layer potential integral operators with singular kernels, and fast algorithms for solving the resulting dense linear systems. Classically, these tools were developed separately. In this work, we present a unified numerical scheme based on coupling Quadrature by Expansion, a recent quadrature method, to a customized Fast Multipole Method (FMM) for the Helmholtz equation in two dimensions. The method allows the evaluation of layer potentials in linear-time complexity, anywhere in space, with a uniform, user-chosen level of accuracy as a black-box computational method. Providing this capability requires geometric and algorithmic considerations beyond the needs of standard FMMs as well as careful consideration of the accuracy of multipole translations. We illustrate the speed and accuracy of our method with various numerical examples. Keywords: Layer Potentials; Singular Integrals; Quadrature; High-order accuracy; Integral equations; Helmholtz equation; Fast multipole method.
References in corpus (2)
Cited by in corpus (20)
- A fast integral equation method for solid particles in viscous flow using quadrature by expansion
- An adaptive fast multipole accelerated Poisson solver for complex geometries
- A Fast Algorithm for Quadrature by Expansion in Three Dimensions
- A local target specific quadrature by expansion method for evaluation of layer potentials in 3D
- Adaptive quadrature by expansion for layer potential evaluation in two dimensions
- Estimation of quadrature errors in layer potential evaluation using quadrature by expansion
- Taylor States in Stellarators: A Fast High-order Boundary Integral Solver
- A Fast Algorithm with Error Bounds for Quadrature by Expansion
- Optimization of Fast Algorithms for Global Quadrature by Expansion Using Target-Specific Expansions
- Efficient high-order singular quadrature schemes in magnetic fusion
- A robust solver for elliptic PDEs in 3D complex geometries
- Highly accurate special quadrature methods for Stokesian particle suspensions in confined geometries
- An Integral Equation Method for the Cahn-Hilliard Equation in the Wetting Problem
- Integral equation formulation of the biharmonic Dirichlet problem
- Conformal Mapping via a Density Correspondence for the Double-Layer Potential
- Quadrature by Two Expansions: Evaluating Laplace Layer Potentials using Complex Polynomial and Plane Wave Expansions
- On the Approximation of Local Expansions of Laplace Potentials by the Fast Multipole Method
- Multiscale Hydrophobic Lipid Dynamics Simulated by Efficient Integral Equation Methods
- Finite elements for Helmholtz equations with a nonlocal boundary condition
- Accurate Derivative Evaluation for any Grad-Shafranov Solver