The restarted and augmented GMRES method is considered. Harmonic Ritz vectors are added to the Krylov subspace.
Estimates for residual norm and sufficient condition for convergence are presented. Numerical results demonstrate that the new version is always better than the standard GMRES.