I find that the norm(info.residual) that GMRES prints at each step (inside a single iteration) is different from the actual residual norm ||(b-Ax)|| that one computes with the current approximate x manually. Are these 2 norms (GMRES and manual) supposed to match at the end of every iteration though?
Related topics
Topic | Replies | Views | Activity | |
---|---|---|---|---|
How to extract intermediate results from GMRES of Krylov.jl | 7 | 594 | September 12, 2022 | |
Custom dot product for GMRES | 3 | 72 | April 23, 2025 | |
Krylov subspace polynomial of one right hand side, applied to another | 0 | 285 | November 16, 2020 | |
Matrix free GMRES with structs? | 3 | 703 | April 15, 2022 | |
Convergence of gmres from IterativeSolvers in minimal examples | 7 | 329 | June 23, 2024 |