Efficient computation of vertical current reactions in layered-media for the mixed-potential integral equation

For method of moments simulations of printed structures in layered media that involve both vertical and horizontal metallization, a two-dimensional interpolation scheme is generally required to evaluate efficiently the Sommerfeld integrals arising from the Green functions. A more effective way is to integrate the vertical bounds of the reaction integrals analytically in the spectral-domain. This is further accelerated by employing the discrete complex image method (DCIM) to approximate the radial variation, thus avoiding an expensive Sommerfeld integral.