Hermitian and Skew-Hermitian Splitting Methods for Non-Hermitian Positive Definite Linear Systems

We study efficient iterative methods for the large sparse non-Hermitian positive definite system of linear equations based on the Hermitian and skew-Hermitian splitting of the coefficient matrix. These methods include a Hermitian/skew-Hermitian splitting (HSS) iteration and its inexact variant, the inexact Hermitian/skew-Hermitian splitting (IHSS) iteration, which employs some Krylov subspace methods as its inner iteration processes at each step of the outer HSS iteration. Theoretical analyses show that the HSS method converges unconditionally to the unique solution of the system of linear equations. Moreover, we derive an upper bound of the contraction factor of the HSS iteration which is dependent solely on the spectrum of the Hermitian part and is independent of the eigenvectors of the matrices involved. Numerical examples are presented to illustrate the effectiveness of both HSS and IHSS iterations. In addition, a model problem of a three-dimensional convection-diffusion equation is used to illustrate the advantages of our methods.

[1]  Gene H. Golub,et al.  Block Triangular and Skew-Hermitian Splitting Methods for Positive-Definite Linear Systems , 2005, SIAM J. Sci. Comput..

[2]  Gene H. Golub,et al.  Preconditioned HSS methods for the solution of non-Hermitian positive definite linear systems and applications to the discrete convection-diffusion equation , 2005, Numerische Mathematik.

[3]  Gene H. Golub,et al.  Preconditioned Hermitian and skew-Hermitian splitting methods for non-Hermitian positive semidefinite linear systems , 2004, Numerische Mathematik.

[4]  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.

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

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

[7]  Bai,et al.  A REGULARIZED CONJUGATE GRADIENT METHOD FOR SYMMETRIC POSITIVE DEFINITE SYSTEM OF LINEAR EQUATIONS , 2002 .

[8]  Z. Bai,et al.  Sufficient conditions for the convergent splittings of non-Hermitian positive definite matrices , 2001 .

[9]  Anne Greenbaum,et al.  Iterative methods for solving linear systems , 1997, Frontiers in applied mathematics.

[10]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[11]  R. S. Varga,et al.  Acceleration of Relaxation Methods for Non-Hermitian Linear Systems , 1992, SIAM J. Matrix Anal. Appl..

[12]  G. Golub.  Matrix computations , 1983 .

[13]  J. Strikwerda.  A Generalized Conjugate Gradient Method for Non-Symmetric Systems of Linear Equations. , 1981 .

[14]  O. Widlund.  A Lanczos Method for a Class of Nonsymmetric Systems of Linear Equations , 1978 .