Modified Gram-Schmidt for solving linear least squares problems is equivalent to Gaussian elimination for the normal equations H. Späth (1988) Applicationes Mathematicae