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?