Related papers: Computing the gravitational potential on nested me…
We present a numerical method for solving the Poisson equation on a nested grid. The nested grid consists of uniform grids having different grid spacing and is designed to cover the space closer to the center with a finer grid. Thus our…
We present a simple and effective multigrid-based Poisson solver of second-order accuracy in both gravitational potential and forces in terms of the one, two and infinity norms. The method is especially suitable for numerical simulations…
We introduce the concept of compressed convolution, a technique to convolve a given data set with a large number of non-orthogonal kernels. In typical applications our technique drastically reduces the effective number of computations. The…
We present a generalized algorithm based on a spherical harmonics expansion method for efficient computation of the three-dimensional gravitational potential on a multi-patch grid in spherical geometry. Instead of solving for the…
The determination of the gravitational potential by the polyhedral method is revisited in the case where the surface of a body is composed of triangular facets. Based upon six test-shapes of astrophysical interest (sphere, spheroid,…
The Poisson equation is commonly encountered in engineering, for instance in computational fluid dynamics (CFD) where it is needed to compute corrections to the pressure field to ensure the incompressibility of the velocity field. In the…
We develop a method for reconstructing the two-dimensional, projected gravitational potential of galaxy clusters from observed line-of-sight velocity dispersions of cluster galaxies. It is the third of an intended series of papers aiming at…
We propose an efficient algorithm for the evaluation of the potential and its gradient of gravitational/electrostatic $N$-body systems, which we call particle mesh multipole method (PMMM or PM$^3$). PMMM can be understood both as an…
We present a coarse-grid projection (CGP) method for accelerating incompressible flow computations, which is applicable to methods involving Poisson equations as incompressibility constraints. The CGP methodology is a modular approach that…
Coarse grid projection (CGP) is a multiresolution technique for accelerating numerical calculations associated with a set of nonlinear evolutionary equations along with the stiff Poisson equations. In this article we use CGP for the first…
I describe a nested-grid particle-mesh (NGPM) code designed to study gravitational instability in three-dimensions. The code is based upon a standard PM code. Within the parent grid I am able to define smaller sub-grids allowing us to…
The gravitational potential is a key function involved in many astrophysical problems. Its evaluation inside continuous media from Newton's law is known to be challenging because of the diverging kernel 1/|r-r'|. This difficulty is…
We extend the work of Yen et al. (2012) and develop 2nd order formulae to accommodate a nested grid discretization for the direct self-gravitational force calculation for infinitesimally thin gaseous disks. This approach uses a…
We calculate ab initio the gravitational potential energy per unit area for a gravitationally coupled multi-component galactic disk of stars and gas, which is given as the integration over vertical density distribution, vertical…
We present a versatile formulation of the convolution operation that we term a "mapped convolution." The standard convolution operation implicitly samples the pixel grid and computes a weighted sum. Our mapped convolution decouples these…
The paper introduces the weighted convolution, a novel approach to the convolution for signals defined on regular grids (e.g., 2D images) through the application of an optimal density function to scale the contribution of neighbouring…
Coarse Grid Projection (CGP) methodology is used to accelerate the computations of sets of decoupled nonlinear evolutionary and linear static equations. In CGP, the linear equations are solved on a coarsened mesh compared to the nonlinear…
We present a practical way of introducing convolutional structure into Gaussian processes, making them more suited to high-dimensional inputs like images. The main contribution of our work is the construction of an inter-domain inducing…
We present a polynomial multigrid method for the nodal interior penalty formulation of the Poisson equation on three-dimensional Cartesian grids. Its key ingredient is a weighted overlapping Schwarz smoother operating on element-centered…
Accurate parameter estimation of gravitational waves from coalescing compact binary sources is a key requirement for gravitational-wave astronomy. Evaluating the posterior probability density function of the binary's parameters (component…