scispace - formally typeset
Search or ask a question

Showing papers by "Zhong-Zhi Bai published in 2007"


01 Jan 2007
TL;DR: An upper bound of the contraction factor of the HSS iteration is derived which is dependent solely on the spectrum of the Hermitian part and is independent of the eigenvectors of the matrices involved.
Abstract: 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.

760 citations


Journal ArticleDOI
TL;DR: These methods involve two iteration parameters whose special choices can recover the known preconditioned HSS iteration methods, as well as yield new ones, and show that the new methods converge unconditionally to the unique solution of the saddle-point problem.
Abstract: We establish a class of accelerated Hermitian and skew-Hermitian splitting (AHSS) iteration methods for large sparse saddle-point problems by making use of the Hermitian and skew-Hermitian splitting (HSS) iteration technique. These methods involve two iteration parameters whose special choices can recover the known preconditioned HSS iteration methods, as well as yield new ones. Theoretical analyses show that the new methods converge unconditionally to the unique solution of the saddle-point problem. Moreover, the optimal choices of the iteration parameters involved and the corresponding asymptotic convergence rates of the new methods are computed exactly. In addition, theoretical properties of the preconditioned Krylov subspace methods such as GMRES are investigated in detail when the AHSS iterations are employed as their preconditioners. Numerical experiments confirm the correctness of the theory and the effectiveness of the methods.

345 citations


Journal ArticleDOI
TL;DR: This work derives necessary and sufficient conditions for guaranteeing the unconditional convergence of the preconditioned Hermitian and skew-Hermitian splitting iteration methods and applies these results to block tridiagonal linear systems in order to obtain convergence conditions for the corresponding block variants of these methods.
Abstract: For the non-Hermitian and positive semidefinite systems of linear equations, we derive necessary and sufficient conditions for guaranteeing the unconditional convergence of the preconditioned Hermitian and skew-Hermitian splitting iteration methods. We then apply these results to block tridiagonal linear systems in order to obtain convergence conditions for the corresponding block variants of the preconditioned Hermitian and skew-Hermitian splitting iteration methods.

178 citations


Journal ArticleDOI
TL;DR: Theoretical analyses show that the NSS method converges unconditionally to the exact solution of the system of linear equations, and an upper bound of the contraction factor is derived which is dependent solely on the spectrum of the normal splitting matrix, and is independent of the eigenvectors of the matrices involved.
Abstract: We further generalize the technique for constructing the Hermitian/skew-Hermitian splitting (HSS) iteration method for solving large sparse non-Hermitian positive definite system of linear equations to the normal/skew-Hermitian (NS) splitting obtaining a class of normal/skew-Hermitian splitting (NSS) iteration methods. Theoretical analyses show that the NSS method converges unconditionally to the exact solution of the system of linear equations. Moreover, we derive an upper bound of the contraction factor of the NSS iteration which is dependent solely on the spectrum of the normal splitting matrix, and is independent of the eigenvectors of the matrices involved. We present a successive-overrelaxation (SOR) acceleration scheme for the NSS iteration, which specifically results in an acceleration scheme for the HSS iteration. Convergence conditions for this SOR scheme are derived under the assumption that the eigenvalues of the corresponding block Jacobi iteration matrix lie in certain regions in the complex plane. A numerical example is used to show that the SOR technique can significantly accelerate the convergence rate of the NSS or the HSS iteration method. Copyright © 2006 John Wiley & Sons, Ltd.

173 citations


Journal ArticleDOI
TL;DR: Global convergence properties of the proposed methods are established and numerical results are provided, showing that the new method, called the Newton-GMRES with quasi-conjugate-gradient backtracking (NGQCGB), is very robust and effective.

82 citations


Journal ArticleDOI
TL;DR: In this article, the authors review the recently developed Hermitian/skew-Hermitian splitting (HSS) iteration, normal/Skew Hermitians splitting (NSS), PSS iteration, and block triangular/Hermitians split (BTSS) iteration.
Abstract: For large sparse system of linear equations with a non-Hermitian positive definite coefficient matrix, we review the recently developed Hermitian/skew-Hermitian splitting (HSS) iteration, normal/skew-Hermitian splitting (NSS) iteration, positive-definite/skew-Hermitian splitting (PSS) iteration, and block triangular/skew-Hermitian splitting (BTSS) iteration. These methods converge unconditionally to the exact solution of the linear system, with the upper bounds of their convergence factors being only dependent on the spectrum of the Hermitian (or normal, or positive-definite) splitting matrix, but independent of the spectrum of the skew-Hermitian splitting matrix as well as the eigenvectors of all matrices involved.

49 citations


Journal ArticleDOI
TL;DR: It is found that a combination of the Newton/fixed-point iteration with the preconditioned GMRES method is quite efficient for the Sinc-Galerkin discretization of the Burgers equations.
Abstract: We study the Newton method and a fixed-point method for solving the system of nonlinear equations arising from the Sinc-Galerkin discretization of the Burgers equations. In each step of the Newton method or the fixed-point method, a structured subsystem of linear equations is obtained and needs to be solved numerically. In this paper, preconditioning techniques are applied to solve such linear subsystems. The bounds for eigenvalues of the preconditioned matrices are derived and numerical examples are given to illustrate the effectiveness of the proposed methods. We also find that a combination of the Newton/fixed-point iteration with the preconditioned GMRES method is quite efficient for the Sinc-Galerkin discretization of the Burgers equations.

43 citations