Abstract
The solution of large linear systems with multiple right-hand sides given simultaneously is required in many largescale scientific and engineering applications modelled by either partial differential or boundary integral equations. In this paper we introduce variants of the block GMRES method that combine initial deflation and eigenvalue recycling strategies to remedy some typical convergence problems of block Krylov solvers.