115 SUBROUTINE cung2r( M, N, K, A, LDA, TAU, WORK, INFO )
123 INTEGER INFO, K, LDA, M, N
126 COMPLEX A( lda, * ), TAU( * ), WORK( * )
133 parameter ( one = ( 1.0e+0, 0.0e+0 ),
134 $ zero = ( 0.0e+0, 0.0e+0 ) )
152 ELSE IF( n.LT.0 .OR. n.GT.m )
THEN
154 ELSE IF( k.LT.0 .OR. k.GT.n )
THEN
156 ELSE IF( lda.LT.max( 1, m ) )
THEN
160 CALL xerbla(
'CUNG2R', -info )
184 CALL clarf(
'Left', m-i+1, n-i, a( i, i ), 1, tau( i ),
185 $ a( i, i+1 ), lda, work )
188 $
CALL cscal( m-i, -tau( i ), a( i+1, i ), 1 )
189 a( i, i ) = one - tau( i )
subroutine cung2r(M, N, K, A, LDA, TAU, WORK, INFO)
CUNG2R
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine cscal(N, CA, CX, INCX)
CSCAL
subroutine clarf(SIDE, M, N, V, INCV, TAU, C, LDC, WORK)
CLARF applies an elementary reflector to a general rectangular matrix.