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 cscal(N, CA, CX, INCX)
CSCAL
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine cung2r(M, N, K, A, LDA, TAU, WORK, INFO)
CUNG2R
subroutine clarf(SIDE, M, N, V, INCV, TAU, C, LDC, WORK)
CLARF applies an elementary reflector to a general rectangular matrix.