Summation of certain locally bilinear forms and its applications to the Fast Multipole Method

TL;DR

This paper introduces a simplified implementation method for the Fast Multipole Method (FMM) that eliminates complex shift processes, making coding easier and reducing computation time while maintaining accuracy.

Contribution

The paper presents a novel approach to implement FMM without shift processes, simplifying coding and improving computational efficiency.

Findings

Simplified FMM implementation reduces coding complexity.

Elimination of shift processes saves computational time.

Accuracy of potential field calculations is maintained.

Abstract

The Fast Multipole Method (FMM) reduces the computation of pairwise two-body interactions among $N$-particles to order $N$, whose computation cost should be of order $N^2$ by brute force. However, its implementation is somewhat complicated and requires a considerable amount of time to write the code. In this paper, I show a method that enables us to implement and write FMM algorithm code simply and briefly. FMM algorithm is composed of several steps. The main steps are Upward Pass and Downward Pass. Both the Upward Pass and Downward Pass include shift processes by which we move the centers of local expansions and multipole expansions. In this paper, I show a method that enables us to get rid of these processes.As a result of this simplification, the coding of FMM becomes much easier, and we can save considerable computation time. I compared the accuracy and time required to calculate potential fields with that of the existing FMM code.