User Manual Privacy Policy Disclaimer Contact us
  Advanced SearchBrowse




Journal Article

An optimizing symbolic algebra approach for generating fast multipole method operators


Bieri,  Rebekka
Galaxy Formation, MPI for Astrophysics, Max Planck Society;

External Ressource
No external resources are shared
Fulltext (public)
There are no public fulltexts stored in PuRe
Supplementary Material (public)
There is no public supplementary material available

Coles, J. P., & Bieri, R. (2020). An optimizing symbolic algebra approach for generating fast multipole method operators. Computer Physics Communications, 251: 107081. doi:10.1016/j.cpc.2019.107081.

Cite as: http://hdl.handle.net/21.11116/0000-0006-8FCB-5
We have developed a symbolic algebra approach to automatically produce, verify, and optimize computer code for the Fast Multipole Method (FMM) operators. This approach allows for flexibility in choosing a basis set and kernel, and can generate computer code for any expansion order in multiple languages. The procedure is implemented in the publicly available Python program Mosaic. Optimizations performed at the symbolic level through algebraic manipulations significantly reduce the number of mathematical operations compared with a straightforward implementation of the equations. We find that the optimizer is able to eliminate 20-80% of the floating-point operations and for the expansion orders p≤10 it changes the observed scaling properties. We present our approach using three variants of the operators with the Cartesian basis set for the harmonic potential kernel 1/r, including the use of totally symmetric and traceless multipole tensors.