Abstract:The quasi-2D electrostatic systems, characterized by periodicity in two dimensions with a free third dimension, have garnered significant interest in many fields. We apply the sum-of-Gaussians (SOG) approximation to the Laplace kernel, dividing the interactions into near-field, mid-range, and long-range components. The near-field component, singular but compactly supported in a local domain, is directly calculated. The mid-range component is managed using a procedure similar to nonuniform fast Fourier transforms in three dimensions. The long-range component, which includes Gaussians of large variance, is treated with polynomial interpolation/anterpolation in the free dimension and Fourier spectral solver in the other two dimensions on proxy points. Unlike the fast Ewald summation, which requires extensive zero padding in the case of high aspect ratios, the separability of Gaussians allows us to handle such case without any zero padding in the free direction. Furthermore, while NUFFTs typically rely on certain upsampling in each dimension, and the truncated kernel method introduces an additional factor of upsampling due to kernel oscillation, our scheme eliminates the need for upsampling in any direction due to the smoothness of Gaussians, significantly reducing computational cost for large-scale problems. Finally, whereas all periodic fast multipole methods require dividing the periodic tiling into a smooth far part and a near part containing its nearest neighboring cells, our scheme operates directly on the fundamental cell, resulting in better performance with simpler implementation. We provide a rigorous error analysis showing that upsampling is not required in NUFFT-like steps, achieving $O(N\log N)$ complexity with a small prefactor. The performance of the scheme is demonstrated via extensive numerical experiments.

An adaptive fast Gauss transform in two dimensions

A new version of the adaptive fast Gauss transform for discrete and continuous sources

A Fast Adaptive Method for the Heat Equation with Moving or Free Boundaries in One Dimension

Parallel Conjugate Gradient Method Based on Spline Difference Method for the One-Dimensional Heat Equation

Fast integral equation methods for linear and semilinear heat equations in moving domains

Fast Adaptive Fourier Integration for Spectral Densities of Gaussian Processes

A fast spectral sum-of-Gaussians method for electrostatic summation in quasi-2D systems

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

A Dual-space Multilevel Kernel-splitting Framework for Discrete and Continuous Convolution

Fast Computation of the Multidimensional Discrete Fourier Transform and Discrete Backward Fourier Transform on Sparse Grids

Efficient Gaussian Process Calculations Using Chebyshev Nodes and Fast Fourier Transform

A Dynamic Low-Rank Fast Gaussian Transform

A Fast Fourier-Galerkin Method Solving Boundary Integral Equations for the Helmholtz Equation with Exponential Convergence

Combining triangle Gaussian integration and modified NUFFT for evaluating two-dimensional Fourier transform integrals

Fast multiscale Gaussian beam methods for wave equations in bounded convex domains

Fast Fourier transforms for the evaluation of convolution products: CPU versus GPU implementation

Highly efficient Gauss's law-preserving spectral algorithms for Maxwell's double-curl source and eigenvalue problems based on eigen-decomposition

Fast computation of general Fourier Transforms on GPUS

Uniform approximation of common Gaussian process kernels using equispaced Fourier grids

Fast Fourier Transforms for Spherical Gauss-Laguerre Basis Functions