Semi-convergence analysis of the Uzawa-SOR methods for singular saddle point problems

Recently, Zhang and Shang proposed a class of Uzawa-SOR methods for solving the nonsingular saddle point problems; see Zhang and Shang (2010). In this paper, we give the semi-convergence analysis of this method when it is applied to solve the singular saddle point problems under some conditions. Finally, numerical experiments are presented to illustrate the feasibility and effectiveness of this method.     

Preconditioned AHSS iteration method for singular saddle point problems

In this paper, for solving the singular saddle point problems, we present a new preconditioned accelerated Hermitian and skew-Hermitian splitting (AHSS) iteration method. The semi-convergence of this method and the eigenvalue distribution of the preconditioned iteration matrix are studied. In addition, we prove that all eigenvalues of the iteration matrix are clustered for any positive iteration parameters α and β. Numerical experiments illustrate the theoretical results and examine the numerical effectiveness of the AHSS iteration method served either as a preconditioner or as a solver.     

A generalized preconditioned HSS method for singular saddle point problems

Li et al. recently studied the generalized HSS (GHSS) method for solving singular linear systems (see Li et al., J. Comput. Appl. Math. 236, 2338–2353 (2012)). In this paper, we generalize the method and present a generalized preconditioned Hermitian and skew-Hermitian splitting method (GPHSS) to solve singular saddle point problems. We prove the semi-convergence of GPHSS under some conditions, and weaken some semi-convergent conditions of GHSS, moreover, we analyze the spectral properties of the corresponding preconditioned matrix. Numerical experiments are given to illustrate the efficiency of GPHSS method with appropriate parameters both as a solver and as a preconditioner.     

A generalized variant of the deteriorated PSS preconditioner for nonsymmetric saddle point problems

Based on the variant of the deteriorated positive-definite and skew-Hermitian splitting (VDPSS) preconditioner developed by Zhang and Gu (BIT Numer. Math. 56:587–604, 2016), a generalized VDPSS (GVDPSS) preconditioner is established in this paper by replacing the parameter α in (2,2)-block of the VDPSS preconditioner by another parameter β. This preconditioner can also be viewed as a generalized form of the VDPSS preconditioner and the new relaxed HSS (NRHSS) preconditioner which has been exhibited by Salkuyeh and Masoudi (Numer. Algorithms, 2016). The convergence properties of the GVDPSS iteration method are derived. Meanwhile, the distribution of eigenvalues and the forms of the eigenvectors of the preconditioned matrix are analyzed in detail. We also study the upper bounds on the degree of the minimum polynomial of the preconditioned matrix. Numerical experiments are implemented to illustrate the effectiveness of the GVDPSS preconditioner and verify that the GVDPSS preconditioned generalized minimal residual method is superior to the DPSS, relaxed DPSS, SIMPLE-like, NRHSS, and VDPSS preconditioned ones for solving saddle point problems in terms of the iterations and computational times.     

An efficient numerical method for preconditioned saddle point problems

In this paper, we consider the solution of linear systems of saddle point type by a preconditioned numerical method. We first transform the original linear system into two sub-systems with small size by a preconditioning strategy, then employ the conjugate gradient (CG) method to solve the linear system with a SPD coefficient matrix, and a splitting iteration method to solve the other sub-system, respectively. Numerical experiments show that the new method can achieve faster convergence than several effective preconditioners published in the recent literature in terms of total runtime and iteration steps.     

On the UPSS method for non-Hermitian singular saddle point problems

Recently, a new Uzawa-type method, referred as the UPSS method, is proposed for solving the non-Hermitian nonsingular saddle point problems, see Dou, Yang and Wu (2017). In this paper, we give the semi-convergence analysis of the UPSS method when it is used to solve non-Hermitian singular saddle point problems. An example is given to verify the effectiveness of this method for solving non-Hermitian singular saddle point problems.     

A generalized modified SOR-like method for the singular saddle point problems

Recently, Guo et al. proposed a modified SOR-like (MSOR-like) iteration method for solving the nonsingular saddle point problem. In this paper, we further prove the semi-convergence of this method when it is applied to solve the singular saddle point problems under suitable conditions on the involved iteration parameters. Moreover, the optimal iteration parameters and the corresponding optimal semi-convergence factor for the MSOR-like method are determined. In addition, numerical experiments are used to show the feasibility and effectiveness of the MSOR-like method for solving singular saddle point problems, arising from the incompressible flow problems.     

鞍点问题的修正对称超松弛迭代算法

为了提高求解鞍点问题的迭代算法的速度,通过设置合适的加速变量,对修正超松弛迭代算法（简记作MSOR-like算法）和广义对称超松弛迭代算法（简记作GSSOR-like算法）进行了修正,给出了修正对称超松弛迭代算法,即MSSOR-like （modified symmetric successiveover-relaxation）算法,并研究了该算法收敛的充分必要条件.最后,通过数值例子表明,选择合适的参数后,新算法的迭代速度和迭代次数均优于MSOR-like （modified successive overrelaxation）和GSSOR-like （generalized symmetric successive over-relaxation）算法,因此,它是一种较好的解决鞍点问题的算法.     

