Combining the preconditioned conjugate gradient method and a matrix iterative method