Reduced forward operator for electromagnetic wave scattering problems

The paper describes a reduced forward operator for solving electromagnetic scattering problems using a volume integral equation in conjunction with a conjugate gradient fast Fourier transform scheme. The reduction is obtained by decoupling of the interaction between the locations in the spatial computational domain at which there is non-zero contrast and those positions at which there is zero contrast. The decoupling is achieved by multiplication of the kernel by a diagonal matrix whose entries reflect the presence or absence of contrast at the associated point. Analysis supported by numerical experiments shows that the conjugate gradient algorithm applied to the reduced system converges more rapidly than when it is applied to the original system