Related papers: Hierarchical Interpolative Factorization for Self …
The Fast Multipole Method (FMM) obeys periodic boundary conditions "natively" if it uses a periodic Green function for computing the multipole expansion in the interaction zone of each FMM oct-tree node. One can define the "optimal" Green…
Neural networks (NNs) have been widely used to solve partial differential equations (PDEs) in the applications of physics, biology, and engineering. One effective approach for solving PDEs with a fixed differential operator is learning…
In this paper, new representations of the Green's function for an acoustic d-dimensional half-space problem with impedance boundary conditions are presented. The main features of the new representation are: a) in addition to additive terms…
In this work we introduce the Dual Boson Diagrammatic Monte Carlo technique for strongly interacting electronic systems. This method combines the strength of dynamical mean-filed theory for non-perturbative description of local correlations…
We introduce Neural Green's Function, a neural solution operator for linear partial differential equations (PDEs) whose differential operators admit eigendecompositions. Inspired by Green's functions, the solution operators of linear PDEs…
We present a fast direct solver for boundary integral equations on complex surfaces in three dimensions using an extension of the recently introduced recursive strong skeletonization scheme. For problems that are not highly oscillatory, our…
The disorder averaged single-particle Green's function of electrons subject to a time-dependent random potential with long-range spatial correlations is calculated by means of bosonization in arbitrary dimensions. For static disorder our…
In this paper, we propose a fast multipole method (FMM) for 3-D linearized Poisson-Boltzmann (PB) equation in layered media. The main framework of the algorithm is analogous to the FMM for Helmholtz and Laplace equation in layered media…
We present a fast direct solver for structured linear systems based on multilevel matrix compression. Using the recently developed interpolative decomposition of a low-rank matrix in a recursive manner, we embed an approximation of the…
A novel orthogonalization-free method together with two specific algorithms are proposed to solve extreme eigenvalue problems. On top of gradient-based algorithms, the proposed algorithms modify the multi-column gradient such that earlier…
High-order derivatives of Green's functions are a key ingredient in Taylor-based fast multipole methods, Barnes-Hut $n$-body algorithms, and quadrature by expansion (QBX). In these settings, derivatives underpin either the formation,…
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…
This note proposes an efficient preconditioner for solving linear and semi-linear parabolic equations. With the Crank-Nicholson time stepping method, the algebraic system of equations at each time step is solved with the conjugate gradient…
We present an accelerated, or 'look-ahead' version of the Newton-Dinkelbach method, a well-known technique for solving fractional and parametric optimization problems. This acceleration halves the Bregman divergence between the current…
This article presents a new high-order accurate algorithm for finding a particular solution to a linear, constant-coefficient partial differential equation (PDE) by means of a convolution of the volumetric source function with the Green's…
We study the Green function of the Poisson equation in two, three and four dimensions. The solution g of the equation nabla^2 g(x - x') = delta^(D)(x - x'), where x and x are D-dimensional position vectors, is customarily expanded into…
The linear Boltzmann equation with constant coefficients in the three-dimensional infinite space is revisited. It is known that the Green's function can be calculated via the Fourier transform in the case of isotropic scattering. In this…
An efficient algorithm for solving Poisson's equation in two and three spatial dimensions is discussed. The algorithm, which is described in detail, is based on the integral form of Poisson's equation and utilizes spherical coordinates and…
Wereportonanewmultiscalemethodapproachforthestudyofsystemswith wide separation of short-range forces acting on short time scales and long-range forces acting on much slower scales. We consider the case of the Poisson-Boltzmann equation that…
The stability of a differentially rotating fluid subject to its own gravity is a problem with applications across wide areas of astrophysics--from protoplanetary discs (PPDs) to entire galaxies. The shearing box formalism offers a…