On Two-Stage Householder Orthogonalization
Abstract
Two-stage orthogonalization is essential in numerical algorithms such as Krylov subspace methods. For this task we need to orthogonalize a matrix against another matrix with orthonormal columns. A common approach is to employ the block Gram--Schmidt algorithm. However, its stability largely depends on the condition number of . While performing a Householder orthogonalization on is unconditionally stable, it does not utilize the knowledge that 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 , our algorithm only needs to orthogonalizes a square submatrix of . Theoretical analysis and numerical experiments demonstrate that our method is also unconditionally stable.
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}
}