A Class of Symmetric Factored Approximate Inverses and Hybrid Two-Level Solver

A new class of symmetric factored approximate inverses is proposed and used in conjunction with the Preconditioned Conjugate Gradient method for solving sparse symmetric linear systems. Additionally, a new hybrid two-level solver is proposed utilizing a block independent set reordering, in order to create the two level hierarchy. The Schur complement is formed explicitly by inverting the blocks created by reordering. Then, the preconditioned conjugate gradient method is used in conjunction with the symmetric factored approximate inverse to solve the reduced order linear system. Furthermore, numerical results on the performance and convergence behavior for solving various model problems are presented.

[1]  Michele Benzi,et al.  Robust Approximate Inverse Preconditioning for the Conjugate Gradient Method , 2000, SIAM J. Sci. Comput..

[2]  Ruipeng Li,et al.  Matrix Reordering Using Multilevel Graph Coarsening for ILU Preconditioning , 2013, SIAM J. Sci. Comput..

[3]  Yousef Saad,et al.  ARMS: an algebraic recursive multilevel solver for general sparse linear systems , 2002, Numer. Linear Algebra Appl..

[4]  C. K. Filelis-Papadopoulos,et al.  Parallel multigrid algorithms based on generic approximate sparse inverses: an SMP approach , 2013, The Journal of Supercomputing.

[5]  George A. Gravvanis,et al.  High Performance Inverse Preconditioning , 2009 .

[6]  Miroslav Tuma,et al.  Factorized Approximate Inverses with Adaptive Dropping , 2016, SIAM J. Sci. Comput..

[7]  Carlo Janna,et al.  FSAIPACK: A Software Package for High-Performance Factored Sparse Approximate Inverse Preconditioning , 2015, ACM Trans. Math. Softw..

[8]  T. Manteuffel.  An incomplete factorization technique for positive definite linear systems , 1980 .

[9]  C. K. Filelis-Papadopoulos,et al.  A class of generic factored and multi-level recursive approximate inverse techniques for solving general sparse systems , 2016 .

[10]  Timothy A. Davis,et al.  The university of Florida sparse matrix collection , 2011, TOMS.

[11]  Carlo Janna,et al.  The Use of Supernodes in Factored Sparse Approximate Inverse Preconditioning , 2015, SIAM J. Sci. Comput..

[12]  L. Kolotilina,et al.  Factorized Sparse Approximate Inverse Preconditionings I. Theory , 1993, SIAM J. Matrix Anal. Appl..

[13]  Olaf Schenk,et al.  Solving unsymmetric sparse systems of linear equations with PARDISO , 2004, Future Gener. Comput. Syst..

[14]  M. Benzi,et al.  A comparative study of sparse approximate inverse preconditioners , 1999 .

[15]  Kengo Nakajima.  SPARSE APPROXIMATE INVERSE PRECONDITIONER FOR CONTACT PROBLEMS ON THE EARTH SIMULATOR USING OPENMP , 2008 .

[16]  Edmond Chow,et al.  A Priori Sparsity Patterns for Parallel Sparse Approximate Inverse Preconditioners , 1999, SIAM J. Sci. Comput..

[17]  Marcus J. Grote,et al.  Parallel Preconditioning with Sparse Approximate Inverses , 1997, SIAM J. Sci. Comput..

[18]  Michele Benzi,et al.  A Sparse Approximate Inverse Preconditioner for the Conjugate Gradient Method , 1996, SIAM J. Sci. Comput..

[19]  Jun Zhang,et al.  BILUTM: A Domain-Based Multilevel Block ILUT Preconditioner for General Sparse Matrices , 1999, SIAM J. Matrix Anal. Appl..

[20]  M. A. Ajiz,et al.  A robust incomplete Choleski‐conjugate gradient algorithm , 1984 .

[21]  Fred Wubs,et al.  Nested grids ILU-decomposition (NGILU) , 1996 .

[22]  Edmond Chow,et al.  Parallel Implementation and Practical Use of Sparse Approximate Inverse Preconditioners with a Priori Sparsity Patterns , 2001, Int. J. High Perform. Comput. Appl..

[23]  Masha Sosonkina,et al.  pARMS: a parallel version of the algebraic recursive multilevel solver , 2003, Numer. Linear Algebra Appl..

[24]  J. Meijerink,et al.  An iterative solution method for linear systems of which the coefficient matrix is a symmetric -matrix , 1977 .

[25]  D. J. Evans,et al.  Explicit semi-direct methods based on approximate inverse matrix techniques for solving boundary-value problems on parallel processors , 1987 .

[26]  Carlo Janna,et al.  Enhanced Block FSAI Preconditioning Using Domain Decomposition Techniques , 2013, SIAM J. Sci. Comput..