Add iteration control mechanism and provide default based on residuals.

This closes #44 and goes part way to closing #43.  We can add more flexibility to norms etc in 0.4.0.
Unfortunately, I'm not very familiar with functional programming. Could you sketch the implementation here?