CONSTRAINT PRECONDITIONERS FOR SYMMETRIC INDEFINITE MATRICES

被引:147
作者
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 条
[21]   Implicit-factorization preconditioning and iterative solvers for regularized saddle-point systems [J].
Dollar, HS ;
Gould, NIM ;
Schilders, WHA ;
Wathen, AJ .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 2006, 28 (01) :170-189
[22]   Approximate factorization constraint preconditioners for saddle-point matrices [J].
Dollar, HS ;
Wathen, AJ .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2006, 27 (05) :1555-1572
[23]   Eigenmodes of isospectral drums [J].
Driscoll, TA .
SIAM REVIEW, 1997, 39 (01) :1-17
[24]   Indefinitely preconditioned conjugate gradient method for large sparse equality and inequality constrained quadratic problems [J].
Durazzi, C ;
Ruggiero, V .
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS, 2003, 10 (08) :673-688
[25]  
ELMAN H, 2005, FINITE ELEMENTS FAST
[26]  
Golub G. H., 2012, Matrix computations, V4th
[27]   On the solution of equality constrained quadratic programming problems arising in optimization [J].
Gould, NIM ;
Hribar, ME ;
Nocedal, J .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2001, 23 (04) :1375-1394
[28]   On the computation of a very large number of eigenvalues for selfadjoint elliptic operators by means of multigrid methods [J].
Heuveline, V .
JOURNAL OF COMPUTATIONAL PHYSICS, 2003, 184 (01) :321-337
[29]  
Heuveline V., 2000, E W J NUMER MATH, V8, P275
[30]   A note on preconditioning nonsymmetric matrices [J].
Ipsen, ICF .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2001, 23 (03) :1050-1051