New preconditioners based on symmetric-triangular decomposition for saddle point problems

In this paper, we construct some new triangular preconditioners for saddle point problems based on the symmetric and triangular (ST) decomposition. Furthermore, we obtain some estimations on the condition number for the preconditioned systems and give the quasi-optimal parameters. Numerical experiments on the Stokes problems are given to illustrate fast convergence of the associated conjugate gradient method.

[1]  Martin Stoll,et al.  All-at-once preconditioning in PDE-constrained optimization , 2010, Kybernetika.

[2]  Jinyun Yuan,et al.  Numerical methods for generalized least squares problems , 1996 .

[3]  Xiaojun Chen,et al.  Validated Solutions of Saddle Point Linear Systems , 2008, SIAM J. Matrix Anal. Appl..

[4]  Xiaojun Chen,et al.  Numerical validation of solutions of saddle point matrix equations , 2003, Numer. Linear Algebra Appl..

[5]  Stephen J. Wright.  Stability of Augmented System Factorizations in Interior-Point Methods , 1997, SIAM J. Matrix Anal. Appl..

[6]  J. Pasciak,et al.  A preconditioning technique for indefinite systems resulting from mixed approximations of elliptic problems , 1988 .

[7]  Martin Stoll,et al.  Combination Preconditioning and the Bramble-Pasciak+ Preconditioner , 2008, SIAM J. Matrix Anal. Appl..

[8]  Zeng-Qi Wang,et al.  On parameterized inexact Uzawa methods for generalized saddle point problems , 2008 .

[9]  G. Golub,et al.  SYMMETRIC-TRIANGULAR DECOMPOSITION AND ITS APPLICATIONS , 2002 .

[10]  Gene H. Golub,et al.  Numerical solution of saddle point problems , 2005, Acta Numerica.

[11]  Nicholas I. M. Gould,et al.  Preconditioning Saddle-Point Systems with Applications in Optimization , 2010, SIAM J. Sci. Comput..

[12]  Gene H. Golub,et al.  Symmetric-Triangular Decomposition and its Applications Part I: Theorems and Algorithms , 2002 .

[13]  Beresford N. Parlett,et al.  On generalized successive overrelaxation methods for augmented linear systems , 2005, Numerische Mathematik.

[14]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[16]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[17]  H. Elman.  Multigrid and Krylov subspace methods for the discrete Stokes equations , 1994 .

[18]  Howard C. Elman,et al.  Algorithm 866: IFISS, a Matlab toolbox for modelling incompressible flow , 2007, TOMS.

[19]  Z. Bai,et al.  Restrictively preconditioned conjugate gradient methods for systems of linear equations , 2003 .

[20]  Jin Yun Yuan,et al.  Row-wise algorithms of ST decomposition , 2006, Appl. Math. Comput..