Related papers: Yet another fast multipole method without multipol…
Particle tracing through numerical integration is a well-known approach to generating pathlines for visualization. However, for particle simulations, the computation of pathlines is expensive, since the interpolation method is complicated…
An algorithm for fast calculation of the Coulombic forces and energies of point particles with free boundary conditions is proposed. Its calculation time scales as N log N for N particles. This novel method has lower crossover point with…
We present a novel approach to the parallelization of the parabolic fast multipole method for a space-time boundary element method for the heat equation. We exploit the special temporal structure of the involved operators to provide an…
This article describes two independent developments aimed at improving the Particle Tracking Method for measurements of flow or particle velocities. First, a stereoscopic multicamera calibration method that does not require any optical…
A new algorithm is proposed to describe the propagation of fronts advected in the normal direction with prescribed speed function F. The assumptions on F are that it does not depend on the front itself, but can depend on space and time.…
Dissipative particle dynamics (DPD) is an effective mesoscopic particle model with a lower computational cost than molecular dynamics because of the soft potentials that it employs. However, the soft potential is not strong enough to…
We present the applications of methods from nonlinear local harmonic analysis in variational framework to calculations of nonlinear motions in polynomial/rational approximations (up to any order) of arbitrary n-pole fields. Our approach is…
We review three different approaches for the calculation of electromagnetic multipoles, namely the Cartesian primitive multipoles, the Cartesian irreducible multipoles and the spherical multipoles. We identify the latter as the best suited…
The accurate and efficient evaluation of Newtonian potentials over general 2-D domains is important for the numerical solution of Poisson's equation and volume integral equations. In this paper, we present a simple and efficient high-order…
A new method of numerical solution for partial differential equations is proposed. The method is based on a fast matrix multiplication algorithm. Two-dimensional Poison equation is used for comparison of the proposed method with…
Nonnegative tensors arise very naturally in many applications that involve large and complex data flows. Due to the relatively small requirement in terms of memory storage and number of operations per step, the (shifted) higher-order power…
Spectral line profiles are powerful diagnostic tools for both laboratory and astrophysical plasmas, as their shape is sensitive to the plasma environment. The low-frequency component of the electric microfield is an important input for…
We present an MPI-parallel algorithm for the in-situ visualization of computational data that is built around a distributed linear forest-of-octrees data structure. Such octrees are frequently used in element-based numerical simulations;…
Tree amplitudes of the production of two kinds of scalar particles at threshold from one virtual particle are calculated in a model of two scalar fields with $O(2)$ symmetric quartic interaction and unequal masses. These amplitudes exhibit…
This paper presents an octree construction method, called Cornerstone, that facilitates global domain decomposition and interactions between particles in mesh-free numerical simulations. Our method is based on algorithms developed for 3D…
Diffusion models have become a leading method for generative modeling of both image and scientific data. As these models are costly to train and \emph{evaluate}, reducing the inference cost for diffusion models remains a major goal.…
We present methods for generating computationally simple parameter-free pair potentials useful for solids, liquids and plasma at arbitrary temperatures. They successfully treat warm-dense matter (WDM) systems like carbon or silicon with…
We consider space-saving versions of several important operations on univariate polynomials, namely power series inversion and division, division with remainder, multi-point evaluation, and interpolation. Now-classical results show that…
A new order-N method for calculating the electronic structure of general (non-tight-binding) potentials is presented. The method uses a combination of the ``purification''-based approaches used by Li, Nunes and Vanderbilt, and Daw, and a…
In this paper, we propose a successive pseudo-convex approximation algorithm to efficiently compute stationary points for a large class of possibly nonconvex optimization problems. The stationary points are obtained by solving a sequence of…