On Fast Tabulation of Green’s Function in Layered Medium Based on the GPU Parallelism

The method of moments for layered media structures requires the establishment of a three-dimensional interpolation table to compute the Green's function of layered media using the interpolation method. For multi-layered large-scale structures, the computation of interpolation tables is very time-consuming. This paper proposes an efficient computation method for three-dimensional interpolation table of dyadic Green’s Function for layered medium based on the graphics processing units. The decoupling of the horizontal and vertical position parameters in the spectral domain Green's function and Bessel function is utilized to reduce the computational burden associated with the integrand in the direct integration method. The numerical results show that the proposed method has good accuracy and reduces the time for the creation of the interpolation table to a few seconds from dozens of minutes.