English

On Two-Stage Householder Orthogonalization

Numerical Analysis 2026-03-24 v2 Numerical Analysis

Abstract

Two-stage orthogonalization is essential in numerical algorithms such as Krylov subspace methods. For this task we need to orthogonalize a matrix AA against another matrix VV with orthonormal columns. A common approach is to employ the block Gram--Schmidt algorithm. However, its stability largely depends on the condition number of [V,A][V,A]. While performing a Householder orthogonalization on [V,A][V,A] is unconditionally stable, it does not utilize the knowledge that VV has orthonormal columns. To address these issues, we propose a two-stage Householder orthogonalization algorithm based on the generalized Householder transformation. Instead of explicitly orthogonalizing the entire VV, our algorithm only needs to orthogonalizes a square submatrix of VV. Theoretical analysis and numerical experiments demonstrate that our method is also unconditionally stable.

Keywords

Cite

@article{arxiv.2602.14449,
  title  = {On Two-Stage Householder Orthogonalization},
  author = {Zhuang-Ao He and Meiyue Shao},
  journal= {arXiv preprint arXiv:2602.14449},
  year   = {2026}
}
R2 v1 2026-07-01T10:37:59.821Z