A preconditioner for generalized saddle point problems

被引:374
作者
Benzi, M [1 ]
Golub, GH
机构
[1] Emory Univ, Dept Math & Comp Sci, Atlanta, GA 30322 USA
[2] Stanford Univ, Sci Comp & Computat Math Program, Stanford, CA 94305 USA
关键词
saddle point problems; matrix splittings; iterative methods; preconditioning;
D O I
10.1137/S0895479802417106
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
this paper we consider the solution of linear systems of saddle point type by preconditioned Krylov subspace methods. A preconditioning strategy based on the symmetric/skew-symmetric splitting of the coefficient matrix is proposed, and some useful properties of the preconditioned matrix are established. The potential of this approach is illustrated by numerical experiments with matrices from various application areas.
引用
收藏
页码:20 / 41
页数:22
相关论文
共 64 条
[1]  
[Anonymous], 1983, STUD MATH APPL
[2]   A null space algorithm for mixed finite-element approximations of Darcy's equation [J].
Arioli, M ;
Manzini, G .
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING, 2002, 18 (09) :645-657
[3]   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
[4]  
BAI ZZ, 2002, SCCM0212 STANF U DEP
[5]  
BANK RE, 1990, NUMER MATH, V56, P645, DOI 10.1007/BF01405194
[6]   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
[7]   Optimization of the hermitian and skew-Hermitian splitting iteration for saddle-point problems [J].
Benzi, M ;
Gander, MJ ;
Golub, GH .
BIT, 2003, 43 (05) :881-900
[8]  
Benzi M., 1993, REND MAT, V13, P275
[9]  
BERTACCINI D, 2002, SCCM0211 STANF U DEP
[10]  
Betts J.T., 2001, ADV DESIGN CONTROL