Related papers: Fast summation of Stokes potentials using a new ke…
Kernel-based methods are heavily used in machine learning. However, they suffer from $O(N^2)$ complexity in the number $N$ of considered data points. In this paper, we propose an approximation procedure, which reduces this complexity to…
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…
We focus on a recently developed generalized pseudospectral method for accurate, efficient treatment of certain central potentials of interest in various branches in quantum mechanics, usually having singularity. Essentially this allows…
A highly accurate method for simulating surfactant-covered droplets in two-dimensional Stokes flow with solid boundaries is presented. The method handles both periodic channel flows of arbitrary shape and stationary solid constrictions. A…
This work describes a new version of the Fast Multipole Method for summing pairwise particle interactions that arise from discretizing integral transforms and convolutions on the sphere. The kernel approximations use barycentric Lagrange…
Quantile computation is a core primitive in large-scale data analytics. In Spark, practitioners typically rely on the Greenwald-Khanna (GK) Sketch, an approximate method. When exact quantiles are required, the default option is an expensive…
With the aim of describing bound and continuum states for diatomic molecules, we develop and implement a spectral method that makes use of Generalized Sturmian Functions (GSF) in prolate spheroidal coordinates. In order to master all…
Kohn-Sham (KS) density functional theory (DFT) is a very efficient method for calculating various properties of solids as, for instance, the total energy, the electron density, or the electronic band structure. The KS-DFT method leads to…
In this paper, a fast multipole method (FMM) is proposed to compute long-range interactions of wave sources embedded in 3-D layered media. The layered media Green's function for the Helmholtz equation, which satisfies the transmission…
The continuum random-phase approximation is extended to the one applicable to deformed nuclei. We propose two different approaches. One is based on the use of the three dimensional (3D) Green's function and the other is the small-amplitude…
We present a fast multipole method (FMM) for solving Maxwell's equations in three-dimensional (3-D) layered media, based on the magnetic vector potential $\boldsymbol A$ under the Lorenz gauge, to derive the layered dyadic Green's function.…
A unified treatment for fast and spectrally accurate evaluation of electrostatic potentials subject to periodic boundary conditions in any or none of the three spatial dimensions is presented. Ewald decomposition is used to split the…
We present a unified treatment of the Fourier spectra of spherically symmetric nonlocal diffusion operators. We develop numerical and analytical results for the class of kernels with weak algebraic singularity as the distance between source…
We utilize the domain integral equation formulation to simulate two-dimensional transverse electric scattering in a homogeneous medium and a summation of modulated Gaussian functions to approximate the dual Gabor window. Then we apply Ewald…
Coulomb interaction, following an inverse-square force-law, quantifies the amount of force between two stationary and electrically charged particles. The long-range nature of Coulomb interactions poses a major challenge to molecular…
Kernel methods are a highly effective and widely used collection of modern machine learning algorithms. A fundamental limitation of virtually all such methods are computations involving the kernel matrix that naively scale quadratically…
The fast multipole method (FMM) performs fast approximate kernel summation to a specified tolerance $\epsilon$ by using a hierarchical division of the domain, which groups source and receiver points into regions that satisfy local…
The integral equation method is widely used in numerical simulations of 2D/3D acoustic and electromagnetic scattering problems, which needs a large number of values of the Green's functions. A significant topic is the scattering problems in…
The Fast Multipole Method (FMM) computes pairwise interactions between particles with an efficiency that scales linearly with the number of particles. The method works by grouping particles based on their spatial distribution and…
This paper introduces a fast algorithm, applicable throughout the electromagnetic spectrum, for the numerical solution of problems of scattering by periodic surfaces in two-dimensional space. The proposed algorithm remains highly accurate…