It would be nice to have a CGResult struct that saves things like the last residual, number of iterations needed to converge, and options. If this is approved in principle, I can make a PR. This will be a breaking change however. A similar thing may also need to be done for the other algorithms but I am mostly interested in CG now.