RECFMM: Recursive Parallelization of the Adaptive Fast Multipole Method for Coulomb and Screened Coulomb Interactions


Journal Article

Copyright © Global-Science Press 2016. We present RECFMM, a program representation and implementation of a recursive scheme for parallelizing the adaptive fast multipole method (FMM) on shared-memory computers. It achieves remarkable high performance while maintaining mathematical clarity and flexibility. The parallelization scheme signifies the recursion feature that is intrinsic to the FMM but was not well exploited. The program modules of RECFMM constitute a map between numerical computation components and advanced architecture mechanisms. The mathematical structure is preserved and exploited, not obscured nor compromised, by parallel rendition of the recursion scheme. Modern software system - CILK in particular, which provides graph-theoretic optimal scheduling in adaptation to the dynamics in parallel execution - is employed. RECFMM supports multiple algorithm variants that mark the major advances with low-frequency interaction kernels, and includes the asymmetrical version where the source particle ensemble is not necessarily the same as the target particle ensemble. We demonstrate parallel performance with Coulomb and screened Coulomb interactions.

Full Text

Duke Authors

Cited Authors

  • Zhang, B; Huang, J; Pitsianis, NP; Sun, X

Published Date

  • August 1, 2016

Published In

Volume / Issue

  • 20 / 2

Start / End Page

  • 534 - 550

Electronic International Standard Serial Number (EISSN)

  • 1991-7120

International Standard Serial Number (ISSN)

  • 1815-2406

Digital Object Identifier (DOI)

  • 10.4208/cicp.230216.140416sw

Citation Source

  • Scopus