The recursive Green's function method for surface integral equation analysis of inhomogeneous media

The numerical simulation of electromagnetic phenomena in continuously varying inhomogeneous domains is an important problem which finds application in numerous areas. Many such simulations involve a moment method solution of the volume integral equation. The difficulty with this approach is the overwhelming computational resources consumed when applied to even moderately-sized domains. This paper introduces a new approach to solving the scalar wave equation in inhomogeneous domains which allows substantial reduction in the computational storage and complexity requirements as compared to conventional techniques. This algorithm, called the recursive Green's function method (Jensen and Freeze, 1996), allows efficient construction of the Green's function for the region which can subsequently be used in surface integral formulations for the fields.