The spectral properties of the preconditioned matrix for nonsymmetric saddle point problems

In this paper, on the basis of matrix splitting, two preconditioners are proposed and analyzed, for nonsymmetric saddle point problems. The spectral property of the preconditioned matrix is studied in detail. When the iteration parameter becomes small enough, the eigenvalues of the preconditioned matrices will gather into two clusters—one is near (0,0) and the other is near (2,0)—for the PPSS preconditioner no matter whether A is Hermitian or non-Hermitian and for the PHSS preconditioner when A is a Hermitian or real normal matrix. Numerical experiments are given, to illustrate the performances of the two preconditioners.     

Eigenvalue estimates for preconditioned saddle point matrices

New accurate eigenvalue bounds for symmetric matrices of saddle point form are derived and applied for both unpreconditioned and preconditioned versions of the matrices. The estimates enable a better understanding of how preconditioners should be chosen. The preconditioners provide efficient iterative solution of the corresponding linear systems with, for some important applications, an optimal order of computational complexity. The methods are applied for Stokes problem and for linear elasticity problems. Copyright © 2005 John Wiley & Sons, Ltd.     

关于广义鞍点问题的HSS迭代方法的收缩因子(英文)

白中治等提出了解非埃尔米特正定线性方程组的埃尔米特和反埃尔米特分裂(HSS)迭代方法(Bai Z Z,Golub G H,Ng M K.Hermitian and skew-Hermitian splitting methodsfor non-Hermitian positive definite linear systems.SIAM J.Matrix Anal.Appl.,2003,24:603-626).本文精确地估计了用HSS迭代方法求解广义鞍点问题时在加权2-范数和2-范数下的收缩因子.在实际的计算中,正是这些收缩因子而不是迭代矩阵的谱半径,本质上控制着HSS迭代方法的实际收敛速度.根据文中的分析,求解广义鞍点问题的HSS迭代方法的收缩因子在加权2-范数下等于1,在2-范数下它会大于等于1,而在某种适当选取的范数之下,它则会小于1.最后,用数值算例说明了理论结果的正确性.     

On lower iteration complexity bounds for the convex concave saddle point problems

Mathematical Programming - In this paper, we study the lower iteration complexity bounds for finding the saddle point of a strongly convex and strongly concave saddle point problem: $$min _xmax...     

On an iterative method for saddle point problems

Convergence behavior of a nested iterative scheme presented in a paper by Bank, Welfert and Yserentant is studied. It is shown that this scheme converges under conditions weaker than that stated in their paper. Received November 18, 1996 / Revised version received March 18, 1997     

An iterative perturbation method for saddle point problems

An iterative perturbation method is presented for saddle pointproblems a as an iterative scheme for mixed finite element methods.It can be viewed as an extension of augmented Lagrangian methodsto second-order elliptic equations with a reaction term. Thisiterative method can lead to positive definite systems and decouplethe flux variable from the pressure variable. Convergence analysesare established at the differential and finite-dimensional levels.Various computational experiments for problems with continuousand discontinuous coefficients and full diffusion tensors areconducted to confirm the theoretical claims. These experimentsare implemented in the object oriented programming paradigmusing C++ to take advantage of object oriented software engineeringtechniques such as encapsulation, inheritance, and polymorphism.     

The alternating-direction iterative method for saddle point problems

In the paper, a new alternating-direction iterative method is proposed based on matrix splittings for solving saddle point problems. The convergence analysis for the new method is given. When the better values of parameters are employed, the proposed method has faster convergence rate and less time cost than the Uzawa algorithm with the optimal parameter and the Hermitian and skew-Hermitian splitting iterative method. Numerical examples further show the effectiveness of the method.     

On local Hermitian and skew-Hermitian splitting iteration methods for generalized saddle point problems

In this paper, we first present a local Hermitian and skew-Hermitian splitting (LHSS) iteration method for solving a class of generalized saddle point problems. The new method converges to the solution under suitable restrictions on the preconditioning matrix. Then we give a modified LHSS (MLHSS) iteration method, and further extend it to the generalized saddle point problems, obtaining the so-called generalized MLHSS (GMLHSS) iteration method. Numerical experiments for a model Navier-Stokes problem are given, and the results show that the new methods outperform the classical Uzawa method and the inexact parameterized Uzawa method.     

Modified PHSS iterative methods for solving nonsingular and singular saddle point problems

Numerical Algorithms - For large sparse saddle point problems, we establish a new version of the preconditioned Hermitian and skew-Hermitian splitting (PHSS) iteration method, called the modified...