Loss and Recapture of Orthogonality in the Modified Gram–Schmidt Algorithm

From MaRDI portal





The authors derive some important numerical properties of the \(QR\) factorization of an \(m\times n\) matrix \(A\) via the modified Gram-Schmidt algorithm (\(MGS\)) in the presence of rounding errors. It is shown that the upper triangular matrix \(R\) obtained from \(MGS\) is numerically as good as that obtained from rotations or Householder transformations. The special structure of the product of the Householder transformation is derived, and then used to explain and bound the loss of orthogonality in \(MGS\). This numerical equivalence is used to show how orthogonality in \(MGS\) can be regained in general. This is illustrated by deriving a backward stable algorithm based on \(MGS\) for a class of problems which includes the solution of nonsingular linear systems, the minimum 2-norm solution of undetermined linear systems, and linear least squares problems. In the final part of the paper there are comparisons of \(MGS\) and Householder transformations. It is explained when \(MGS\) may be used in preference to the Householder \(QR\) factorization of \(A\).




Cited in
(63)








This page was built for publication: Loss and Recapture of Orthogonality in the Modified Gram–Schmidt Algorithm

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q3988983)