A fast multipole method for Maxwell equations stable at all frequencies.

The solution of Helmholtz and Maxwell equations by integral formulations (kernel in exp(i kr)/r) leads to large dense linear systems. Using direct solvers requires large computational costs in O(N(3)). Using iterative solvers, the computational cost is reduced to large matrix-vector products.