Preconditioners for nonsymmetric indefinite linear systems

Abstract In this paper, we develop algorithms for solving nonsymmetric indefinite linear systems by considering the augmented linear systems resulting from a weighted linear least squares problem. Even though the augmented system is more ill-conditioned than the original linear system, one can construct preconditioned GMRES methods for solving these augmented systems capable of obtaining reasonable approximation of the solution in fewer iterations than the classical ILU preconditioned GMRES method for solving the original linear system. More specifically, we present two different preconditioners for these augmented systems, examine the spectral properties of these preconditioned augmented systems, and report numerical results to illustrate the effectiveness of these preconditioners.

[1]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[2]  M. Benzi.  Preconditioning techniques for large linear systems: a survey , 2002 .

[3]  Gene H. Golub,et al.  A Preconditioner for Generalized Saddle Point Problems , 2004, SIAM J. Matrix Anal. Appl..

[4]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[5]  O. Axelsson.  Preconditioning of Indefinite Problems by Regularization , 1979 .

[6]  M. Stynes,et al.  Numerical methods for singularly perturbed differential equations : convection-diffusion and flow problems , 1996 .

[7]  Zhong-Zhi Bai.  Quasi-HSS iteration methods for non-Hermitian positive definite linear systems of strong skew-Hermitian parts , 2018, Numer. Linear Algebra Appl..

[8]  Louis A. Hageman,et al.  Iterative Solution of Large Linear Systems. , 1971 .

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

[10]  Jejoong Yoo,et al.  Modeling and simulation of ion channels. , 2012, Chemical reviews.

[11]  A. Wathen,et al.  Minimum residual methods for augmented systems , 1998 .

[12]  Y. Saad.  Preconditioning techniques for nonsymmetric and indefinite linear systems , 1988 .

[13]  Gene H. Golub,et al.  Hermitian and Skew-Hermitian Splitting Methods for Non-Hermitian Positive Definite Linear Systems , 2002, SIAM J. Matrix Anal. Appl..

[14]  T. F. Russell,et al.  NUMERICAL METHODS FOR CONVECTION-DOMINATED DIFFUSION PROBLEMS BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OR FINITE DIFFERENCE PROCEDURES* , 1982 .

[15]  Y. Saad,et al.  Experimental study of ILU preconditioners for indefinite matrices , 1997 .

[16]  Ahmed H. Sameh,et al.  A parallel sparse linear system solver based on Hermitian/skew-Hermitian splitting , 2016, Comput. Math. Appl..

[17]  Michele Benzi,et al.  Preconditioning Highly Indefinite and Nonsymmetric Matrices , 2000, SIAM J. Sci. Comput..

[18]  Owe Axelsson,et al.  A Class of Nested Iteration Schemes for Linear Systems with a Coefficient Matrix with a Dominant Positive Definite Symmetric Part , 2004, Numerical Algorithms.

[19]  Yin Zhang,et al.  Limited Memory Block Krylov Subspace Optimization for Computing Dominant Singular Value Decompositions , 2013, SIAM J. Sci. Comput..

[20]  Gene H. Golub,et al.  Numerical solution of saddle point problems , 2005, Acta Numerica.

[21]  Michele Benzi,et al.  Spectral Properties of the Hermitian and Skew-Hermitian Splitting Preconditioner for Saddle Point Problems , 2005, SIAM J. Matrix Anal. Appl..

[22]  Owe Axelsson,et al.  Preconditioning methods for linear systems arising in constrained optimization problems , 2003, Numer. Linear Algebra Appl..

[23]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[24]  G. Golub,et al.  Optimization of the Hermitian and Skew-Hermitian Splitting Iteration for Saddle-Point Problems , 2003 .