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

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.