Bi-Newton's method for computing spectral projectors

An efficient Newton-like method for computing the spectral projector associated with a separated group of eigenvalues near a specified shift of a large sparse matrix is proposed and justified. A number of numerical experiments with a discrete analogue of the non-Hermitian elliptic operator are discussed.