159 SUBROUTINE cunml2( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC,
168 CHARACTER SIDE, TRANS
169 INTEGER INFO, K, LDA, LDC, M, N
172 COMPLEX A( lda, * ), C( ldc, * ), TAU( * ), WORK( * )
179 parameter ( one = ( 1.0e+0, 0.0e+0 ) )
183 INTEGER I, I1, I2, I3, IC, JC, MI, NI, NQ
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(
'CUNML2', -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 = conjg( tau( i ) )
277 $
CALL clacgv( nq-i, a( i, i+1 ), lda )
280 CALL clarf( side, mi, ni, a( i, i ), lda, taui, c( ic, jc ),
284 $
CALL clacgv( nq-i, a( i, i+1 ), lda )
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine clarf(SIDE, M, N, V, INCV, TAU, C, LDC, WORK)
CLARF applies an elementary reflector to a general rectangular matrix.
subroutine cunml2(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, INFO)
CUNML2 multiplies a general matrix by the unitary matrix from a LQ factorization determined by cgelqf...
subroutine clacgv(N, X, INCX)
CLACGV conjugates a complex vector.