Scattering from quasi-planar and moderate rough surfaces: Efficient method to fill the EFIE-Galerkin MoM impedance matrix and to solve the linear system

International audienceFirst, an acceleration to compute the impedance matrix obtained from the EFIE (electric field integral equation) discretized by the Galerkin MoM (method of moments) with Rao-Wilton-Glisson basis functions, is addressed. It is based on a far-field approximation and makes it possible to avoid to loop on the source and observation triangles. Next, the impedance matrix is split into strong and weak interactions; the latter is compressed by expressing it from Toeplitz sub-matrices. Then, the linear system is efficiently solved by a bi-iterative scheme. For a given order, the LSQR (least squares QR) algorithm is applied to solve the sparse linear system related to the strong interactions, while the matrix-vector products, related to the weak interactions, are accelerated by using FFTs. Numerical results of the field scattered by perfectly-conducting paraboloid-shape object and Gaussian rough surface are shown