Fast convolution with free-space Green's functions

Abstract

We introduce a fast algorithm for computing volume potentials - that is, the convolution of a translation invariant, free-space Green's function with a compactly supported source distribution defined on a uniform grid. The algorithm relies on regularizing the Fourier transform of the Green's function by cutting off the interaction in physical space beyond the domain of interest. This permits the straightforward application of trapezoidal quadrature and the standard FFT, with superalgebraic convergence for smooth data. Moreover, the method can be interpreted as employing a Nystrom discretization of the corresponding integral operator, with matrix entries which can be obtained explicitly and rapidly. This is of use in the design of preconditioners or fast direct solvers for a variety of volume integral equations. The method proposed permits the computation of any derivative of the potential, at the cost of an additional FFT.

21 pages, 10 Figures, 7 Tables

Cited by in corpus (29)

• An adaptive fast multipole accelerated Poisson solver for complex geometries
• Fast Ewald summation for free-space Stokes potentials
• Partition of Unity Extension of Functions on Complex Domains
• Predicting nonlinear optical scattering with physics-driven neural networks
• FLUPS: a Fourier-based Library of Unbounded Poisson Solvers
• High-precision real-space simulation of electrostatically-confined few-electron states
• The Spectral Ewald method for singly periodic domains
• Fast Algorithm for Quasi-2D Coulomb Systems
• Fast Ewald summation for Stokes flow with arbitrary periodicity
• Non-singular Green's functions for the unbounded Poisson equation in one, two and three dimensions
• Fast Ewald summation for Green's functions of Stokes flow in a half-space
• Improved convergence of fast integral equation solvers for acoustic scattering by inhomogeneous penetrable media with discontinuous material interface
• Simulations of Future Particle Accelerators: Issues and Mitigations
• An iterative domain decomposition, spectral finite element method on non-conforming meshes suitable for high frequency Helmholtz problems
• Unified treatment of exact and approximate scalar electromagnetic wave scattering
• Fast Ewald summation for electrostatic potentials with arbitrary periodicity
• Computation of Volume Potentials on Structured Grids Using the Method of Local Corrections
• Optimal control of the self-bound dipolar droplet formation process
• Fast accurate approximation of convolutions with weakly singular kernel and its applications
• Sparsify and sweep: an efficient preconditioner for the Lippmann-Schwinger equation
• An integral equation method for the advection-diffusion equation on time-dependent domains in the plane
• Reduction of multivariate mixtures and its applications
• A fast adaptive algorithm for scattering from a two dimensional radially-symmetric potential
• An efficient spectral Poisson solver for the nirvana-III code: the shearing-box case with vertical vacuum boundary conditions
• Diffraction Tomography with Helmholtz Equation: Efficient and Robust Multigrid-Based Solver
• Spectrally accurate Ewald summation for the Yukawa potential in two dimensions
• Domain decomposition of the modified Born series approach for large-scale wave propagation simulations
• Vico-Greengard-Ferrando quadratures in the tensor solver for integral equations
• A Massively Parallel Performance Portable Free-space Spectral Poisson Solver