The semi-convergence properties of MHSS method for a class of complex nonsymmetric singular linear systems

We use the modified Hermitian and skew-Hermitian splitting (MHSS) iteration method to solve a class of complex nonsymmetric singular linear systems. The semi-convergence properties of the MHSS method are studied by analyzing the spectrum of the iteration matrix. Moreover, after investigating the semi-convergence factor and estimating its upper bound for the MHSS iteration method, an optimal iteration parameter that minimizes the upper bound of the semi-convergence factor is obtained. Numerical experiments are used to illustrate the theoretical results and examine the effectiveness of the MHSS method served both as a preconditioner for GMRES method and as a solver.     

On semi-convergence of modified HSS method for a class of complex singular linear systems

In this paper, we discuss the semi-convergence of the modified Hermitian and skew-Hermitian splitting (MHSS) iteration method for solving a broad class of complex singular linear systems. Some semi-convergence theories of the MHSS iteration method are established and are weaker than those appeared in previously published works.     

Complex-extrapolated MHSS iteration method for singular complex symmetric linear systems

In this paper, by extrapolating the modified Hermitian and skew-Hermitian splitting (MHSS) iteration method with a complex relaxation parameter, a complex-extrapolated MHSS (CMHSS) iteration method is present for solving a class of complex singular symmetric of linear equations. We study the semi-convergence properties of the CMHSS iteration method and the extent of the optimal iterative parameters. Furthermore, the convergence conditions also hold for solving nonsingular complex systems. Numerical experiments are given to verify the effectiveness of the CMHSS iteration method for solving both singular and nonsingular complex symmetric systems.     

一类弱非线性方程组的Picard-MHSS迭代方法

修正的Hermite/反Hermite分裂(MHSS)迭代方法是一类求解大型稀疏复对称线性代数方程组的无条件收敛的迭代算法.基于非线性代数方程组的特殊结构和性质,我们选取Picard迭代为外迭代方法,MHSS迭代作为内迭代方法,构造了求解大型稀疏弱非线性代数方程组的Picard-MHSS和非线性MHSS-like方法.这两类方法的优点是不需要在每次迭代时均精确计算和存储Jacobi矩阵,仅需要在迭代过程中求解两个常系数实对称正定子线性方程组.除此之外,在一定条件下,给出了两类方法的局部收敛性定理.数值结果证明了这两类方法是可行、有效和稳健的.     

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.     

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.     

Semi-convergence analysis of preconditioned deteriorated PSS iteration method for singular saddle point problems

In this paper, we propose a two-parameter preconditioned variant of the deteriorated PSS iteration method (J. Comput. Appl. Math., 273, 41–60 (2015)) for solving singular saddle point problems. Semi-convergence analysis shows that the new iteration method is convergent unconditionally. The new iteration method can also be regarded as a preconditioner to accelerate the convergence of Krylov subspace methods. Eigenvalue distribution of the corresponding preconditioned matrix is presented, which is instructive for the Krylov subspace acceleration. Note that, when the leading block of the saddle point matrix is symmetric, the new iteration method will reduce to the preconditioned accelerated HSS iteration method (Numer. Algor., 63 (3), 521–535 2013), the semi-convergence conditions of which can be simplified by the results in this paper. To further improve the effectiveness of the new iteration method, a relaxed variant is given, which has much better convergence and spectral properties. Numerical experiments are presented to investigate the performance of the new iteration methods for solving singular saddle point problems.     

Accelerated PMHSS iteration methods for complex symmetric linear systems

In this paper, we introduce and analyze an accelerated preconditioning modification of the Hermitian and skew-Hermitian splitting (APMHSS) iteration method for solving a broad class of complex symmetric linear systems. This accelerated PMHSS algorithm involves two iteration parameters α,β and two preconditioned matrices whose special choices can recover the known PMHSS (preconditioned modification of the Hermitian and skew-Hermitian splitting) iteration method which includes the MHSS method, as well as yield new ones. The convergence theory of this class of APMHSS iteration methods is established under suitable conditions. Each iteration of this method requires the solution of two linear systems with real symmetric positive definite coefficient matrices. Theoretical analyses show that the upper bound σ₁(α,β) of the asymptotic convergence rate of the APMHSS method is smaller than that of the PMHSS iteration method. This implies that the APMHSS method may converge faster than the PMHSS method. Numerical experiments on a few model problems are presented to illustrate the theoretical results and examine the numerical effectiveness of the new method.     

The generalized HSS method for solving singular linear systems

For the singular, non-Hermitian, and positive semidefinite linear systems, we propose an alternating-direction iterative method with two parameters based on the Hermitian and skew-Hermitian splitting. The semi-convergence analysis and the quasi-optimal parameters of the proposed method are discussed. Moreover, the corresponding preconditioner based on the splitting is given to improve the semi-convergence rate of the GMRES method. Numerical examples are given to illustrate the theoretical results and the efficiency of the generalized HSS method either as a solver or a preconditioner for GMRES.     

An Iterative Multigrid Regularization Method for Toeplitz Discrete Ill-Posed Problems

