Numerical accuracy of multipole expansion for 2-D MLFMA

Numerical study of the multipole expansion for the multilevel fast multipole algorithm (MLFMA) is presented. In the numerical implementation of MLFMA, the error comes from three sources: the truncation of the addition theorem; the approximation of the integration; and the aggregation and disaggregation process. These errors are due to the factorization of the Green's function which is the mathematical core of this algorithm. Among the three error sources, we focus on the truncation error in this paper and a new approach of selecting truncation numbers for the addition theorem is proposed. Using this approach, the error prediction and control can be improved for the small buffer sizes and high accuracy requirements.