Perform warm-start outside the Krylov methods

[amontoison]

I should investigate whether we can implement warm-start directly in krylov_solve.jl.
This approach would be more robust because it would shift the right-hand side before invoking the Krylov method on the shifted linear system.
After solving, the solution would then be updated by adding the initial solution to the correction.

To be investigated...