An automatic procedure for updating the block size in the block conjugate gradient method for solving linear systems (Q1405653)

The paper considers the problem of constructing an efficient automatic procedure for reducing the block size in the block conjugate gradient method insuring that the resulting rate of convergence is comparable with that of the block conjugate gradient method with constant block size. The numerical results provided show that, independently of the type of distribution of the smallest eigenvalues of the preconditioned matrix, the procedure suggested always leads to a decrease of the arithmetic costs with respect to those of the block method with constant block size.