The computation of elementary unitary matrices. The construction of elementary unitary matrices that transform a complex vector to a multiple of e 1 , the first column of the identity matrix, is studied. We present four variants and their software implementation, including a discussion on the LAPACK subroutine CLARFG. Comparisons are also given.

This software is also peer reviewed by journal TOMS.