Layered H-Matrix Based Inverse and LU Algorithms for Fast Direct Finite-Element-Based Computation of Electromagnetic Problems

作者：Liu Haixin^*; Jiao Dan

来源：IEEE Transactions on Antennas and Propagation, 2013, 61(3): 1273-1284.

DOI：10.1109/TAP.2012.2230236

• H-matrix
• direct matrix solution
• finite element methods
• electromagnetic analysis

A layered H-matrix based inverse algorithm and a layered H-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 systemmatrix 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.

• 出版日期2013-3