141 SUBROUTINE clatrz( M, N, L, A, LDA, TAU, WORK )
152 COMPLEX A( lda, * ), TAU( * ), WORK( * )
159 parameter ( zero = ( 0.0e+0, 0.0e+0 ) )
177 ELSE IF( m.EQ.n )
THEN
189 CALL clacgv( l, a( i, n-l+1 ), lda )
190 alpha = conjg( a( i, i ) )
191 CALL clarfg( l+1, alpha, a( i, n-l+1 ), lda, tau( i ) )
192 tau( i ) = conjg( tau( i ) )
196 CALL clarz(
'Right', i-1, n-i+1, l, a( i, n-l+1 ), lda,
197 $ conjg( tau( i ) ), a( 1, i ), lda, work )
198 a( i, i ) = conjg( alpha )
subroutine clarfg(N, ALPHA, X, INCX, TAU)
CLARFG generates an elementary reflector (Householder matrix).
subroutine clatrz(M, N, L, A, LDA, TAU, WORK)
CLATRZ factors an upper trapezoidal matrix by means of unitary transformations.
subroutine clacgv(N, X, INCX)
CLACGV conjugates a complex vector.
subroutine clarz(SIDE, M, N, L, V, INCV, TAU, C, LDC, WORK)
CLARZ applies an elementary reflector (as returned by stzrzf) to a general matrix.