110 SUBROUTINE zungl2( M, N, K, A, LDA, TAU, WORK, INFO )
117 INTEGER INFO, K, LDA, M, N
120 COMPLEX*16 A( LDA, * ), TAU( * ), WORK( * )
127 parameter( one = ( 1.0d+0, 0.0d+0 ),
128 $ zero = ( 0.0d+0, 0.0d+0 ) )
137 INTRINSIC dconjg, max
146 ELSE IF( n.LT.m )
THEN
148 ELSE IF( k.LT.0 .OR. k.GT.m )
THEN
150 ELSE IF( lda.LT.max( 1, m ) )
THEN
154 CALL xerbla(
'ZUNGL2', -info )
171 IF( j.GT.k .AND. j.LE.m )
181 CALL zlacgv( n-i, a( i, i+1 ), lda )
183 CALL zlarf1f(
'Right', m-i, n-i+1, a( i, i ), lda,
184 $ conjg( tau( i ) ), a( i+1, i ), lda,
187 CALL zscal( n-i, -tau( i ), a( i, i+1 ), lda )
188 CALL zlacgv( n-i, a( i, i+1 ), lda )
190 a( i, i ) = one - dconjg( tau( i ) )
subroutine zlarf1f(side, m, n, v, incv, tau, c, ldc, work)
ZLARF1F applies an elementary reflector to a general rectangular
subroutine zungl2(m, n, k, a, lda, tau, work, info)
ZUNGL2 generates all or part of the unitary matrix Q from an LQ factorization determined by cgelqf (u...