Hermitian and skew-Hermitian splitting methods for non-hermitian positive definite linear systems

被引:956
作者
Bai, ZZ
Golub, GH
Ng, MK
机构
[1] Chinese Acad Sci, State Key Lab Sci Engn Comp, Inst Computat Math & Sci Engn Comp, Acad Math & Syst Sci, Beijing 100080, Peoples R China
[2] Stanford Univ, Dept Comp Sci, Stanford, CA 94305 USA
[3] Univ Hong Kong, Dept Math, Hong Kong, Hong Kong, Peoples R China
关键词
non-Hermitian matrix; splitting; Hermitian matrix; skew-Hermitian matrix; iterative methods;
D O I
10.1137/S0895479801395458
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
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.
引用
收藏
页码:603 / 626
页数:24
相关论文
共 27 条
[1]  
AXELSSON O, IN PRESS NUMER ALGOR
[2]   Sharp error bounds of some Krylov subspace methods for non-Hermitian linear systems [J].
Bai, ZZ .
APPLIED MATHEMATICS AND COMPUTATION, 2000, 109 (2-3) :273-285
[3]   A class of incomplete orthogonal factorization methods. I: Methods and theories [J].
Bai, ZZ ;
Duff, IS ;
Wathen, AJ .
BIT, 2001, 41 (01) :53-70
[4]   Existence and uniqueness of splittings for stationary iterative methods with applications to alternating methods [J].
Benzi, M ;
Szyld, DB .
NUMERISCHE MATHEMATIK, 1997, 76 (03) :309-321
[5]   Conjugate gradient methods for toeplitz systems [J].
Chan, RH ;
Ng, MK .
SIAM REVIEW, 1996, 38 (03) :427-482
[6]   Block-circulant preconditioners for systems arising from discretization of the three-dimensional convection-diffusion equation [J].
Cheung, WM ;
Ng, MK .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2002, 140 (1-2) :143-158
[7]  
CONCUS P, 1976, STANCS76535 STANF U
[8]  
Douglas Jr J., 1962, Numer. Math., V4, P41, DOI DOI 10.1007/BF01386295
[9]   ACCELERATION OF RELAXATION METHODS FOR NON-HERMITIAN LINEAR-SYSTEMS [J].
EIERMANN, M ;
NIETHAMMER, W ;
VARGA, RS .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 1992, 13 (03) :979-991
[10]   ITERATIVE METHODS FOR CYCLICALLY REDUCED NON-SELF-ADJOINT LINEAR-SYSTEMS [J].
ELMAN, HC ;
GOLUB, GH .
MATHEMATICS OF COMPUTATION, 1990, 54 (190) :671-700