Make sure rho is updated prior to return
Check that the reduced discriminant is strictly positive in truncated conjugate gradient