A generalization of -step variants of gradient methods

详细信息	查看全文 | 推荐本文 |

• 作者：J.A. Alvarez-Dios^a ; ^{joseantonio.alvarez.dios@usc.es} ; J.C. Cabaleiro^b ; ^{jc.cabaleiro@usc.es} ; G. Casal^c ; ^{xerardo.casal@usc.es}
• 关键词：Iterative methods ; s-step ; Large linear systems ; Krylov subspace
• 刊名：Journal of Computational and Applied Mathematics
• 出版年：2012
• 期刊代码：72_03770427
• 类别：cp
• 出版时间：June, 2012
• 卷：236
• 期：12
• 页码：2938-2953
• 文件大小：348 K

摘要

The -step methods were proposed by Chronopoulos to gain efficiency in parallel programming of iterative methods for linear systems. They are variants of classical iterative methods based on the construction of a Krylov subspace basis on each iteration. These s-step methods were inferred from algorithms like the Conjugate Gradient, Generalized Conjugate Residual or the Minimal Residual. They converge for all symmetric, nonsymmetric definite and some nonsymmetric indefinite matrices. In this paper, we introduce an -step variant of a General Orthogonalization Algorithm, that is, a generalization of s-step variants of gradient methods. We prove convergence and obtain error estimates. We also describe an Orthomin variant, together with a convergence theorem. From this we derive the well known -step methods as particular cases, and some which are newfound to our knowledge. This provides a unified framework to derive and study -step methods. Some of the methods obtained are convergent for every nonsingular matrix. Finally, we give some numerical results for the new proposed methods, showing that the parallel implementation of these overcomes the original ones.