Layered ${\cal H}$ -Matrix Based Inverse and LU Algorithms for Fast Direct Finite-Element-Based Computation of Electromagnetic Problems

A layered -matrix based inverse algorithm and a layered -matrix based LU factorization and solution algorithm are developed to accelerate the direct solution of the finite element matrix resulting from electromagnetics-based analysis of general 3-D problems. In these algorithms, the direct solution of the original 3-D system matrix is transformed to the direct solution of multiple 2-D problems. The size of the matrix to be computed is thus reduced from a 3-D size to a 2-D size. Moreover, the growth rate of the rank of the inverse finite element matrix with electric size is reduced from a 3-D based growth rate to a 2-D based growth rate. Numerical experiments have demonstrated the clear advantages of the proposed direct solvers over a state-of-the-art multifrontal based direct sparse solver as well as existing fast H-matrix based direct solvers, in both CPU time and memory consumption, for finite-element analysis of general 3-D electromagnetic problems.

[1]  Robert J. Adams,et al.  Sparse factorization of finite element matrices using overlapped localizing solution modes , 2008 .

[2]  A. George.  Nested Dissection of a Regular Finite Element Mesh , 1973 .

[3]  W. Hackbusch.  A Sparse Matrix Arithmetic Based on $\Cal H$-Matrices. Part I: Introduction to ${\Cal H}$-Matrices , 1999, Computing.

[4]  Changhong Dai,et al.  A Layered Finite Element Method for Electromagnetic Analysis of Large-Scale High-Frequency Integrated Circuits , 2007, IEEE Transactions on Antennas and Propagation.

[5]  Boris N. Khoromskij,et al.  A Sparse H-Matrix Arithmetic. Part II: Application to Multi-Dimensional Problems , 2000, Computing.

[6]  Wolfgang Hackbusch,et al.  Construction and Arithmetics of H-Matrices , 2003, Computing.

[7]  Dan Jiao,et al.  A theoretical study on the rank's dependence with electric size of the inverse finite element matrix for large-scale electrodynamic analysis , 2012, Proceedings of the 2012 IEEE International Symposium on Antennas and Propagation.

[8]  Dan Jiao,et al.  A direct finite-element-based solver of significantly reduced complexity for solving large-scale electromagnetic problems , 2009, 2009 IEEE MTT-S International Microwave Symposium Digest.

[9]  W. Hackbusch,et al.  Introduction to Hierarchical Matrices with Applications , 2003 .

[10]  D. Jiao,et al.  Existence of ${\cal H}$-Matrix Representations of the Inverse Finite-Element Matrix of Electrodynamic Problems and ${\cal H}$-Based Fast Direct Finite-Element Solvers , 2010, IEEE Transactions on Microwave Theory and Techniques.