Two-step lagrange interpolation method for the multilevel fast multipole algorithm

We present a two-step Lagrange interpolation method for the efficient solution of large-scale electromagnetics problems with the multilevel fast multipole algorithm (MLFMA). Local interpolations are required during aggregation and disaggregation stages of MLFMA in order to match the different sampling rates for the radiated and incoming fields in consecutive levels. The conventional one-step method is decomposed into two one-dimensional interpolations, applied successively. As it provides a significant acceleration in processing time, the proposed two-step method is especially useful for problems involving large-scale objects discretized with millions of unknowns.