159 SUBROUTINE zunm2l( 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, nq ) )
THEN
223 ELSE IF( ldc.LT.max( 1, m ) )
THEN
227 CALL xerbla(
'ZUNM2L', -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
270 taui = dconjg( tau( i ) )
274 CALL zlarf( side, mi, ni, a( 1, i ), 1, taui, c, ldc, work )
subroutine zunm2l(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, INFO)
ZUNM2L multiplies a general matrix by the unitary matrix from a QL factorization determined by cgeqlf...
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.