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 | |
---|---|---|---|---|
Unbounded IterativeSolvers.Gmres! different output to MATLAB, slows convergence significantly | 13 | 1876 | March 7, 2019 | |
Linear system problems | 35 | 2493 | February 4, 2021 | |
Custom dot product for GMRES | 14 | 260 | April 28, 2025 | |
How to extract intermediate results from GMRES of Krylov.jl | 7 | 622 | September 12, 2022 | |
Convergence of gmres from IterativeSolvers in minimal examples | 7 | 395 | June 23, 2024 |