English

A stable one-synchronization variant of reorthogonalized block classical Gram--Schmidt

Numerical Analysis 2025-06-06 v2 Numerical Analysis

Abstract

The block classical Gram--Schmidt (BCGS) algorithm and its reorthogonalized variant are widely-used methods for computing the economic QR factorization of block columns XX due to their lower communication cost compared to other approaches such as modified Gram--Schmidt and Householder QR. To further reduce communication, i.e., synchronization, there has been a long ongoing search for a variant of reorthogonalized BCGS variant that achieves O(u)O(u) loss of orthogonality while requiring only \emph{one} synchronization point per block column, where uu represents the unit roundoff. Utilizing Pythagorean inner products and delayed normalization techniques, we propose the first provably stable one-synchronization reorthogonalized BCGS variant, demonstrating that it has O(u)O(u) loss of orthogonality under the condition O(u)κ2(X)1/2O(u) \kappa^2(X) \leq 1/2, where κ()\kappa(\cdot) represents the condition number. By incorporating one additional synchronization point, we develop a two-synchronization reorthogonalized BCGS variant which maintains O(u)O(u) loss of orthogonality under the improved condition O(u)κ(X)1/2O(u) \kappa(X) \leq 1/2. An adaptive strategy is then proposed to combine these two variants, ensuring O(u)O(u) loss of orthogonality while using as few synchronization points as possible under the less restrictive condition O(u)κ(X)1/2O(u) \kappa(X) \leq 1/2. As an example of where this adaptive approach is beneficial, we show that using the adaptive orthogonalization variant, ss-step GMRES achieves a backward error comparable to ss-step GMRES with BCGSI+, also known as BCGS2, both theoretically and numerically, but requires fewer synchronization points.

Cite

@article{arxiv.2411.07077,
  title  = {A stable one-synchronization variant of reorthogonalized block classical Gram--Schmidt},
  author = {Erin Carson and Yuxin Ma},
  journal= {arXiv preprint arXiv:2411.07077},
  year   = {2025}
}

Comments

25 pages, 8. figures

R2 v1 2026-06-28T19:55:41.682Z