Reorthogonalized block classical Gram-Schmidt

Jesse L. Barlow, Alicja Smoktunowicz

Research output: Contribution to journalArticlepeer-review

11 Scopus citations

Abstract

A reorthogonalized block classical Gram-Schmidt algorithm is proposed that factors a full column rank matrix A into A = QR where Q is left orthogonal (has orthonormal columns) and R is upper triangular and nonsingular. This block Gram-Schmidt algorithm can be implemented using matrix-matrix operations making it more efficient on modern architectures than orthogonal factorization algorithms based upon matrix-vector operations and purely vector operations. Gram-Schmidt orthogonal factorizations are important in the stable implementation of Krylov space methods such as GMRES and in approaches to modifying orthogonal factorizations when columns and rows are added or deleted from a matrix. With appropriate assumptions about the diagonal blocks of R, the algorithm, when implemented in floating point arithmetic with machine unit ε M, produces Q and R such that {double pipe}I - QT Q{double pipe} = O(εM) and {double pipe}A - QR{double pipe} = O(εM{double pipe}A{double pipe}). The first of these bounds has not been shown for a block Gram-Schmidt procedure before. As consequence of these results, we provide a different analysis, with a slightly different assumption, that re-establishes a bound of Giraud et al. (Num Math, 101(1):87-100, 2005) for the CGS2 algorithm.

Original languageEnglish (US)
Pages (from-to)395-423
Number of pages29
JournalNumerische Mathematik
Volume123
Issue number3
DOIs
StatePublished - Mar 2013

All Science Journal Classification (ASJC) codes

  • Computational Mathematics
  • Applied Mathematics

Fingerprint Dive into the research topics of 'Reorthogonalized block classical Gram-Schmidt'. Together they form a unique fingerprint.

Cite this