Volume 26, pp. 285-298, 2007.

Joint domain-decomposition H-LU preconditioners for saddle point problems

Sabine Le Borne and Suely Oliveira

Abstract

For saddle point problems in fluid dynamics, several popular preconditioners exploit the block structure of the problem to construct block triangular preconditioners. The performance of such preconditioners depends on whether fast, approximate solvers for the linear systems on the block diagonal (representing convection-diffusion problems) as well as for the Schur complement (in the pressure variables) are available. In this paper, we will introduce a completely different approach in which we ignore this given block structure. We will instead compute an approximate LU-factorization of the complete system matrix using hierarchical matrix techniques. In particular, we will use domain-decomposition clustering with an additional local pivoting strategy to order the complete index set. As a result, we obtain an H-matrix structure in which an H-LU factorization is computed more efficiently and with higher accuracy than for the corresponding block structure based clustering. H-LU preconditioners resulting from the block and joint approaches will be discussed and compared through numerical results.

Full Text (PDF) [723 KB], BibTeX

Key words

hierarchical matrices, data-sparse approximation, Oseen equations, preconditioning, factorization

AMS subject classifications

65F05, 65F30, 65F50

Links to the cited ETNA articles

[20] Vol. 23 (2006), pp. 158-179 Olaf Schenk and Klaus Gärtner: On fast factorization pivoting methods for sparse symmetric indefinite systems