RECFMM: recursive parallelization of the adaptive fast multipole method for Coulomb and screened Coulomb interactions. (English) Zbl 1373.78001

Zhang, Bo; Huang, Jingfang; Pitsianis, Nikos P.; Sun, Xiaobai

Commun. Comput. Phys. 20, No. 2, 534-550 (2016).

Summary: 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 systems – CILK in particular, which provides graph-theoretic optimal scheduling in adaptation to the dynamics in parallel execution – are 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.

MSC:

78-04	Software, source code, etc. for problems pertaining to optics and electromagnetic theory
68W10	Parallel algorithms in computer science

Software:

RECFMM; PAFMPB; Cilk