Stochastic Preconditioning for Diagonally Dominant Matrices

Abstract. This paper presents a new stochastic preconditioning approach for large sparse matrices. For the class of matrices that are row-wise and column-wise irreducibly diagonally dominant, we prove that an incomplete LDLT factorization in a symmetric case or an incomplete LDU factorization in an asymmetric case can be obtained from random walks, and used as a preconditioner. It is argued that our factor matrices have better quality, i.e., better accuracy-size tradeoffs, than preconditioners produced by existing incomplete factorization methods. Therefore a resulting preconditioned Krylov-subspace iterative solver requires less computation than traditional methods to solve a set of linear equations with the same error tolerance. The advantage increases for larger and denser matrices. These claims are verified by numerical tests, and we provide techniques that can potentially extend the theory to non-diagonally-dominant matrices. Key words. incomplete factorization, iterative solver, ...

• Type: Journal
• Year: 2008
• Where: SIAMSC
• Authors: Haifeng Qian, Sachin S. Sapatnekar
• Keywords: Factorization, Incomplete Factorization, Incomplete Ldlt Factorization, SIAMSC 2008