Preconditioned conjugate gradient method