广义鞍点问题的松弛维数分解预条件子

本文将Benzi等提出的松弛维数分解(Relaxed dimensionalfactorization, RDF)预条件子进一步推广到广义鞍点问题上,并称为GRDF(Generalized RDF)预条件子.该预条件子可看做是用维数分裂迭代法求解广义鞍点问题而导出的改进维数分裂(Modified dimensional split, MDS)预条件子的松弛形式, 它相比MDS预条件子更接近于系数矩阵, 因而结合Krylov子空间方法(如GMRES)有更快的收敛速度.文中分析了GRDF预处理矩阵特征值的一些性质,并用数值算例验证了新预条件子的有效性.     

求解增广线性系统的局部多分裂迭代法

为了在高性能计算机上求解增广线性系统,基于并行多分裂的两种技巧,本文提出一种局部多分裂迭代格式,给出当增广线性系统的矩阵为M-矩阵和H-矩阵时新方法的收敛性理论.并讨论预条件矩阵的特征值情形.     

鞍点问题的广义位移分裂预条件子

对于大型稀疏非Hermitian正定线性方程组,Bai等人提出了一种位移分裂预条件子(J.Comput.Math.,24(2006)539-552).本文将这种思想用到鞍点问题上并提出了一种广义位移分裂(Generalized Shift Splitting,GSS)预条件子,同时证明了该预条件子所对应分裂迭代法的无条件收敛性.最后用数值算例验证了新预条件子的有效性.     

一种求解加权Toeplitz最小二乘问题的新预条件子（英文）

本文研究加权Toeplitz最小二乘问题的快速求解算法.首先,在增广线性系统的基础上,设计了一种用于求解此类线性系统的新型简单预条件子.其次,研究了迭代法的收敛性,并证明了预条件矩阵的所有特征值均是实数且非单位特征值位于某正区间.再次,研究了预条件矩阵的特征向量分布和最小多项式的维数.最后,相关数值实验表明新型预条件子比一些已有的预条件子更有效.     

速度追踪问题中鞍点系统的新分裂迭代

 有限元离散一类速度追踪问题后得到具有鞍点结构的线性系统,针对该鞍点系统,本文提出了一种新的分裂迭代技术.证明了新的分裂迭代方法的无条件收敛性,详细分析了新的分裂预条件子对应的预处理矩阵的谱性质.数值结果验证了对于大范围的网格参数和正则参数,新的分裂预条件子在求解有限元离散速度追踪问题得到的鞍点系统时的可行性和有效性.     

多项式分裂可行问题

本文研究多项式分裂可行问题,即由多项式不等式定义的分裂可行问题,包括凸与非凸、可行与不可行的问题;给出多项式分裂可行问题解集的半定松弛表示;研究其半定松弛化问题的性质;并基于这些性质建立求解多项式分裂可行问题的半定松弛算法.本文在较为一般的条件下证明了,如果分裂可行问题有解,则可通过本文建立的算法求得一个解点;如果问题无解,则该算法能够判别问题不可行.最后通过数值实验对算法进行验证.     

速度追踪问题产生的鞍点系统的新的分裂迭代技术

有限元离散一类速度追踪问题后得到具有鞍点结构的线性系统,针对该鞍点系统,本文提出了一种新的分裂迭代技术.证明了新的分裂迭代方法的无条件收敛性,详细分析了新的分裂预条件子对应的预处理矩阵的谱性质.数值结果验证了对于大范围的网格参数和正则参数,新的分裂预条件子在求解有限元离散速度追踪问题得到的鞍点系统时的可行性和有效性.     

广义鞍点问题基于PSS的约束预条件子

对于(1,1)块为非Hermitian阵的广义鞍点问题,本文给出了一种基于正定和反对称分裂(Positive definite andskew-Hermitian splitting, PSS)的约束预条件子.该预条件子的(1,1)块由求解非Hermitian正定线性方程组时的PSS迭代法所构造得到.文中分析了PSS约束预条件子的一些性质并证明了预处理迭代法的收敛性.最后用数值算例验证了该预条件子的有效性.     

H-矩阵线性方程组的一类预条件并行多分裂SOR迭代法

基于并行多分裂算法的思想及SOR迭代格式, 本文提出一种求解H-矩阵线性方程组新的并行多分裂SOR迭代法, 新方法某种程度上避免了SOR迭代法中选取最优参数的困难. 同时, 选取Kohno等(1997)提出的预条件子$P=I+S_{\alpha}$对原始线性方程组进行预处理, 进而给出了一种实用的预条件并行多分裂SOR迭代法. 理论分析和数值实验均表明, 新算法是实用而有效的.     

一类预条件AOR迭代法的比较定理(英文)

本文研究了当线性方程组的系数矩阵是严格对角占优L-矩阵时带有预条件子P1→kα的预条件AOR迭代方法.利用矩阵分裂的相关理论,获得了预条件AOR迭代法的收敛性结论以及参数α和k对收敛速度影响的比较定理.结果表明当α和k取值较大时这类预条件方法更加有效.文中的结论推广了Li等人关于预条件Gauss-Seidel迭代法的相关结论.最后,用数值例子进一步验证了这些结果.     

