Efficient fast multipole method for low-frequency scattering (Q598150): Difference between revisions

The solution of the Helmholtz and Maxwell equations using integral formulations requires to solve large complex linear systems. A direct solution of those problems using a Gauss elimination is practical only for very small systems with few unknowns. The use of an iterative method such as GMRES can reduce the computational expense. Most of the expense is then computing large complex matrix vector products. The cost can be further reduced by using the fast multipole method which accelerates the matrix vector product. For a linear system of size \(N\), the use of an iterative method combined with the fast multipole method reduces the total expense of the computation to \(N\log N\). There exist two versions of the fast multipole method: one which is based on a multipole expansion of the interaction kernel \(\exp\iota kr/r\) and which was first proposed by \textit{V. Rokhlin} [ibid. 60, 187--207 (1985; Zbl 0629.65122)] and another based on a plane wave expansion of the kernel, first proposed by \textit{W. C. Chew, J. M. Jin, C. C. Lu, E. Michielssen} and \textit{J. M. M. Song} [Fast solution emthods in electromagnetics, IEEE Trans. Antenn. Progag. 45, No. 3, 533--543 (1997)]. In this paper, the authors propose a third approach, the stable plane wave expansion, which has a lower computational expense than the multipole expansion and does not have the accuracy and stability problems of the plane wave expansion. The computational complexity is \(N\log N\) as with the other methods.