Singularity subtraction for evaluation of Green´s functions for multilayer media

This paper presents an efficient method to evaluate the two- and three-dimensional multilayered medium Green´s functions for general electric and magnetic sources. Without finding any surface poles or steepest descent path, a special subtraction procedure is applied to each term of the Sommerfeld integrands to make them rapidly decreasing functions of k_rho. The contributions of the subtracted terms are calculated analytically. The remaining integrals are computed adaptively by using Gaussian quadratures. The accuracy of the method has been confirmed by comparison with many examples in literature, and the high efficiency has been verified