Breakdown-free GMRES for Singular Systems

From MaRDI portal





The generalized minimal residual (GMRES) method is based on the Arnoldi process, which may suffer from benign or hard breakdown. First the authors investigate in which space the solution of the system (or the least squares solution in the case of an inconsistent system) \(Ax=b\) can be found. If \({\mathcal K}_{n-1}\) is the Krylov subspace at hard breakdown, then the solution is in \({\mathcal K}_{n-1}+{\mathcal N}(A^p)\), which requires the inclusion of an eigenvector of \(A^p\) corresponding to the eigenvalue 0, to extend the Krylov sequence. This is the basis of their algorithm to overcome hard breakdown. Another breakdown-free variant of the range restricted version was given by \textit{D. Calvetti, B. Lewis}, and \textit{L. Reichel} [Linear Algebra Appl. 316, No.~1--3, 157--169 (2000; Zbl 0963.65042)]. The present method is also related to, but different from the method proposed by \textit{Q. Ye} [Math. Comput. 62, No.~205, 179--207 (1994; Zbl 0796.65046)].




Cited in
(47)








This page was built for publication: Breakdown-free GMRES for Singular Systems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5317638)