Note on a Householder implementation of the GMRES method

A GMRES method of solution of a linear system, Ax = b, is given. The implementation uses Householder transformations which include normalization of the Krylov vectors. This method is compared to a technique using the Gram-Schmidt process. The tests indicate that the Householder version of GMRES is susceptible to numerical error. 6 refs., 1 fig., 2 tabs. (WRF)