凸二次函数优化问题在快速多极边界元法中的应用

利用快速多极边界元法(FMM-BEM)求解大规模工程问题最终结为稀疏线性方程组的求解,因此,采用更好的方法求解线性方程组可以提高边界元法的计算效率.本文利用最优化数值技术处理,将稀疏线性方程组的求解等价为求解一个凸二次函数极小化的问题,并利用最优化理论及相关数学理论证明了其解的存在唯一性,为该理论的形成和发展奠定了理论基础.     

不等式约束优化一个可行序列线性方程组算法

提出了求解非线性不等式约束优化问题的一个可行序列线性方程组算法.在每次迭代中,可行下降方向通过求解两个线性方程组产生,系数矩阵具有较好的稀疏性.在较为温和的条件下,算法具有全局收敛性和强收敛性,数值试验表明算法是有效的.     

求解非对称线性方程组Lanczos算法的改进

考虑线性方程组求解问题这里A是大型稀疏、非对称和不定的可逆阵。求解问题(1)的双边Lanczos算法为算法1     

反五对角与拟反五对角方程组的追赶法

本文研究了反五对角和拟五对角线性方程组的求解问题.利用矩阵分解方法以及将系数矩阵A分解成三个简单矩阵的乘积A=LUD,获得了反五对角线性方程组以及拟反五对角线性方程组的追赶法,从而推广了对角型线性方程组追赶法.     

一种带参数的不完全Cholesky分解共轭梯度法

共轭梯度法在解高阶稀疏线性方程组方面有许多其它经典的迭代法所没有的优点,但当线性方程组相当病态、系数矩阵条件数很坏时,共轭梯度法的收敛速度很慢.因此,又产生了预条件处理共轭梯度法. 我们用预条件处理共轭梯度法求解线性方程组Ax=b(这里A是对称正定稀疏阵且条件数很大).预条件处理共轭梯度法旨在寻找一适当的正定矩阵C,C通常写成     

三维泊松方程边值问题3次Lagrange有限元方程的病态结构和最优预条件子

正1引言在科学计算和工程应用中,偏微分方程大规模数值求解问题通常转化为病态(高条件数)的大规模稀疏线性方程组的求解问题,其条件数(病态)经常随着问题规模的增加而增加[1],成为影响求解效率和精度的瓶颈因素,因此,在求解之前,使用预处理技术来减少方程组的病态,成为提高求解效率和精度的必要措施.所谓"预处理技术"是指在求解方程组     

关于非Hermitian正定线性代数方程组的超松弛HSS方法

本文针对求解大型稀疏非Hermitian正定线性方程组的HSS迭代方法,利用迭代法的松弛技术进行加速,提出了一种具有三个参数的超松弛HSS方法(SAHSS)和不精确的SAHSS方法(ISAHSS),它采用CG和一些Krylov子空间方法作为其内部过程,并研究了SAHSS和ISAHSS方法的收敛性.数值例子验证了新方法的有效性.     

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

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

稀疏线性方程组的一种新解法及其应用

本方给出了求解稀疏线性方程组的一种新方法,称作松弛求解法。该法是化学方程式的一种新的高效全能配平法的理论基础。     

求解稀疏分裂可行问题的一种投影算法

本文研究了稀疏分裂可行问题.通过将分裂可行问题转化为一个目标函数为凸函数的稀疏约束优化问题,设计一种梯度投影算法来求解此问题,获得了算法产生的点列可以收敛到稀疏分裂可行问题的一个解.用数值例子说明了算法的有效性.     

Preconditioning of two-by-two block matrix systems with square matrix blocks,with applications

Two-by-two block matrices of special form with square matrix blocks arise in important applications, such as in optimal control of partial differential equations and in high order time integration methods.Two solution methods involving very efficient preconditioned matrices, one based on a Schur complement reduction of the given system and one based on a transformation matrix with a perturbation of one of the given matrix blocks are presented. The first method involves an additional inner solution with the pivot matrix block but gives a very tight condition number bound when applied for a time integration method. The second method does not involve this matrix block but only inner solutions with a linear combination of the pivot block and the off-diagonal matrix blocks.Both the methods give small condition number bounds that hold uniformly in all parameters involved in the problem, i.e. are fully robust. The paper presents shorter proofs, extended and new results compared to earlier publications.     

Structured preconditioners for nonsingular matrices of block two-by-two structures

For the large sparse block two-by-two real nonsingular matrices, we establish a general framework of practical and efficient structured preconditioners through matrix transformation and matrix approximations. For the specific versions such as modified block Jacobi-type, modified block Gauss-Seidel-type, and modified block unsymmetric (symmetric) Gauss-Seidel-type preconditioners, we precisely describe their concrete expressions and deliberately analyze eigenvalue distributions and positive definiteness of the preconditioned matrices. Also, we show that when these structured preconditioners are employed to precondition the Krylov subspace methods such as GMRES and restarted GMRES, fast and effective iteration solvers can be obtained for the large sparse systems of linear equations with block two-by-two coefficient matrices. In particular, these structured preconditioners can lead to efficient and high-quality preconditioning matrices for some typical matrices from the real-world applications.

     

