Change `itmax` behaviour of non-linear solver

[MarcoArtiano]

At the moment the krylov solver prescribe a default itmax = 2n where n is the dimension of the right hand side vector. We should inform DiffEq of the convergence rather than setting itmax.

See #31 comments