Reza Gholami, Jamiu Mojolagbe, Anton Menshov, Farhad Sheikh Hosseini Lori, Vladimir Okhmatovski
R. Gholami, J. Mojolagbe, A. Menshov, F. Sheikh Hosseini Lori, and V. Okhmatovski, “An H‑matrix accelerated solution of a new single-source integral equation for scattering on penetrable objects,” in USNC-URSI Radio Sci. Meet., San-Diego, CA, Jul. 2017.
Publication year: 2017


Recently, a novel single-source integral equation (SSIE) has been proposed for 2D (A. Menshov,, IEEE T-MTT, no. 1, vol. 61, pp. 341–350, 2013) and 3D (F. Sheikh,, submitted to IEEE T-APS, 2016) scattering problems on homogeneous dielectric objects. Unlike the previously known SSIEs, the new equation involves only one product of integral operators and requires only electric-field-type of Green’s function while still reducing the number of unknowns by half compared to traditional surface integral equation formulations. However, memory requirements and computational complexity of dense matrix operations and storage for naive MoM solution of the novel SSIE limit its application to small problems. In this work, we investigate a fast direct and iterative solutions for the new SSIE based on the hierarchical matrix (H-matrix) framework (W. Hackbusch, Computing, 62, 89–108, 1999).

At first, H-matrix based solver hierarchically divides the structure of interest into sub-domains such that the matrix is decomposed into sub-blocks each representing the interaction between a pair of sub-domains. Then, for each well-separated pair of sub-domains τ and σ, the corresponding block Zτ×σ is expected to be low rank and is represented in the Zτ×σ ≈ AB^H form. The accuracy of the approximation is dependent on the low-rank factorization strategy, such as singular value decomposition, adaptive cross approximation, and grid-based algorithms. For the new SSIE, each of the entering integral operators (surface-to-surface, surface-to-volume, and volume-to-surface) is approximated using H-matrices, and their combination is performed in the already compressed form. Acceleration of the MoM solution by H-matrices results in the significant reduction of storage requirements and asymptotic complexity of matrix-matrix products, LU-decomposition, and back-substitution.

At the conference, the performance of the algorithm in terms of the accuracy, computational time, and memory will be discussed and finally the optimized method will be applied to accelerate the solution of the new SSIE for the large-scale 3D scattering problem.