CONSTRAINT PRECONDITIONERS FOR SYMMETRIC INDEFINITE MATRICES

被引:146
作者
Bai, Zhong-Zhi [1 ]
Ng, Michael K. [2 ]
Wang, Zeng-Qi [3 ]
机构
[1] Chinese Acad Sci, Inst Computat Math & Sci Engn Comp, Acad Math & Syst Sci, State Key Lab Sci Engn Comp, Beijing 100080, Peoples R China
[2] Hong Kong Baptist Univ, Dept Math, Kowloon Tong, Hong Kong, Peoples R China
[3] Shanghai Jiao Tong Univ, Dept Math, Shanghai 200240, Peoples R China
关键词
symmetric indefinite systems; constraint preconditioners; FAST ITERATIVE SOLUTION; LARGE SPARSE EQUALITY; POINT; EIGENVALUES; ALGORITHM; SYSTEMS; PART;
D O I
10.1137/080720243
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We study the eigenvalue bounds of block two-by-two nonsingular and symmetric indeifnite matrices whose (1, 1) block is symmetric positive definite and Schur complement with respect to its (2, 2) block is symmetric indefinite. A constraint preconditioner for this matrix is constructed by simply replacing the (1, 1) block by a symmetric and positive definite approximation, and the spectral properties of the preconditioned matrix are discussed. Numerical results show that, for a suitably chosen (1, 1) block-matrix, this constraint preconditioner outperforms the block-diagonal and the block-tridiagonal ones in iteration step and computing time when they are used to accelerate the GMRES method for solving these block two-by-two symmetric positive indefinite linear systems. The new results extend the existing ones about block two-by-two matrices of symmetric negative semidefinite (2, 2) blocks to those of general symmetric (2, 2) blocks.
引用
收藏
页码:410 / 433
页数:24
相关论文
共 52 条
[1]   Multilevel preconditioned iterative eigensolvers for Maxwell eigenvalue problems [J].
Arbenz, P ;
Geus, R .
APPLIED NUMERICAL MATHEMATICS, 2005, 54 (02) :107-121
[2]   Eigenvalue estimates for preconditioned saddle point matrices [J].
Axelsson, O ;
Neytcheva, M .
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS, 2006, 13 (04) :339-360
[3]  
Babuka I., 1991, Handbook of Numerical Analysis, VII, P641
[4]  
Bai Z-Z., 2004, J. of Shanghai Univ, V8, P397, DOI [10.1007/s11741-004-0050-2, DOI 10.1007/S11741-004-0050-2]
[5]  
Bai ZZ, 2007, IMA J NUMER ANAL, V27, P1, DOI [10.1093/imanum/dr1017, 10.1093/imanum/drl017]
[6]  
Bai ZZ, 2006, MATH COMPUT, V75, P791, DOI 10.1090/S0025-5718-05-01801-6
[7]   On inexact preconditioners for nonsymmetric matrices [J].
Bai, ZZ ;
Ng, MK .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2005, 26 (05) :1710-1724
[8]   Hermitian and skew-Hermitian splitting methods for non-hermitian positive definite linear systems [J].
Bai, ZZ ;
Golub, GH ;
Ng, MK .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 2003, 24 (03) :603-626
[9]   A fast algorithm for the electromagnetic scattering from a large cavity [J].
Bao, G ;
Sun, WW .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2005, 27 (02) :553-574
[10]   CONVERGENCE OF THE FRANCIS SHIFTED QR ALGORITHM ON NORMAL MATRICES [J].
BATTERSON, S .
LINEAR ALGEBRA AND ITS APPLICATIONS, 1994, 207 :181-195