The three-dimensional algorithm of solving the electric field integral equation using face-centered node points, conjugate gradient method, and FFT

It has been known for a long time that the accuracy of solving the scattering by a dielectric body using the electric field integral equation (EFIE) is poor when the permittivity of the scatterer becomes large. Recently, this problem has been settled by using a procedure involving face-centered node points. Such a procedure is efficient, since it preserves the convolution property in the EFIE and, hence, the applicability of the fast Fourier transform (FFT). This procedure is generalized to the three-dimensional and anisotropic case. The generalization is quite straightforward in both the formulation and the programming. A calculation for a scatterer with a relative permittivity as high as 100 indicates that the proposed procedure converges quite rapidly, whereas the conventional using the conjugate gradient method approach fails to converge.>