This class implements a conjugate gradient matrix solver. Preconditioning is implemented using a Jacobi (diagonal) preconditioner.