In this paper we derive convergence estimates for the iterative solution of nonsymmetric linear systems by GMRES. We work in the context of strongly convergent-collectively compact sequences of approximations to linear compact xed point problems. Our estimates are intended to explain the observations that the performance of GMRES is independent of the discretization if the resolution of the dis...

Over the past several years a number of hybrid algorithms have been proposed for solving large sparse systems of linear algebraic equations. In this paper we consider the hybrid GMRES algorithm by Nachtigal, Reichel and Trefethen (1992) and show that in the floating-point arithmetic there exist some cases in which the properties of this algorithm are lost, e.g., the result is false, or the coef...

In this work, a two-layer coupled fluid-structure-structure interaction model is considered, which incorporates an anisotropic structure model into the fluid-structure interaction problems. We propose two domain decomposition solvers for such a class of coupled problems: a Robin-Robin preconditioned GMRES solver combined with an inner Dirichlet-Neumann iterative solver, and a Robin-Robin precon...

The GMRES method for solving nonsymmetric linear equations is generally used with restarting to reduce storage and orthogonalization costs. Restarting slows down the convergence. However, it is possible to save some important information at the time of the restart. It is proposed that approximate eigenvectors corresponding to a few of the smallest eigenvalues be formed and added to the subspace...

We look at multigrid methods for unsteady viscous compressible flows. We specifically target smoothers that can be used in parallel and without computation of a Jacobian, which are particlarly attractive candidates in the context of Discontinuous Galerkin discretizations. In CFD, a plethora of nonlinear smoothers have been suggested which are hard to analyze. Our methodology is to use a linear ...

and Kk(BA,Br) = span{Br, (BA)Br, . . . , (BA)k−1Br}, (3) where B ∈ Rn×m is the mapping and preconditioning matrix, and apply Krylov subspace iteration methods on these subspaces. For overdetermined problems, applying the standard CG method to Kk(BA,Br) leads to the preconditioned CGLS [3] or CGNR [9] method while for underdetermined problems it leads to preconditioned CGNE [9] method. The GMRES...

In this paper, we propose two new algorithms based on modified global Arnoldi algorithm for solving large Sylvester matrix equations AX + XB = C where A ∈ Rn×n, B ∈ Rs×s, X and C ∈ Rn×s. These algorithms are based on the global FOM and GMRES algorithms and we call them by Global FOM-SylvesterLike(GFSL) and Global GMRES-Sylvester-Like(GGSL) algorithms, respectively. Some theoretical results and ...