159 SUBROUTINE zunml2( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC,
168 CHARACTER SIDE, TRANS
169 INTEGER INFO, K, LDA, LDC, M, N
172 COMPLEX*16 A( lda, * ), C( ldc, * ), TAU( * ), WORK( * )
179 parameter ( one = ( 1.0d+0, 0.0d+0 ) )
183 INTEGER I, I1, I2, I3, IC, JC, MI, NI, NQ
194 INTRINSIC dconjg, max
201 left = lsame( side,
'L' )
202 notran = lsame( trans,
'N' )
211 IF( .NOT.left .AND. .NOT.lsame( side,
'R' ) )
THEN
213 ELSE IF( .NOT.notran .AND. .NOT.lsame( trans,
'C' ) )
THEN
215 ELSE IF( m.LT.0 )
THEN
217 ELSE IF( n.LT.0 )
THEN
219 ELSE IF( k.LT.0 .OR. k.GT.nq )
THEN
221 ELSE IF( lda.LT.max( 1, k ) )
THEN
223 ELSE IF( ldc.LT.max( 1, m ) )
THEN
227 CALL xerbla(
'ZUNML2', -info )
233 IF( m.EQ.0 .OR. n.EQ.0 .OR. k.EQ.0 )
236 IF( ( left .AND. notran .OR. .NOT.left .AND. .NOT.notran ) )
THEN
272 taui = dconjg( tau( i ) )
277 $
CALL zlacgv( nq-i, a( i, i+1 ), lda )
280 CALL zlarf( side, mi, ni, a( i, i ), lda, taui, c( ic, jc ),
284 $
CALL zlacgv( nq-i, a( i, i+1 ), lda )
subroutine zlarf(SIDE, M, N, V, INCV, TAU, C, LDC, WORK)
ZLARF applies an elementary reflector to a general rectangular matrix.
subroutine zunml2(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, INFO)
ZUNML2 multiplies a general matrix by the unitary matrix from a LQ factorization determined by cgelqf...
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine zlacgv(N, X, INCX)
ZLACGV conjugates a complex vector.