Broadband multilevel fast multipole algorithm for acoustic scattering problems

A broadband multilevel fast multipole algorithm (MLFMA) for the acoustic scattering from a sound-hard obstacle is presented. The formulation is based on the Burton-Miller boundary integral equation and Galerkin's method, avoiding any hypersingular integral operators. The resulting matrix equation has good iterative properties for all frequencies and avoids the interior resonance problem. The main novel feature is the use of a broadband MLFMA to accelerate the iterative generalized minimal residual (GMRES) solver. The algorithm is based on a combination of Rokhlin's translation formula for large division cubes and the spectral representation of the Green's function for cubes smaller than one half wavelength, thereby avoiding the sub-wavelength breakdown of the high-frequency MLFMA.