Using approximate inverses in algebraic multilevel methods (Q1272813)

This paper deals with the iterative solution of large sparse symmetric positive definite systems. The author investigates preconditioning techniques of the two-level type that are based on a block factorization of the system matrix. Whereas the basic scheme assumes an exact inversion of the submatrix related to the first block of unknowns, the author analyzes the effect of using an approximate inverse instead. Derived condition number estimates are valid for any type of approximation of the Schur complement and do not assume the use of the hierarchical basis. They show that the two-level methods are stable when using approximate inverses based on modified ILU techniques, or explicit inverses that meet some row-sum criterion. On the other hand, the author brings to the light that use of standard approximate inverses based on convergent splittings can have a dramatic effect on the convergence rate. These conclusions are numerically illustrated on some examples.