Is there a way to extract the Krylov Space?

[Veenty]

Given a linear operator A, right hand side b, and initial guess x0. I would like to extract the Krylov space that is generated in gmres(A, b, x0). Is there any easy way of doing this?

[amontoison]

Hi @Veenty ! Yes, you can recover the Krylov basis with the in-place version of gmresand workspace.V. You can also call the Arnoldi process with V, β, H = arnoldi(A, b - A * x0, k) (without a preconditioner).

[amontoison]

@Veenty May I close the issue? Should I add something in the documentation to specify the connection between the method and the Krylov processes / subspaces?