IMA Journal of Numerical Analysis Advance Access published online on March 14, 2008
IMA Journal of Numerical Analysis, doi:10.1093/imanum/drn005
Numerically stable LDLT-factorization of 
-type saddle point matrices
Research Institute of Mathematics and Computing Science, University of Groningen, Nijenborgh 9, 9747 AG Groningen, The Netherlands, and Witteveen+Bos, Consulting Engineers, van Twickelostraat 2, 7411 SC Deventer, The Netherlands
Research Institute of Mathematics and Computing Science, University of Groningen, Nijenborgh 9, 9747 AG Groningen, The Netherlands
Email: wubs{at}math.rug.nl
Received on 7 December 2006. Accepted for publication 4 December 2007.
 |
Abstract |
|---|
We present a new algorithm that constructs a fill-reducing ordering for a special class of saddle point matrices: the -matrices. This class contains the matrix occurring after discretization of the Stokes equation on a C-grid. The commonly used approach is to construct a fill-reducing ordering for the whole matrix followed by an adaptation of the ordering such that it becomes feasible. We propose to compute first a fill-reducing ordering for an extension of the definite submatrix. This ordering can be easily extended to an ordering for the whole matrix. In this manner, the construction of the ordering is straightforward and it can be computed efficiently. We show that much of the structure of the matrix is preserved during Gaussian elimination. For an -matrix, the preserved structure allows us to prove that any feasible ordering obtained in this way is numerically stable. The growth factor of this factorization is much smaller than the one for general indefinite matrices and is bounded by a number that depends linearly on the number of indefinite nodes. The algorithm allows for generalization to saddle point problems that are not of -type and are nonsymmetric, e.g. the incompressible Navier–Stokes equations (with Coriolis force) on a C-grid. Numerical results for -matrices show that the algorithm is able to produce a factorization with low fill.
Key Words: saddle point problem; indefinite matrix; f-matrix; factorization; numerical stability; growth factor; C-grid; (Navier–)Stokes equations; electrical networks