189 SUBROUTINE dtplqt( M, N, L, MB, A, LDA, B, LDB, T, LDT, WORK,
198 INTEGER INFO, LDA, LDB, LDT, N, M, L, MB
201 DOUBLE PRECISION A( lda, * ), B( ldb, * ), T( ldt, * ), WORK( * )
208 INTEGER I, IB, LB, NB, IINFO
220 ELSE IF( n.LT.0 )
THEN
222 ELSE IF( l.LT.0 .OR. (l.GT.min(m,n) .AND. min(m,n).GE.0))
THEN
224 ELSE IF( mb.LT.1 .OR. (mb.GT.m .AND. m.GT.0))
THEN
226 ELSE IF( lda.LT.max( 1, m ) )
THEN
228 ELSE IF( ldb.LT.max( 1, m ) )
THEN
230 ELSE IF( ldt.LT.mb )
THEN
234 CALL xerbla(
'DTPLQT', -info )
240 IF( m.EQ.0 .OR. n.EQ.0 )
RETURN
246 ib = min( m-i+1, mb )
247 nb = min( n-l+i+ib-1, n )
254 CALL dtplqt2( ib, nb, lb, a(i,i), lda, b( i, 1 ), ldb,
255 $ t(1, i ), ldt, iinfo )
260 CALL dtprfb(
'R',
'N',
'F',
'R', m-i-ib+1, nb, ib, lb,
261 $ b( i, 1 ), ldb, t( 1, i ), ldt,
262 $ a( i+ib, i ), lda, b( i+ib, 1 ), ldb,
subroutine dtprfb(SIDE, TRANS, DIRECT, STOREV, M, N, K, L, V, LDV, T, LDT, A, LDA, B, LDB, WORK, LDWORK)
DTPRFB applies a real or complex "triangular-pentagonal" blocked reflector to a real or complex matri...
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine dtplqt2(M, N, L, A, LDA, B, LDB, T, LDT, INFO)
DTPLQT2 computes a LQ factorization of a real or complex "triangular-pentagonal" matrix, which is composed of a triangular block and a pentagonal block, using the compact WY representation for Q.
subroutine dtplqt(M, N, L, MB, A, LDA, B, LDB, T, LDT, WORK, INFO)
DTPLQT