algorithm. The first result states that, provided the normal equations associated with the
initial vectors are numerically nonsingular, the loss of orthogonality of the vectors computed
by the classical Gram-Schmidt algorithm depends quadratically on the condition number of
the initial vectors. The second result states that, provided the initial set of vectors has
numerical full rank, two iterations of the classical Gram-Schmidt algorithm are enough for …