Add Jacobian printing at end of solve
Add parameters to handle the cholesky correction