COMBINATIVE PRECONDITIONERS OF MODIFIED INCOMPLETE CHOLESKY FACTORIZATION AND SHERMAN-MORRISON-WOODBURY UPDATE FOR SELF-ADJOINT ELLIPTIC DIRICHLET-PERIODIC BOUNDARY VALUE PROBLEMS

For the system of linear equations arising from discretization of the second-order self-adjoint elliptic Dirichlet-periodic boundary value problems,by making use of the specialstructure of the coefficient matrix we present a class of combinative preconditioners whichare technical combinations of modified incomplete Cholesky factorizations and Sherman-Morrison-Woodbury update.Theoretical analyses show that the condition numbers of thepreconditioned matrices can be reduced to(?)(h~(-1)),one order smaller than the conditionnumber(?)(h~(-2))of the original matrix.Numerical implementations show that the resultingpreconditioned conjugate gradient methods are feasible,robust and efficient for solving thisclass of linear systems.     

Rotated block triangular preconditioning based on PMHSS

Based on the PMHSS preconditioning matrix, we construct a class of rotated block triangular preconditioners for block two-by-two matrices of real square blocks, and analyze the eigen-properties of the corresponding preconditioned matrices. Numerical experiments show that these rotated block triangular preconditioners can be competitive to and even more efficient than the PMHSS pre-conditioner when they are used to accelerate Krylov subspace iteration methods for solving block two-by-two linear systems with coefficient matrices possibly of nonsymmetric sub-blocks.     

Kronecker product approximation preconditioners for convection–diffusion model problems

We consider the iterative solution of linear systems arising from four convection–diffusion model problems: scalar convection–diffusion problem, Stokes problem, Oseen problem and Navier–Stokes problem. We design preconditioners for these model problems that are based on Kronecker product approximations (KPAs). For this we first identify explicit Kronecker product structure of the coefficient matrices, in particular for the convection term. For the latter three model cases, the coefficient matrices have a 2 × 2 block structure, where each block is a Kronecker product or a summation of several Kronecker products. We then use this structure to design a block diagonal preconditioner, a block triangular preconditioner and a constraint preconditioner. Numerical experiments show the efficiency of the three KPA preconditioners, and in particular of the constraint preconditioner that usually outperforms the other two. This can be explained by the relationship that exists between these three preconditioners: the constraint preconditioner can be regarded as a modification of the block triangular preconditioner, which at its turn is a modification of the block diagonal preconditioner based on the cell Reynolds number. Copyright © 2009 John Wiley & Sons, Ltd.     

Additive block diagonal preconditioning for block two-by-two linear systems of skew-Hamiltonian coefficient matrices

For a class of block two-by-two systems of linear equations with certain skew-Hamiltonian coefficient matrices, we construct additive block diagonal preconditioning matrices and discuss the eigen-properties of the corresponding preconditioned matrices. The additive block diagonal preconditioners can be employed to accelerate the convergence rates of Krylov subspace iteration methods such as MINRES and GMRES. Numerical experiments show that MINRES preconditioned by the exact and the inexact additive block diagonal preconditioners are effective, robust and scalable solvers for the block two-by-two linear systems arising from the Galerkin finite-element discretizations of a class of distributed control problems.     

Block ILU Preconditioners for a Nonsymmetric Block-Tridiagonal M-Matrix

We propose block ILU (incomplete LU) factorization preconditioners for a nonsymmetric block-tridiagonal M-matrix whose computation can be done in parallel based on matrix blocks. Some theoretical properties for these block ILU factorization preconditioners are studied and then we describe how to construct them effectively for a special type of matrix. We also discuss a parallelization of the preconditioner solver step used in nonstationary iterative methods with the block ILU preconditioners. Numerical results of the right preconditioned BiCGSTAB method using the block ILU preconditioners are compared with those of the right preconditioned BiCGSTAB using a standard ILU factorization preconditioner to see how effective the block ILU preconditioners are.     

无界分块算子矩阵的谱分布

本文首先给出次对角元有界的2×2阶无界算子矩阵的Gershgorin定理,然后利用主对角元算子的谱和数值域刻画整个算子矩阵的谱分布.特别地,当次对角元算子互为共轭(反共轭)算子时,结合二次数值域和Gershgorin定理对谱分布给出更精细的描述.     

Band Toeplitz preconditioners for block Toeplitz systems

We consider the solutions of block Toeplitz systems with Toeplitz blocks by the preconditioned conjugate gradient (PCG) method. Here the block Toeplitz matrices are generated by nonnegative functions f(x,y). We use band Toeplitz matrices as preconditioners. The generating functions g(x,y) of the preconditioners are trigonometric polynomials of fixed degree and are determined by minimizing (f − g)/f∞. We prove that the condition number of the preconditioned system is O(1). An a priori bound on the number of iterations for convergence is obtained.     

Preconditioning techniques for an image deblurring problem

In this paper, we consider the solution of a large linear system of equations, which is obtained from discretizing the Euler–Lagrange equations associated with the image deblurring problem. The coefficient matrix of this system is of the generalized saddle point form with high condition number. One of the blocks of this matrix has the block Toeplitz with Toeplitz block structure. This system can be efficiently solved using the minimal residual iteration method with preconditioners based on the fast Fourier transform. Eigenvalue bounds for the preconditioner matrix are obtained. Numerical results are presented. Copyright © 2016 John Wiley & Sons, Ltd.