Set default initial guess equal to the right-hand side
Replace recomputation of dot product with stored value
Extend resnorm to work with a system matrix passed as
a sparse matrix or as a function performing a linear
operation