159 SUBROUTINE zunmr2( 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, 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(
'ZUNMR2', -info )
233 IF( m.EQ.0 .OR. n.EQ.0 .OR. k.EQ.0 )
236 IF( ( left .AND. .NOT.notran .OR. .NOT.left .AND. notran ) )
THEN
268 taui = dconjg( tau( i ) )
272 CALL zlacgv( nq-k+i-1, a( i, 1 ), lda )
275 CALL zlarf( side, mi, ni, a( i, 1 ), lda, taui, c, ldc, work )
277 CALL zlacgv( nq-k+i-1, a( i, 1 ), lda )
subroutine zunmr2(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, INFO)
ZUNMR2 multiplies a general matrix by the unitary matrix from a RQ factorization determined by cgerqf...
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 zlacgv(N, X, INCX)
ZLACGV conjugates a complex vector.