115 SUBROUTINE zung2r( M, N, K, A, LDA, TAU, WORK, INFO )
123 INTEGER INFO, K, LDA, M, N
126 COMPLEX*16 A( lda, * ), TAU( * ), WORK( * )
133 parameter ( one = ( 1.0d+0, 0.0d+0 ),
134 $ zero = ( 0.0d+0, 0.0d+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(
'ZUNG2R', -info )
184 CALL zlarf(
'Left', m-i+1, n-i, a( i, i ), 1, tau( i ),
185 $ a( i, i+1 ), lda, work )
188 $
CALL zscal( m-i, -tau( i ), a( i+1, i ), 1 )
189 a( i, i ) = one - tau( i )
subroutine zung2r(M, N, K, A, LDA, TAU, WORK, INFO)
ZUNG2R
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine zlarf(SIDE, M, N, V, INCV, TAU, C, LDC, WORK)
ZLARF applies an elementary reflector to a general rectangular matrix.
subroutine zscal(N, ZA, ZX, INCX)
ZSCAL