ИСТИНА

New numerical method for 3d integral equations of electrodynamics was devised. The most challenging aspect of using these equations in forward problems of EM sounding is associated with a need to solve a high-order system of linear equations with a dense matrix. The main distinctive feature of the proposed method is a double integration of the electric Green's tensor in the process of algebraization of the original equation to decrease the condition number of the matrix. Another computational challenge arises from the fact that the Green's tensors elements are integrals containing the Bessel functions. In order to solve this problem we developed the new method motivated with W. Anderson's ideas on Hankel transforms calculations. The method of integral averages was implemented as parallel program using both MPI and OpenMP technologies. The computational experiments were performed by Bluegene/P and “Lomonosov” supercomputers from MSU. These experiments showed the efficiency of the proposed method in high-contrast electromagnetic sounding problems.