A Deflation Technique for Linear Systems of Equations

Iterative methods for solving linear systems of equations can be very eecient if the structure of the coeecient matrix can be exploited to accelerate the convergence of the iterative process. However, for classes of problems for which suitable preconditioners cannot be found or for which the iteration scheme does not converge, iterative techniques may be inappropriate. This paper proposes a technique for deeating the eigenvalues, and associated eigenvectors, of the iteration matrix which either slow down convergence or cause divergence. This process is completely general and works by approximating the eigenspace I P corresponding to the unstable or slowly converging modes and then applying a coupled iteration scheme on I P and its orthogonal complementQ: