Retooling the method of block conjugate gradients
Electronic transactions on numerical analysis, Tome 12 (2001), pp. 216-233
Block implementations of the conjugate-gradients method for the solution of a linear system must deal with linear dependences that may appear in the descent or residual blocks in the course of the iteration. New algorithms presented here avoid rank estimation and deflation through the use of changes of bases and algorithmic reformulations that eliminate rank near defects. The transformations include a robust process of nonunitary orthogonalisation in the metric of a symmetric positive-definite matrix.
Classification :
65F10, 65F25
Keywords: linear equations, conjugate gradients, orthogonalisation
Keywords: linear equations, conjugate gradients, orthogonalisation
@article{ETNA_2001__12__a1,
author = {Dubrulle, A.A.},
title = {Retooling the method of block conjugate gradients},
journal = {Electronic transactions on numerical analysis},
pages = {216--233},
year = {2001},
volume = {12},
zbl = {0985.65021},
language = {en},
url = {http://geodesic.mathdoc.fr/item/ETNA_2001__12__a1/}
}
Dubrulle, A.A. Retooling the method of block conjugate gradients. Electronic transactions on numerical analysis, Tome 12 (2001), pp. 216-233. http://geodesic.mathdoc.fr/item/ETNA_2001__12__a1/