Incomplete LU factorization

In numerical linear algebra, an incomplete LU factorization (abbreviated as ILU) of a matrix is a sparse approximation of the LU factorization often used as a preconditioner.

Introduction

edit

Consider a sparse linear system  . These are often solved by computing the factorization  , with L lower unitriangular and U upper triangular. One then solves  ,  , which can be done efficiently because the matrices are triangular.

For a typical sparse matrix, the LU factors can be much less sparse than the original matrix — a phenomenon called fill-in. The memory requirements for using a direct solver can then become a bottleneck in solving linear systems. One can combat this problem by using fill-reducing reorderings of the matrix's unknowns, such as the Minimum degree algorithm.

An incomplete factorization instead seeks triangular matrices L, U such that   rather than  . Solving for   can be done quickly but does not yield the exact solution to  . So, we instead use the matrix   as a preconditioner in another iterative solution algorithm such as the conjugate gradient method or GMRES.

Definition

edit

For a given matrix   one defines the graph   as

 

which is used to define the conditions a sparsity patterns   needs to fulfill

 

A decomposition of the form   where the following hold

  •   is a lower unitriangular matrix
  •   is an upper triangular matrix
  •   are zero outside of the sparsity pattern:  
  •   is zero within the sparsity pattern:  

is called an incomplete LU decomposition (with respect to the sparsity pattern  ).

The sparsity pattern of L and U is often chosen to be the same as the sparsity pattern of the original matrix A. If the underlying matrix structure can be referenced by pointers instead of copied, the only extra memory required is for the entries of L and U. This preconditioner is called ILU(0).

Stability

edit

Concerning the stability of the ILU the following theorem was proven by Meijerink and van der Vorst.[1]

Let   be an M-matrix, the (complete) LU decomposition given by  , and the ILU by  . Then

 

holds. Thus, the ILU is at least as stable as the (complete) LU decomposition.

Generalizations

edit

One can obtain a more accurate preconditioner by allowing some level of extra fill in the factorization. A common choice is to use the sparsity pattern of A2 instead of A; this matrix is appreciably more dense than A, but still sparse over all. This preconditioner is called ILU(1). One can then generalize this procedure; the ILU(k) preconditioner of a matrix A is the incomplete LU factorization with the sparsity pattern of the matrix Ak+1.

More accurate ILU preconditioners require more memory, to such an extent that eventually the running time of the algorithm increases even though the total number of iterations decreases. Consequently, there is a cost/accuracy trade-off that users must evaluate, typically on a case-by-case basis depending on the family of linear systems to be solved.

An approximation to the ILU factorization can be performed as a fixed-point iteration in a highly parallel way.[2]

See also

edit

References

edit
  • Saad, Yousef (1996), Iterative methods for sparse linear systems (1st ed.), Boston: PWS, ISBN 978-0-534-94776-7. See Section 10.3 and further.
  1. ^ Meijerink, J. A.; Vorst, Van Der; A, H. (1977). "An iterative solution method for linear systems of which the coefficient matrix is a symmetric 𝑀-matrix". Mathematics of Computation. 31 (137): 148–162. doi:10.1090/S0025-5718-1977-0438681-4. ISSN 0025-5718.
  2. ^ Chow, Edmond; Patel, Aftab (2015). "Fine-grained parallel incomplete LU factorization". SIAM Journal on Scientific Computing. 37 (2): C169-C193. doi:10.1137/140968896.


edit