An optimizing symbolic algebra approach for generating fast multipole method 
operators

Coles, Jonathan P.; Bieri, Rebekka

doi:10.1016/j.cpc.2019.107081

Item

ITEM ACTIONSEXPORT

Add to Basket

Local TagsRelease HistoryDetailsSummary

Released

Journal Article

An optimizing symbolic algebra approach for generating fast multipole method operators

MPS-Authors

/persons/resource/persons222724

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

External Resource

No external resources are shared

Fulltext (restricted access)

There are currently no full texts shared for your IP range.

Fulltext (public)

There are no public fulltexts stored in PuRe

Supplementary Material (public)

There is no public supplementary material available

Citation

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: https://hdl.handle.net/21.11116/0000-0006-8FCB-5

Abstract

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.