Optimization of Fast Algorithms for Global Quadrature by Expansion Using Target-Specific Expansions

Abstract

We develop an algorithm for the asymptotically fast evaluation of layer potentials close to and on the source geometry, combining Geometric Global Accelerated QBX (`GIGAQBX') and target-specific expansions. GIGAQBX is a fast high-order scheme for evaluation of layer potentials based on Quadrature by Expansion (`QBX') using local expansions formed via the Fast Multipole Method (FMM). Target-specific expansions serve to lower the cost of the formation and evaluation of QBX local expansions, reducing the associated computational effort from $O((p+1{)}^2)$ to $O(p+1)$ in three dimensions, without any accuracy loss compared with conventional expansions, but with the loss of source/target separation in the expansion coefficients. GIGAQBX is a `global' QBX scheme, meaning that the potential is mediated entirely through expansions for points close to or on the boundary. In our scheme, this single global expansion is decomposed into two parts that are evaluated separately: one part incorporating near-field contributions using target-specific expansions, and one part using conventional spherical harmonic expansions of far-field contributions, noting that convergence guarantees only exist for the sum of the two sub-expansions. By contrast, target-specific expansions were originally introduced as an acceleration mechanism for `local' QBX schemes, in which the far-field does not contribute to the QBX expansion. Compared with the unmodified GIGAQBX algorithm, we show through a reproducible, time-calibrated cost model that the combined scheme yields a considerable cost reduction for the near-field evaluation part of the computation. We support the effectiveness of our scheme through numerical results demonstrating performance improvements for Laplace and Helmholtz kernels.

References in corpus (7)

• Quadrature by Expansion: A New Method for the Evaluation of Layer Potentials
• Fast algorithms for Quadrature by Expansion I: Globally valid expansions
• 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
• Estimation of quadrature errors in layer potential evaluation using quadrature by expansion
• A Fast Algorithm with Error Bounds for Quadrature by Expansion
• Optimizing the adaptive fast multipole method for fractal sets

Cited by in corpus (8)

• A Fast Integral Equation Method for the Two-Dimensional Navier-Stokes Equations
• Scalable Simulation of Realistic Volume Fraction Red Blood Cell Flows through Vascular Networks
• A robust solver for elliptic PDEs in 3D complex geometries
• Fast multipole methods for evaluation of layer potentials with locally-corrected quadratures
• Highly accurate special quadrature methods for Stokesian particle suspensions in confined geometries
• A fast, high-order scheme for evaluating volume potentials on complex 2D geometries via area-to-line integral conversion and domain mappings
• 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