APPLICATION OF A BLOCK MODIFIED CHEBYSHEV ALGORITHM TO THE ITERATIVE SOLUTION OF SYMMETRICAL LINEAR-SYSTEMS WITH MULTIPLE RIGHT-HAND SIDE VECTORS

An adaptive Richardson iteration method is described for the solution of large sparse symmetric positive definite linear systems of equation s with multiple right-hand side vectors. This scheme ''learns'' about the linear system to be solved by computing inner products of residual matrices during the iterations. These inner products are interpreted as block modified moments. A block version of the modified Chebyshev a lgorithm is presented which yields a block tridiagonal matrix from the block modified moments and the recursion coefficients of the residual polynomials. The eigenvalues of this block tridiagonal matrix define an interval, which determines the choice of relaxation parameters for Richardson iteration. Only minor modifications are necessary in order to obtain a scheme for the solution of symmetric indefinite linear sys tems with multiple right-hand side vectors. We outline the changes req uired.