稀疏近似逆与多层块ILU预条件技术

设计了一种求解一般稀疏线性方程组的健壮且有效的可并行化预条件子,这种预条件子涉及在多层块ILU预条件子(BILUM)中使用稀疏近似逆(AINV)技术．所得的预条件子保持了BILUM的健壮性,它比标准的BILUM预条件子有两点优势:控制稀疏性的能力和增强了并行性．数值例子显示了新预条件子的有效性和效率．     

Bilus: A block version of ilus factorization

ILUS factorization has many desirable properties such as its amenability to the skyline format, the ease with which stability may be monitored, and the possibility of constructing a preconditioner with symmetric structure. In this paper we introduce a new preconditioning technique for general sparse linear systems based on the ILUS factorization strategy. The resulting preconditioner has the same properties as the ILUS preconditioner. Some theoretical properties of the new preconditioner are discussed and numerical experiments on test matrices from the Harwell-Boeing collection are tested. Our results indicate that the new preconditioner is cheaper to construct than the ILUS preconditioner.     

A new relaxed HSS preconditioner for saddle point problems

We present a preconditioner for saddle point problems. The proposed preconditioner is extracted from a stationary iterative method which is convergent under a mild condition. Some properties of the preconditioner as well as the eigenvalues distribution of the preconditioned matrix are presented. The preconditioned system is solved by a Krylov subspace method like restarted GMRES. Finally, some numerical experiments on test problems arisen from finite element discretization of the Stokes problem are given to show the effectiveness of the preconditioner.     

An alternating preconditioner for saddle point problems

Based on matrix splittings, a new alternating preconditioner with two parameters is proposed for solving saddle point problems. Some theoretical analyses for the eigenvalues of the associated preconditioned matrix are given. The choice of the parameters is considered and the quasi-optimal parameters are obtained. The new preconditioner with these quasi-optimal parameters significantly improves the convergence rate of the generalized minimal residual (GMRES) iteration. Numerical experiments from the linearized Navier-Stokes equations demonstrate the efficiency of the new preconditioner, especially on the larger viscosity parameter ν. Further extensions of the preconditioner to generalized saddle point matrices are also checked.     

On the preconditioning of the block tridiagonal linear system of equations

Two algorithms for computing the inverse factors of general tridiagonal and pentadiagonal matrices are obtained. Then, these algorithms are used for computing a block ILU preconditioner for the block tridiagonal linear system of equations. Some numerical results are given to show the robustness and efficiency of the preconditioner. The performance of the proposed preconditioner is compared with a recently proposed method.     

Uniform preconditioners for a parameter dependent saddle point problem with application to generalized Stokes interface equations

We consider an abstract parameter dependent saddle-point problem and present a general framework for analyzing robust Schur complement preconditioners. The abstract analysis is applied to a generalized Stokes problem, which yields robustness of the Cahouet-Chabard preconditioner. Motivated by models for two-phase incompressible flows we consider a generalized Stokes interface problem. Application of the general theory results in a new Schur complement preconditioner for this class of problems. The robustness of this preconditioner with respect to several parameters is treated. Results of numerical experiments are given that illustrate robustness properties of the preconditioner.     

A relaxed positive semi-definite and skew-Hermitian splitting preconditioner for non-Hermitian generalized saddle point problems

For non-Hermitian saddle point linear systems, Pan, Ng and Bai presented a positive semi-definite and skew-Hermitian splitting (PSS) preconditioner (Pan et al. Appl. Math. Comput. 172, 762–771 2006), to accelerate the convergence rate of the Krylov subspace iteration methods like the GMRES method. In this paper, a relaxed positive semi-definite and skew-Hermitian (RPSS) splitting preconditioner based on the PSS preconditioner for the non-Hermitian generalized saddle point problems is considered. The distribution of eigenvalues and the form of the eigenvectors of the preconditioned matrix are analyzed. Moreover, an upper bound on the degree of the minimal polynomial is also studied. Finally, numerical experiments of a model Navier-Stokes equation are presented to illustrate the efficiency of the RPSS preconditioner compared to the PSS preconditioner, the block diagonal preconditioner (BD), and the block triangular preconditioner (BT) in terms of the number of iteration and computational time.     

BTTB preconditioners for BTTB least squares problems

In this paper, we consider solving the least squares problem min_x‖b-Tx‖₂ by using preconditioned conjugate gradient (PCG) methods, where T is a large rectangular matrix which consists of several square block-Toeplitz-Toeplitz-block (BTTB) matrices and b is a column vector. We propose a BTTB preconditioner to speed up the PCG method and prove that the BTTB preconditioner is a good preconditioner. We then discuss the construction of the BTTB preconditioner. Numerical examples, including image restoration problems, are given to illustrate the efficiency of our BTTB preconditioner. Numerical results show that our BTTB preconditioner is more efficient than the well-known Level-1 and Level-2 circulant preconditioners.