Generalized block tuned preconditioners for SPD eigensolvers

Given an n x n symmetric positive definite (SPD) matrix A and an SPD preconditioner P we propose a new class of generalized block tuned (GBT) preconditioners. These are defined as a p-rank correction of P with the property that arbitrary (positive) parameters are eigenvalues of the preconditioned matrix. We propose to employ these GBT preconditioners to accelerate the iterative solution of linear systems like (A - c I) s = r in the framework of iterative eigensolvers. We give theoretical evidence that a suitable, and effective, choice of the scalars is able to shift p eigenvalues of P (A - c I) very close to one. Numerical experiments on various matrices of very large size show that the proposed preconditioner is able to yield an almost constant number of iterations, for different eigenpairs, irrespective of the relative separation between consecutive eigenvalues. We also give numerical evidence that the GBT preconditioner is always far superior to the spectral preconditioner on matrices with highly clustered eigenvalues.