Preprint 62/2004

H-matrix preconditioners in convection-dominated problems

Sabine Le Borne and Lars Grasedyck

Contact the author: Please use for correspondence this email.
Submission date: 28. Sep. 2004
Pages: 11
published in: SIAM journal on matrix analysis and applications, 27 (2006) 4, p. 1172-1183 
DOI number (of the published article): 10.1137/040615845
MSC-Numbers: 65F05, 65F30, 65F50
Keywords and phrases: hierarchical matrices, preconditioning, convection-dominant problems
Download full preprint: PDF (1662 kB), PS ziped (2385 kB)

Hierarchical matrices provide a data-sparse way to approximate fully populated matrices. In this paper we exploit formula8-matrix techniques to approximate the LU-decompositions of stiffness matrices as they appear in (finite element or finite difference) discretizations of convection-dominated elliptic partial differential equations. These sparse formula8-matrix approximations may then be used as preconditioners in iterative methods. Whereas the approximation of the matrix inverse by an formula8-matrix requires some modification in the underlying index clustering when applied to convection-dominant problems, the formula8-LU-decomposition works well in the standard formula8-matrix setting even in the convection dominant case. We will complement our theoretical analysis with some numerical examples.

03.07.2017, 01:41