Iterative regularization multigrid methods have been successfully applied to signal/image deblurring problems. When zero-Dirichlet boundary conditions are imposed the deblurring matrix has a Toeplitz structure and it is potentially full. A crucial task of a multilevel strategy is to preserve the Toeplitz structure at the coarse levels which can be exploited to obtain fast computations. The smoother has to be an iterative regularization method. The grid transfer operator should preserve the regularization property of the smoother. This paper improves the iterative multigrid method proposed in [11] introducing a wavelet soft-thresholding denoising post-smoother. Such post-smoother avoids the noise amplification that is the cause of the semi-convergence of iterative regularization methods and reduces ringing effects. The resulting iterative multigrid regularization method stabilizes the iterations so that the imprecise (over) estimate of the stopping iteration does not have a deleterious effect on the computed solution. Numerical examples of signal and image deblurring problems confirm the effectiveness of the proposed method.     

AIR Tools — A MATLAB package of algebraic iterative reconstruction methods

We present a MATLAB package with implementations of several algebraic iterative reconstruction methods for discretizations of inverse problems. These so-called row action methods rely on semi-convergence for achieving the necessary regularization of the problem. Two classes of methods are implemented: Algebraic Reconstruction Techniques (ART) and Simultaneous Iterative Reconstruction Techniques (SIRT). In addition we provide a few simplified test problems from medical and seismic tomography. For each iterative method, a number of strategies are available for choosing the relaxation parameter and the stopping rule. The relaxation parameter can be fixed, or chosen adaptively in each iteration; in the former case we provide a new “training” algorithm that finds the optimal parameter for a given test problem. The stopping rules provided are the discrepancy principle, the monotone error rule, and the NCP criterion; for the first two methods “training” can be used to find the optimal discrepancy parameter.     

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.     

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.     

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.     

Minimum residual modified HSS iteration method for a class of complex symmetric linear systems

Numerical Algorithms - For solving a class of complex symmetric system of linear equations, we apply the minimum residual technique to the modified Hermitian and skew-Hermitian splitting (MHSS)...     

极大极小随机规划逼近问题最优解集和最优值的稳定性

研究了特殊的二层极大极小随机规划逼近收敛问题. 首先将下层初始随机规划最优解集拓展到非单点集情形, 且可行集正则的条件下, 讨论了下层随机规划逼近问题最优解集关于上层决策变量参数的上半收敛性和最优值函数的连续性. 然后把下层随机规划的epsilon-最优解向量函数反馈到上层随机规划的目标函数中, 得到了上层随机规划逼近问题的最优解集关于最小信息概率度量收敛的上半收敛性和最优值的连续性.     

On a Chebyshev accelerated splitting iteration method with application to two‐by‐two block linear systems

The Chebyshev accelerated preconditioned modified Hermitian and skew‐Hermitian splitting (CAPMHSS) iteration method is presented for solving the linear systems of equations, which have two‐by‐two block coefficient matrices. We derive an iteration error bound to show that the new method is convergent as long as the eigenvalue bounds are not underestimated. Even when the spectral information is lacking, the CAPMHSS iteration method could be considered as an exponentially converging iterative scheme for certain choices of the method parameters. In this case, the convergence rate is independent of the parameters. Besides, the linear subsystems in each iteration can be solved inexactly, which leads to the inexact CAPMHSS iteration method. The iteration error bound of the inexact method is derived also. We discuss in detail the implementation of CAPMHSS for solving two models arising from the Galerkin finite‐element discretizations of distributed control problems and complex symmetric linear systems. The numerical results show the robustness and the efficiency of the new methods.     

函数方程求根的一种新型大范围收敛迭代法

1引言实际解函数方程f(x)=0(z∈[α,β],f(α)f(β)＜0,f(x~*)=0)时,人们常常希望选用那些仅计算函数值,具有大范围收敛性且效率较高的方法,特别对那些表示式复杂的函数以及病态函数．例如,那种仅在的某个充分小邻域内连续,而在该邻域之外光滑性很差的函数；那种在初始含根区间(α,β)上起伏多变的函数；那种|f(α)|和|f(β)|差别甚大,而x~*又十分靠近绝对值较大者一端的函数等等,这种欲望就更加强烈．     

A modified generalized shift-splitting preconditioner for nonsymmetric saddle point problems

For the nonsymmetric saddle point problems with nonsymmetric positive definite (1,1) parts, the modified generalized shift-splitting (MGSS) preconditioner as well as the MGSS iteration method is derived in this paper, which generalize the modified shift-splitting (MSS) preconditioner and the MSS iteration method newly developed by Huang and Su (J. Comput. Appl. Math. 317:535–546, 2017), respectively. The convergent and semi-convergent analyses of the MGSS iteration method are presented, and we prove that this method is unconditionally convergent and semi-convergent. Meanwhile, some spectral properties of the preconditioned matrix are carefully analyzed. Numerical results demonstrate the robustness and effectiveness of the MGSS preconditioner and the MGSS iteration method and also illustrate that the MGSS iteration method outperforms the generalized shift-splitting (GSS) and the generalized modified shift-splitting (GMSS) iteration methods, and the MGSS preconditioner is superior to the shift-splitting (SS), GSS, modified SS (M-SS), GMSS and MSS preconditioners for the generalized minimal residual (GMRES) method for solving the nonsymmetric saddle point problems.     

The approximate and exact solutions of the space- and time-fractional Burgers equations with initial conditions by variational iteration method

In this paper, a scheme is developed to study numerical solution of the space- and time-fractional Burgers equations with initial conditions by the variational iteration method (VIM). The exact and numerical solutions obtained by the variational iteration method are compared with that obtained by Adomian decomposition method (ADM). The results show that the variational iteration method is much easier, more convenient, and more stable and efficient than Adomian decomposition method. Numerical solutions are calculated for the fractional Burgers equation to show the nature of solution as the fractional derivative parameter is changed.