169 SUBROUTINE zunmql( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC,
170 $ work, lwork, info )
178 CHARACTER SIDE, TRANS
179 INTEGER INFO, K, LDA, LDC, LWORK, M, N
182 COMPLEX*16 A( lda, * ), C( ldc, * ), TAU( * ), WORK( * )
189 parameter( nbmax = 64, ldt = nbmax+1 )
192 LOGICAL LEFT, LQUERY, NOTRAN
193 INTEGER I, I1, I2, I3, IB, IINFO, IWS, LDWORK, LWKOPT,
194 $ mi, nb, nbmin, ni, nq, nw
197 COMPLEX*16 T( ldt, nbmax )
202 EXTERNAL lsame, ilaenv
215 left = lsame( side,
'L' )
216 notran = lsame( trans,
'N' )
217 lquery = ( lwork.EQ.-1 )
228 IF( .NOT.left .AND. .NOT.lsame( side,
'R' ) )
THEN
230 ELSE IF( .NOT.notran .AND. .NOT.lsame( trans,
'C' ) )
THEN
232 ELSE IF( m.LT.0 )
THEN
234 ELSE IF( n.LT.0 )
THEN
236 ELSE IF( k.LT.0 .OR. k.GT.nq )
THEN
238 ELSE IF( lda.LT.max( 1, nq ) )
THEN
240 ELSE IF( ldc.LT.max( 1, m ) )
THEN
245 IF( m.EQ.0 .OR. n.EQ.0 )
THEN
252 nb = min( nbmax, ilaenv( 1,
'ZUNMQL', side // trans, m, n,
258 IF( lwork.LT.nw .AND. .NOT.lquery )
THEN
264 CALL
xerbla(
'ZUNMQL', -info )
266 ELSE IF( lquery )
THEN
272 IF( m.EQ.0 .OR. n.EQ.0 )
THEN
278 IF( nb.GT.1 .AND. nb.LT.k )
THEN
280 IF( lwork.LT.iws )
THEN
282 nbmin = max( 2, ilaenv( 2,
'ZUNMQL', side // trans, m, n, k,
289 IF( nb.LT.nbmin .OR. nb.GE.k )
THEN
293 CALL
zunm2l( side, trans, m, n, k, a, lda, tau, c, ldc, work,
299 IF( ( left .AND. notran ) .OR.
300 $ ( .NOT.left .AND. .NOT.notran ) )
THEN
305 i1 = ( ( k-1 ) / nb )*nb + 1
317 ib = min( nb, k-i+1 )
322 CALL
zlarft(
'Backward',
'Columnwise', nq-k+i+ib-1, ib,
323 $ a( 1, i ), lda, tau( i ), t, ldt )
328 mi = m - k + i + ib - 1
333 ni = n - k + i + ib - 1
338 CALL
zlarfb( side, trans,
'Backward',
'Columnwise', mi, ni,
339 $ ib, a( 1, i ), lda, t, ldt, c, ldc, work,
subroutine zunmql(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK, LWORK, INFO)
ZUNMQL
subroutine zlarfb(SIDE, TRANS, DIRECT, STOREV, M, N, K, V, LDV, T, LDT, C, LDC, WORK, LDWORK)
ZLARFB applies a block reflector or its conjugate-transpose to a general rectangular matrix...
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine zlarft(DIRECT, STOREV, N, K, V, LDV, TAU, T, LDT)
ZLARFT forms the triangular factor T of a block reflector H = I - vtvH
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...