An augmented subspace projection method for eigenvalue problems