Hi list, I was wondering the reason why classical gram schmidt is used in the orthogonalizations in the gmres implementation as default? As far as I remember, this was unstable numerically. Best, Umut