Another orthogonal matrix, revisited (Q1654412)

In [J. Comput. Phys. 8, 292--294 (1971; Zbl 0221.65069)], \textit{A. S. Householder} and \textit{K. Fox} introduced a method for computing an orthonormal basis for the range of a projection. However, this method can lead to numerical inaccuracies in some situations, which can be avoided by a method decribed in [\textit{B. N. Parlett} and \textit{E. Barszcz}, Linear Algebra Appl. 417, No. 2--3, 342--346 (2006; Zbl 1101.15026)]. The main result of this paper is an extension of the method of the paper of Parlett and Barszcz [loc. cit.] to ``tall thin'' matrices. In other words, a solution is given to the problem of completing the unitary \(k\)-Hessenberg matrix given its first \(k\) columns (an \(n \times n\) matrix \(A\) is called unitary \(k\)-Hessenberg if it is unitary and the entries above it's \(k\)-th superdiagonal are all zero).