A kernel-independent adaptive fast multipole algorithm in two and three dimensions

Lexing Ying, George Biros, Denis Zorin

Research output: Contribution to journalArticle

Abstract

We present a new fast multipole method for particle simulations. The main feature of our algorithm is that it does not require the implementation of multipole expansions of the underlying kernel, and it is based only on kernel evaluations. Instead of using analytic expansions to represent the potential generated by sources inside a box of the hierarchical FMM tree, we use a continuous distribution of an equivalent density on a surface enclosing the box. To find this equivalent density, we match its potential to the potential of the original sources at a surface, in the far field, by solving local Dirichlet-type boundary value problems. The far-field evaluations are sparsified with singular value decomposition in 2D or fast Fourier transforms in 3D. We have tested the new method on the single and double layer operators for the Laplacian, the modified Laplacian, the Stokes, the modified Stokes, the Navier, and the modified Navier operators in two and three dimensions. Our numerical results indicate that our method compares very well with the best known implementations of the analytic FMM method for both the Laplacian and modified Laplacian kernels. Its advantage is the (relative) simplicity of the implementation and its immediate extension to more general kernels.

Original languageEnglish (US)
Pages (from-to)591-626
Number of pages36
JournalJournal of Computational Physics
Volume196
Issue number2
DOIs
StatePublished - May 20 2004

Fingerprint

multipoles
boxes
far fields
Singular value decomposition
Fast Fourier transforms
Boundary value problems
operators
expansion
evaluation
boundary value problems
decomposition
simulation

Keywords

  • Double-layer potential
  • Fast multipole methods
  • Fast solvers
  • Integral equations
  • N-body problems
  • Particle methods
  • Single-layer potential

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy(all)

Cite this

A kernel-independent adaptive fast multipole algorithm in two and three dimensions. / Ying, Lexing; Biros, George; Zorin, Denis.

In: Journal of Computational Physics, Vol. 196, No. 2, 20.05.2004, p. 591-626.

Research output: Contribution to journalArticle

@article{2efce92aa9ba42bd813b92f2659c4446,
title = "A kernel-independent adaptive fast multipole algorithm in two and three dimensions",
abstract = "We present a new fast multipole method for particle simulations. The main feature of our algorithm is that it does not require the implementation of multipole expansions of the underlying kernel, and it is based only on kernel evaluations. Instead of using analytic expansions to represent the potential generated by sources inside a box of the hierarchical FMM tree, we use a continuous distribution of an equivalent density on a surface enclosing the box. To find this equivalent density, we match its potential to the potential of the original sources at a surface, in the far field, by solving local Dirichlet-type boundary value problems. The far-field evaluations are sparsified with singular value decomposition in 2D or fast Fourier transforms in 3D. We have tested the new method on the single and double layer operators for the Laplacian, the modified Laplacian, the Stokes, the modified Stokes, the Navier, and the modified Navier operators in two and three dimensions. Our numerical results indicate that our method compares very well with the best known implementations of the analytic FMM method for both the Laplacian and modified Laplacian kernels. Its advantage is the (relative) simplicity of the implementation and its immediate extension to more general kernels.",
keywords = "Double-layer potential, Fast multipole methods, Fast solvers, Integral equations, N-body problems, Particle methods, Single-layer potential",
author = "Lexing Ying and George Biros and Denis Zorin",
year = "2004",
month = "5",
day = "20",
doi = "10.1016/j.jcp.2003.11.021",
language = "English (US)",
volume = "196",
pages = "591--626",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "2",

}

TY - JOUR

T1 - A kernel-independent adaptive fast multipole algorithm in two and three dimensions

AU - Ying, Lexing

AU - Biros, George

AU - Zorin, Denis

PY - 2004/5/20

Y1 - 2004/5/20

N2 - We present a new fast multipole method for particle simulations. The main feature of our algorithm is that it does not require the implementation of multipole expansions of the underlying kernel, and it is based only on kernel evaluations. Instead of using analytic expansions to represent the potential generated by sources inside a box of the hierarchical FMM tree, we use a continuous distribution of an equivalent density on a surface enclosing the box. To find this equivalent density, we match its potential to the potential of the original sources at a surface, in the far field, by solving local Dirichlet-type boundary value problems. The far-field evaluations are sparsified with singular value decomposition in 2D or fast Fourier transforms in 3D. We have tested the new method on the single and double layer operators for the Laplacian, the modified Laplacian, the Stokes, the modified Stokes, the Navier, and the modified Navier operators in two and three dimensions. Our numerical results indicate that our method compares very well with the best known implementations of the analytic FMM method for both the Laplacian and modified Laplacian kernels. Its advantage is the (relative) simplicity of the implementation and its immediate extension to more general kernels.

AB - We present a new fast multipole method for particle simulations. The main feature of our algorithm is that it does not require the implementation of multipole expansions of the underlying kernel, and it is based only on kernel evaluations. Instead of using analytic expansions to represent the potential generated by sources inside a box of the hierarchical FMM tree, we use a continuous distribution of an equivalent density on a surface enclosing the box. To find this equivalent density, we match its potential to the potential of the original sources at a surface, in the far field, by solving local Dirichlet-type boundary value problems. The far-field evaluations are sparsified with singular value decomposition in 2D or fast Fourier transforms in 3D. We have tested the new method on the single and double layer operators for the Laplacian, the modified Laplacian, the Stokes, the modified Stokes, the Navier, and the modified Navier operators in two and three dimensions. Our numerical results indicate that our method compares very well with the best known implementations of the analytic FMM method for both the Laplacian and modified Laplacian kernels. Its advantage is the (relative) simplicity of the implementation and its immediate extension to more general kernels.

KW - Double-layer potential

KW - Fast multipole methods

KW - Fast solvers

KW - Integral equations

KW - N-body problems

KW - Particle methods

KW - Single-layer potential

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

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

U2 - 10.1016/j.jcp.2003.11.021

DO - 10.1016/j.jcp.2003.11.021

M3 - Article

AN - SCOPUS:2442446356

VL - 196

SP - 591

EP - 626

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 2

ER -