166 SUBROUTINE zgemlqt( SIDE, TRANS, M, N, K, MB, V, LDV, T, LDT,
167 $ C, LDC, WORK, INFO )
174 CHARACTER SIDE, TRANS
175 INTEGER INFO, K, LDV, LDC, M, N, MB, LDT
178 COMPLEX*16 V( LDV, * ), C( LDC, * ), T( LDT, * ), WORK( * )
185 LOGICAL LEFT, RIGHT, TRAN, NOTRAN
186 INTEGER I, IB, LDWORK, KF, Q
203 left = lsame( side,
'L' )
204 right = lsame( side,
'R' )
205 tran = lsame( trans,
'C' )
206 notran = lsame( trans,
'N' )
211 ELSE IF ( right )
THEN
215 IF( .NOT.left .AND. .NOT.right )
THEN
217 ELSE IF( .NOT.tran .AND. .NOT.notran )
THEN
219 ELSE IF( m.LT.0 )
THEN
221 ELSE IF( n.LT.0 )
THEN
223 ELSE IF( k.LT.0 .OR. k.GT.q )
THEN
225 ELSE IF( mb.LT.1 .OR. (mb.GT.k .AND. k.GT.0))
THEN
227 ELSE IF( ldv.LT.max( 1, k ) )
THEN
229 ELSE IF( ldt.LT.mb )
THEN
231 ELSE IF( ldc.LT.max( 1, m ) )
THEN
236 CALL xerbla(
'ZGEMLQT', -info )
242 IF( m.EQ.0 .OR. n.EQ.0 .OR. k.EQ.0 )
RETURN
244 IF( left .AND. notran )
THEN
247 ib = min( mb, k-i+1 )
248 CALL zlarfb(
'L',
'C',
'F',
'R', m-i+1, n, ib,
249 $ v( i, i ), ldv, t( 1, i ), ldt,
250 $ c( i, 1 ), ldc, work, ldwork )
253 ELSE IF( right .AND. tran )
THEN
256 ib = min( mb, k-i+1 )
257 CALL zlarfb(
'R',
'N',
'F',
'R', m, n-i+1, ib,
258 $ v( i, i ), ldv, t( 1, i ), ldt,
259 $ c( 1, i ), ldc, work, ldwork )
262 ELSE IF( left .AND. tran )
THEN
266 ib = min( mb, k-i+1 )
267 CALL zlarfb(
'L',
'N',
'F',
'R', m-i+1, n, ib,
268 $ v( i, i ), ldv, t( 1, i ), ldt,
269 $ c( i, 1 ), ldc, work, ldwork )
272 ELSE IF( right .AND. notran )
THEN
276 ib = min( mb, k-i+1 )
277 CALL zlarfb(
'R',
'C',
'F',
'R', m, n-i+1, ib,
278 $ v( i, i ), ldv, t( 1, i ), ldt,
279 $ c( 1, i ), ldc, work, ldwork )
subroutine xerbla(SRNAME, INFO)
XERBLA
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 zgemlqt(SIDE, TRANS, M, N, K, MB, V, LDV, T, LDT, C, LDC, WORK, INFO)
ZGEMLQT