On the Approximation of Local Expansions of Laplace Potentials by the Fast Multipole Method

by   Matt Wala, et al.

In this paper, we present a generalization of the classical error bounds of Greengard-Rokhlin for the Fast Multipole Method (FMM) for Laplace potentials in three dimensions, extended to the case of local expansion (instead of point) targets. We also present a complementary, less sharp error bound proven via approximation theory whose applicability is not restricted to Laplace potentials. Our study is motivated by the GIGAQBX FMM, an algorithm for the fast, high-order accurate evaluation of layer potentials near and on the source layer. GIGAQBX is based on the FMM, but unlike a conventional FMM, which is designed to evaluate potentials at point-shaped targets, GIGAQBX evaluates local expansions of potentials at ball-shaped targets. Although the accuracy (or the acceleration error, i.e., error due to the approximation of the potential by the fast algorithm) of the conventional FMM is well understood, the acceleration error of FMM-based algorithms applied to the evaluation of local expansions has not been as well studied. The main contribution of this paper is a proof of a set of hypotheses first demonstrated numerically in the paper "A Fast Algorithm for Quadrature by Expansion in Three Dimensions," which pertain to the accuracy of FMM approximation of local expansions of Laplace potentials in three dimensions. These hypotheses are also essential to the three-dimensional error bound for GIGAQBX, which was previously stated conditionally on their truth and can now be stated unconditionally.



There are no comments yet.


page 1

page 2

page 3

page 4


Accurate quadrature of nearly singular line integrals in two and three dimensions by singularity swapping

The method of Helsing and co-workers evaluates Laplace and related layer...

Quadrature by Two Expansions: Evaluating Laplace Layer Potentials using Complex Polynomial and Plane Wave Expansions

The recently developed quadrature by expansion (QBX) technique accuratel...

High Performance Evaluation of Helmholtz Potentials usingthe Multi-Level Fast Multipole Algorithm

Evaluation of pair potentials is critical in a number of areas of physic...

High Performance Evaluation of Helmholtz Potentials using the Multi-Level Fast Multipole Algorithm

Evaluation of pair potentials is critical in a number of areas of physic...

Quadrature error estimates for layer potentials evaluated near curved surfaces in three dimensions

The quadrature error associated with a regular quadrature rule for evalu...

Piece-wise quadratic approximations of arbitrary error functions for fast and robust machine learning

Most of machine learning approaches have stemmed from the application of...

Fast localization of eigenfunctions via smoothed potentials

We study the problem of predicting highly localized low-lying eigenfunct...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.