Matrix equilibration in method of moment solutions of surface integral equations

Abstract

Basic theory of matrix equilibration is presented, relating it to other techniques for decreasing the condition number of matrix equations obtained by the method of moments (MOM) applied to surface integral equations (SIEs). It is shown that matrix equilibration is a general technique that can be used for both (1) balancing field and source quantities in SIEs, which is used to decrease the condition number in the case of SIEs of mixed type and high contrast in material properties, and (2) scaling basis and test functions in MOM, which is used to decrease the condition number in the case of higher-order bases and patches of different sizes. In particular, it is demonstrated that a combination of such balancing and scaling can be performed using simple matrix equilibration based on magnitudes of diagonal elements and 2-norms of rows/columns of the MOM matrix.

On Sparsity Related to the Integral Equations of 2D Electromagnetic Scattering

“…Major advances in the numerical solution of the integral equations of scattering have been made with the Fast Multipole Method (Chew et al, 2001;Hansen, 2013Hansen, , 2014Kolundzija & Kostic, 2014;Rokhlin, 1990;Ylä-Oijala et al, 2005). The progress mainly relates to the efficient handling of many basis functions and large systems of equations.…”

Fast Integral Equation Solution of Scattering of Multiscale Objects by Domain Decomposition Method with Mixed Basis Functions

“…With the HOHV basis function, a well-conditioned matrix can be obtained for conventional structures [9]. However, when modeling multiscale structures, the high contrast (electrically large and electrically small) sized meshes exist together, so the proper equilibration technique of MoM matrix equation is needed to decrease the condition number 2 International Journal of Antennas and Propagation [11]. At the same time, ill-shaped meshes are inevitable and the orthogonality of HOHV basis function is affected.…”