English

A fast multipole method for stellar dynamics

Instrumentation and Methods for Astrophysics 2014-05-12 v1 Computational Physics

Abstract

The approximate computation of all gravitational forces between NN interacting particles via the fast multipole method (FMM) can be made as accurate as direct summation, but requires less than O(N)\mathcal{O}(N) operations. FMM groups particles into spatially bounded cells and uses cell-cell interactions to approximate the force at any position within the sink cell by a Taylor expansion obtained from the multipole expansion of the source cell. By employing a novel estimate for the errors incurred in this process, I minimise the computational effort required for a given accuracy and obtain a well-behaved distribution of force errors. For relative force errors of 107\sim10^{-7}, the computational costs exhibit an empirical scaling of N0.87\propto N^{0.87}. My implementation (running on a 16 core node) out-performs a GPU-based direct summation with comparable force errors for N105N\gtrsim10^5.

Keywords

Cite

@article{arxiv.1405.2255,
  title  = {A fast multipole method for stellar dynamics},
  author = {Walter Dehnen},
  journal= {arXiv preprint arXiv:1405.2255},
  year   = {2014}
}

Comments

21 pages, 15 figures, accepted for publication in Journal for Computational Astrophysics and Cosmology

R2 v1 2026-06-22T04:10:10.225Z