Abstract:We demonstrate a new, hybrid symbolic-numerical method for the automatic synthesis of all families of translation operators required for the execution of the Fast Multipole Method (FMM). Our method is applicable in any dimensionality and to any translation-invariant kernel. The Fast Multipole Method, of course, is the leading approach for attaining linear complexity in the evaluation of long-range (e.g. Coulomb) many-body interactions. Low complexity in translation operators for the Fast Multipole Method (FMM) is usually achieved by algorithms specialized for a potential obeying a specific partial differential equation (PDE). Absent a PDE or specialized algorithms, Taylor series based FMMs or kernel-independent FMM have been used, at asymptotically higher expense. When symbolically provided with a constant-coefficient elliptic PDE obeyed by the potential, our algorithm can automatically synthesize translation operators requiring $\mathrm{O}(p^d)$ operations, where $p$ is the expansion order and $d$ is dimension, compared with $\mathrm{O}(p^{2d})$ operations in a naive approach carried out on (Cartesian) Taylor expansions. This is achieved by using a compression scheme that asymptotically reduces the number of terms in the Taylor expansion and then operating directly on this ``compressed'' representation. Judicious exploitation of shared subexpressions permits formation, translation, and evaluation of local and multipole expansions to be performed in $\mathrm{O}(p^{d})$ operations, while an FFT-based scheme permits multipole-to-local translations in $\mathrm{O}(p^{d-1}\log(p))$ operations. We demonstrate computational scaling of code generation and evaluation as well as numerical accuracy through numerical experiments on a number of potentials from classical physics.

One-dimension local interpolation algorithm for fast computation of translation operators in MLFMA

An Efficient Modified Interpolation Technique For The Translation Operators In Mlfma

Modified Lagrange Interpolation Technique for Translation Operators in Multilevel Fast Multipole Algorithm (MLFMA)

Efficient Interpolation Scheme for Multilevel Fast Multipole Algorithm

A Directional Interpolation Pattern for Multilevel Green's Function Interpolation Method

An Implementation of Parallel MLFMA on a Cluster of Computers with Distributed Memory

A Directional Multi-Level Complex-Space Fast Multipole Algorithm for Electrically Large Problems of Various Dimensions

Efficient Vector Interpolation Method Based on Lebedev Quadrature for MLFMA

Three-dimensional Local Multilevel Fast Multipole Algorithm

An Improved Implementation of Multilevel Green's Function Interpolation Method for Electric Field Integral Equation

A Local Multilevel Fast Multipole Algorithm for 3D Electromagnetic Scattering

Analysis of Multiscale Problems Using the MLFMA with the Assistance of the FFT-Based Method

Modified Multipole Number Technique in MLFMA

An Efficient Solution for the Parallelization of the Multilevel Fast Multipole Algorithm

A Parallel Multilevel Fast Multipole Algorithm Based on OpenMP

An Improved Full-Wave Multilevel Green’s Function Interpolation Method with RBF-QR Technique for Fast Field Evaluation

A New Preconditioning Scheme for MLFMA Based on Null-Field Generation Technique

Partly Approximate Iteration Technique For Accelerating Cg-Mlfma Solution

Solving 3d Electromagnetic Scattering And Radiation By Local Multilevel Fast Multipole Algorithm

Hybrid IE-DDM-MLFMA with Forward-Backward Method for Conducting Body of Translation

Automatic Synthesis of Low-Complexity Translation Operators for the Fast Multipole Method