Reorthogonalized block classical Gram-Schmidt

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 - Q^T 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.