Adapting free-space fast multipole method for layered media Green's function: Algorithm and analysis

We present a numerical algorithm for an accurate and efficient computation of the convolution of the frequency domain layered media Green's function with a given density function. Instead of compressing the convolution matrix directly as in the classical fast multipole, fast direct solver, and H-matrix algorithms, this new algorithm considers a translated form of the original matrix so existing blocks from the highly optimized free-space fast multipole method can be easily adapted to the layered media Green's function. An asymptotic analysis is performed on the Sommerfeld integrals to provide an estimate of the decay rate in the new “multipole” and “local” expansions. To avoid the highly oscillatory integrand in the original integral representations when the source and target are close to each other, mathematically equivalent alternative direction integral representations are introduced. The convergence of the new expansions and quadrature rules for the original and alternative direction representations are numerically validated.

Cho, M. H., & Huang, J. (2021). Adapting free-space fast multipole method for layered media Green’s function: Algorithm and analysis. Applied and Computational Harmonic Analysis, 51, 414–436. doi:10.1016/j.acha.2019.10.001

J. Huang was supported by the NSF grant DMS1821093, and the work was finished while he was visiting professors at the King Abdullah University of Science and Technology, National Center for Theoretical Sciences (NCTS) in Taiwan, Mathematical Center for Interdisciplinary Research of Soochow University, and Institute for Mathematical Sciences of the National University of Singapore. M.H. Cho was supported by a grant from the Simons Foundation (No. 404499).

Elsevier BV

Applied and Computational Harmonic Analysis



Additional Links

Permanent link to this record

Version History

Now showing 1 - 2 of 2
2022-06-12 12:57:16
Published in journal
2019-12-29 12:05:00
* Selected version