Preconditioned iterative methods for a class of nonlinear eigenvalue problems

This paper proposes new iterative methods for the efficient computation of the smallest eigenvalue of symmetric nonlinear matrix eigenvalue problems of large order with a monotone dependence on the spectral parameter. Monotone nonlinear eigenvalue problems for differential equations have important applications in mechanics and physics. The discretization of these eigenvalue problems leads to nonlinear eigenvalue problems with very large sparse ill-conditioned matrices monotonically depending on the spectral parameter. To compute the smallest eigenvalue of large-scale matrix nonlinear eigenvalue problems, we suggest preconditioned iterative methods: preconditioned simple iteration method, preconditioned steepest descent method, and preconditioned conjugate gradient method. These methods use only matrix-vector multiplications, preconditioner-vector multiplications, linear operations with vectors, and inner products of vectors. We investigate the convergence and derive grid-independent error estimates for these methods. Numerical experiments demonstrate the practical effectiveness of the proposed methods for a model problem.     

关于鞍点问题的预处理HSS-SOR交替分裂迭代方法

为了高效地求解大型稀疏鞍点问题,在白中治,Golub和潘建瑜提出的预处理对称/反对称分裂(PHss)迭代法的基础上,通过结合SOR-like迭代格式对原有迭代算法进行加速,提出了一种预处理HSS-SOR交替分裂迭代方法,并研究了该算法的收敛性.数值例子表明:通过参数值的选择,新算法比SOR-like和PHSS算法都具有更快的收敛速度和更少的迭代次数,选择了合适的参数值后,可以提高算法的收敛效率.     

Multilevel iterative methods for mixed finite element discretizations of elliptic problems

Summary For solving second order elliptic problems discretized on a sequence of nested mixed finite element spaces nearly optimal iterative methods are proposed. The methods are within the general framework of the product (multiplicative) scheme for operators in a Hilbert space, proposed recently by Bramble, Pasciak, Wang, and Xu [5,6,26,27] and make use of certain multilevel decomposition of the corresponding spaces for the flux variable.     

Preconditioned conjugate gradient methods for large-scale fluid flow applications

The simulation of large-scale fluid flow applications often requires the efficient solution of extremely large nonsymmetric linear and nonlinear sparse systems of equations arising from the discretization of systems of partial differential equations. While preconditioned conjugate gradient methods work well for symmetric, positive-definite matrices, other methods are necessary to treat large, nonsymmetric matrices. The applications may also involve highly localized phenomena which can be addressed via local and adaptive grid refinement techniques. These local refinement methods usually cause non-standard grid connections which destroy the bandedness of the matrices and the associated ease of solution and vectorization of the algorithms. The use of preconditioned conjugate gradient or conjugate-gradient-like iterative methods in large-scale reservoir simulation applications is briefly surveyed. Then, some block preconditioning methods for adaptive grid refinement via domain decomposition techniques are presented and compared. These techniques are being used efficiently in existing large-scale simulation codes.     

Preconditioned iterative methods for diffusion problems with high‐contrast inclusions

This paper is concerned with robust numerical treatment of an elliptic PDE with high‐contrast coefficients, for which classical finite‐element discretizations yield ill‐conditioned linear systems. This paper introduces a procedure by which the discrete system obtained from a linear finite element discretization of the given continuum problem is converted into an equivalent linear system of the saddle‐point type. Three preconditioned iterative procedures—preconditioned Uzawa, preconditioned Lanczos, and preconditioned conjugate gradient for the square of the matrix—are discussed for a special type of the application, namely, highly conducting particles distributed in the domain. Robust preconditioners for solving the derived saddle‐point problem are proposed and investigated. Robustness with respect to the contrast parameter and the discretization scale is also justified. Numerical examples support theoretical results and demonstrate independence of the number of iterations of the proposed iterative schemes on the contrast in parameters of the problem and the mesh size.     

Preconditioned iterative methods for linear discrete ill-posed problems from a Bayesian inversion perspective

In this paper we revisit the solution of ill-posed problems by preconditioned iterative methods from a Bayesian statistical inversion perspective. After a brief review of the most popular Krylov subspace iterative methods for the solution of linear discrete ill-posed problems and some basic statistics results, we analyze the statistical meaning of left and right preconditioners, as well as projected-restarted strategies. Computed examples illustrating the interplay between statistics and preconditioning are also presented.     

Preconditioned iterative algorithms for large sparse unsymmetric problems

A numerical study of the efficiency of the generalized conjugate residual methods (GCR) is performed using three different preconditioners all based upon an incomplete LU factorization. The GCR behavior is evaluated in connection with the solution of large, sparse unsymmetric systems of equations, arising from the finite element integration of the diffusion-convection equation for 2-dimensional (2-D) and 3-D problems with different Peclet and Courant numbers. The order of the test matrices ranges from 450 to 1700. Results from a set of numerical experiments are presented and comparisons with preconditioned GCR methods and with direct method are carried out.     

Preconditioned iterative methods on sparse subspaces

When some rows of the system matrix and a preconditioner coincide, preconditioned iterations can be reduced to a sparse subspace. Taking advantage of this property can lead to considerable memory and computational savings. This is particularly useful with the GMRES method. We consider the iterative solution of a discretized partial differential equation on this sparse subspace. With a domain decomposition method and a fictitious domain method the subspace corresponds a small neighborhood of an interface. As numerical examples we solve the Helmholtz equation using a fictitious domain method and an elliptic equation with a jump in the diffusion coefficient using a separable preconditioner.     

Error estimates for discretizations of coupled electromagnetic-mechanical problems

In this note, a method is outlined to obtain a priori error estimates for the finite element discretization of coupled electromagnetic mechanical problems as arise, e.g., in electromagnetic metal forming. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Preconditioned iterative methods for fractional diffusion models in finance

