The fast multipole method for gridless particle simulation

John Ambrosiano, Leslie Greengard, Vladimir Rokhlin

Research output: Contribution to journalArticle

Abstract

The numerical solution to N-body problems in gravitation or electrostatics has traditionally been obtained via particle-in-cell methods (PIC) since direct evaluation of all pairwise interparticle forces, requiring O(N2) operations, is too expensive. Recently, hierarchical solvers, which use tree data structures and lumped-force approximations, have made gridless simulations feasible in O(N log N) operations. In this paper, we explore the use of the fast multipole method (FMM) - a highly accurate order O(N) algorithm - in particle simulations. We briefly describe the FMM and its relation to other methods. Technical considerations of gridless simulations such as discrete particle fluctuations, sampling errors and boundary conditions are discussed and compared with PIC methodology. Examples of electrostatic simulations in plasma physics are presented.

Original languageEnglish (US)
Pages (from-to)117-125
Number of pages9
JournalComputer Physics Communications
Volume48
Issue number1
DOIs
StatePublished - 1988

Fingerprint

multipoles
Electrostatics
Data structures
Gravitation
Physics
simulation
Boundary conditions
electrostatics
Sampling
Plasmas
many body problem
data structures
plasma physics
cells
sampling
methodology
boundary conditions
gravitation
evaluation
approximation

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy(all)

Cite this

The fast multipole method for gridless particle simulation. / Ambrosiano, John; Greengard, Leslie; Rokhlin, Vladimir.

In: Computer Physics Communications, Vol. 48, No. 1, 1988, p. 117-125.

Research output: Contribution to journalArticle

Ambrosiano, John ; Greengard, Leslie ; Rokhlin, Vladimir. / The fast multipole method for gridless particle simulation. In: Computer Physics Communications. 1988 ; Vol. 48, No. 1. pp. 117-125.
@article{c08a508593ff468d82e757fee6749a0b,
title = "The fast multipole method for gridless particle simulation",
abstract = "The numerical solution to N-body problems in gravitation or electrostatics has traditionally been obtained via particle-in-cell methods (PIC) since direct evaluation of all pairwise interparticle forces, requiring O(N2) operations, is too expensive. Recently, hierarchical solvers, which use tree data structures and lumped-force approximations, have made gridless simulations feasible in O(N log N) operations. In this paper, we explore the use of the fast multipole method (FMM) - a highly accurate order O(N) algorithm - in particle simulations. We briefly describe the FMM and its relation to other methods. Technical considerations of gridless simulations such as discrete particle fluctuations, sampling errors and boundary conditions are discussed and compared with PIC methodology. Examples of electrostatic simulations in plasma physics are presented.",
author = "John Ambrosiano and Leslie Greengard and Vladimir Rokhlin",
year = "1988",
doi = "10.1016/0010-4655(88)90029-X",
language = "English (US)",
volume = "48",
pages = "117--125",
journal = "Computer Physics Communications",
issn = "0010-4655",
publisher = "Elsevier",
number = "1",

}

TY - JOUR

T1 - The fast multipole method for gridless particle simulation

AU - Ambrosiano, John

AU - Greengard, Leslie

AU - Rokhlin, Vladimir

PY - 1988

Y1 - 1988

N2 - The numerical solution to N-body problems in gravitation or electrostatics has traditionally been obtained via particle-in-cell methods (PIC) since direct evaluation of all pairwise interparticle forces, requiring O(N2) operations, is too expensive. Recently, hierarchical solvers, which use tree data structures and lumped-force approximations, have made gridless simulations feasible in O(N log N) operations. In this paper, we explore the use of the fast multipole method (FMM) - a highly accurate order O(N) algorithm - in particle simulations. We briefly describe the FMM and its relation to other methods. Technical considerations of gridless simulations such as discrete particle fluctuations, sampling errors and boundary conditions are discussed and compared with PIC methodology. Examples of electrostatic simulations in plasma physics are presented.

AB - The numerical solution to N-body problems in gravitation or electrostatics has traditionally been obtained via particle-in-cell methods (PIC) since direct evaluation of all pairwise interparticle forces, requiring O(N2) operations, is too expensive. Recently, hierarchical solvers, which use tree data structures and lumped-force approximations, have made gridless simulations feasible in O(N log N) operations. In this paper, we explore the use of the fast multipole method (FMM) - a highly accurate order O(N) algorithm - in particle simulations. We briefly describe the FMM and its relation to other methods. Technical considerations of gridless simulations such as discrete particle fluctuations, sampling errors and boundary conditions are discussed and compared with PIC methodology. Examples of electrostatic simulations in plasma physics are presented.

UR - http://www.scopus.com/inward/record.url?scp=0023166090&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0023166090&partnerID=8YFLogxK

U2 - 10.1016/0010-4655(88)90029-X

DO - 10.1016/0010-4655(88)90029-X

M3 - Article

AN - SCOPUS:0023166090

VL - 48

SP - 117

EP - 125

JO - Computer Physics Communications

JF - Computer Physics Communications

SN - 0010-4655

IS - 1

ER -