An H-matrix direct method with improved solution for finite element analysis of scattering problems

An efficient direct method based on hierarchical matrices (H-matrices) was developed for the solution of large sparse finite element systems. H-matrix technique provides a data-sparse way to approximate the inverse of a finite element matrix which is dense originally. The approximate inverse can be computed and stored in the algebra of hierarchical matrices with almost linear complexity. The accuracy of this approximation is controllable with different choices of the related parameters. A modification algorithm was introduced to iteratively improve the accuracy of an approximate solution generated from H-matrix inversion. Some numerical examples are provided to illustrate the accuracy and efficiency of this method.

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

[2]  J. Pasciak,et al.  Computer solution of large sparse positive definite systems , 1982 .

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

[4]  Wolfgang Hackbusch,et al.  A Sparse Matrix Arithmetic Based on H-Matrices. Part I: Introduction to H-Matrices , 1999, Computing.

[5]  Carretera de Valencia,et al.  The finite element method in electromagnetics , 2000 .

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

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

[8]  E. Yung,et al.  SSOR preconditioned GMRES for the FEM analysis of waveguide discontinuities with anisotropic dielectric , 2004 .

[9]  Sabine Le Borne,et al.  H-matrix Preconditioners in Convection-Dominated Problems , 2005, SIAM J. Matrix Anal. Appl..

[10]  Z. Nie,et al.  Application of diagonally perturbed incomplete factorization preconditioned conjugate gradient algorithms for edge finite-element analysis of Helmholtz equations , 2006, IEEE Transactions on Antennas and Propagation.

[11]  Mario Bebendorf,et al.  Why Finite Element Discretizations Can Be Factored by Triangular Hierarchical Matrices , 2007, SIAM J. Numer. Anal..

[12]  D. Jiao,et al.  Performance analysis of the H-matrix-based Fast Direct Solver for finite-element-based analysis of electromagnetic problems , 2009, 2009 IEEE Antennas and Propagation Society International Symposium.

[13]  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.