Most recent qualitative models for financial assets assume that the dynamics of underlying equity prices follows a jump or Lévy process. It has been evident that some most intricate characteristics of such dynamics can be captured by CGMY and KoBoL procedures. The prices of financial derivatives with such models satisfy fractional partial differential equations or partial integro‐differential equations. This study focuses at aforementioned fractional equations and discretizes them via a monotone Crank–Nicolson procedure. A spatial extrapolation strategy is introduced to ensure an overall second‐order accuracy in approximations. Preconditioned conjugate gradient normal residual methods are incorporated for solving resulted linear systems. Numerical examples are given to illustrate the accuracy and efficiency of the novel computational approaches implemented. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1382–1395, 2015     

Preconditioned sor methods for generalized least-squares problems

1.IntroductionThegeneralizedleastsquaresproblem,definedasmin(Ax--b)"W--'(Ax--b),(1.1)xacwhereAERm",m>n,bERm,andWERm'misasymmetricandpositivedefinitematrix,isfrequentlyfoundwhensolvingproblemsinstatistics,engineeringandeconomics.Forexample,wegetgeneralizedleastsquaresproblemswhensolvingnonlinearregressionanalysisbyquasi-likelihoodestimation,imagereconstructionproblemsandeconomicmodelsobtainedbythemaximumlikelihoodmethod(of.[1,21).Paige[3,4]investigatestheproblemexplicitly.Hechangestheorig…     

求解分布控制问题的预处理最小残量方法(英文)

通过分析Bai(Bai Z Z.Block preconditioners for elliptic PDE-constrained optimization problems.Computing,2011,91:379-395)给出的离散分布控制问题的块反对角预处理线性系统,提出了该问题的一个等价线性系统,并且运用带有预处理子的最小残量方法对该系统进行求解.理论分析和数值实验结果表明,所提出的预处理最小残量方法对于求解该类椭圆型偏微分方程约束最优分布控制问题非常有效,尤其当正则参数适当小的时候.     

On iterative methods for equilibrium problems

In this paper, we introduce a hybrid iterative scheme for finding a common element of the set of solutions of an equilibrium problem and the set of fixed points of finitely many nonexpansive mappings. We prove that the approximate solution converges strongly to a solution of a class of variational inequalities under some mild conditions, which is the optimality condition for some minimization problem. We also give some comments on the results of Plubtieng and Punpaeng [S. Plubtieng, R. Punpaeng, A general iterative method for equilibrium problems and fixed point problems in Hilbert spaces, J. Math. Anal. Appl. 336 (2007) 455–469]. Results obtained in this paper may be viewed as an improvement and refinement of the previously known results in this area.     

Preconditioned GAOR methods for solving weighted linear least squares problems

In this paper, we present the preconditioned generalized accelerated overrelaxation (GAOR) method for solving linear systems based on a class of weighted linear least square problems. Two kinds of preconditioning are proposed, and each one contains three preconditioners. We compare the spectral radii of the iteration matrices of the preconditioned and the original methods. The comparison results show that the convergence rate of the preconditioned GAOR methods is indeed better than the rate of the original method, whenever the original method is convergent. Finally, a numerical example is presented in order to confirm these theoretical results.     

On multilevel iterative methods for optimization problems

This paper is concerned with multilevel iterative methods which combine a descent scheme with a hierarchy of auxiliary problems in lower dimensional subspaces. The construction of auxiliary problems as well as applications to elasto-plastic model and linear programming are described. The auxiliary problem for the dual of a perturbed linear program is interpreted as a dual of perturbed aggregated linear program. Coercivity of the objective function over the feasible set is sufficient for the boundedness of the iterates. Equivalents of this condition are presented in special cases.Supported by NSF under grant DMS-8704169, AFOSR under grant 86-0126, and ONR under Contract N00014-83-K-0104. Consulting for American Airlines Decision Technologies, MD 2C55, P.O. Box 619616, DFW, TX 75261-9616, USA.Supported by NSF under grant DMS-8704169 and AFOSR under grant 86-0126.     

Preconditioned CG-type methods for solving the coupled system of fundamental semiconductor equations

This paper presents some of the authors' experimental results in applying Preconditioned CG-type methods to nonsymmetric systems of linear equations arising in the numerical solution of the coupled system of fundamental stationary semiconductor equations. For this type of problem it is shown that these iterative methods are efficient both in computation times and in storage requirements. All results have been obtained on an HP 350 computer.     

Time discretizations for evolution problems

The aim of this work is to give an introductory survey on time discretizations for liner parabolic problems. The theory of stability for stiff ordinary differential equations is explained on this problem and applied to Runge-Kutta and multi-step discretizations. Moreover, a natural connection between Galerkin time discretizations and Runge-Kutta methods together with order reduction phenomenon is discussed.     

Unified finite element discretizations of coupled Darcy–Stokes flow

In this article, we discuss some new finite element methods for flows which are governed by the linear stationary Stokes system on one part of the domain and by a second order elliptic equation derived from Darcy's law in the rest of the domain, and where the solutions in the two domains are coupled by proper interface conditions. All the methods proposed here utilize the same finite element spaces on the entire domain. In particular, we show how the coupled problem can be solved by using standard Stokes elements like the MINI element or the Taylor–Hood element in the entire domain. Furthermore, for all the methods the handling of the interface conditions are straightforward. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Chaotic iterative methods for the linear complementarity problems

A class of parallel multisplitting chaotic relaxation methods is established for the large sparse linear complementarity problems, and the global and monotone convergence is proved for the H-matrix and the L-matrix classes, respectively. Moreover, comparison theorem is given, which describes the influences of the parameters and the multiple splittings upon the monotone convergence rates of the new methods.