Solving a system with an incomplete factorization preconditioner



next up previous contents index
Next: Point incomplete factorizations Up: Creating an incomplete Previous: Creating an incomplete

Solving a system with an incomplete factorization preconditioner

 

Incomplete factorizations can be given in various forms. If (with and nonsingular triangular matrices), solving a system proceeds in the usual way (figure gif),

  
Figure: Preconditioner solve of a system , with

but often incomplete factorizations are given as (with diagonal, and and now strictly triangular matrices, determined through the factorization process). In that case, one could use either of the following equivalent formulations for :

or

In either case, the diagonal elements are used twice (not three times as the formula for would lead one to expect), and since only divisions with are performed, storing explicitly is the practical thing to do.

  
Figure: Preconditioner solve of a system , with .

At the cost of some extra storage, one could store or , thereby saving some computation. Solving a system using the first formulation is outlined in figure gif. The second formulation is slightly harder to implement.



Jack Dongarra
Mon Nov 20 08:52:54 EST 1995