185 SUBROUTINE stplqt( M, N, L, MB, A, LDA, B, LDB, T, LDT, WORK,
193 INTEGER INFO, LDA, LDB, LDT, N, M, L, MB
196 REAL A( LDA, * ), B( LDB, * ), T( LDT, * ), WORK( * )
203 INTEGER I, IB, LB, NB, IINFO
215 ELSE IF( n.LT.0 )
THEN
217 ELSE IF( l.LT.0 .OR. (l.GT.min(m,n) .AND. min(m,n).GE.0))
THEN
219 ELSE IF( mb.LT.1 .OR. (mb.GT.m .AND. m.GT.0))
THEN
221 ELSE IF( lda.LT.max( 1, m ) )
THEN
223 ELSE IF( ldb.LT.max( 1, m ) )
THEN
225 ELSE IF( ldt.LT.mb )
THEN
229 CALL xerbla(
'STPLQT', -info )
235 IF( m.EQ.0 .OR. n.EQ.0 )
RETURN
241 ib = min( m-i+1, mb )
242 nb = min( n-l+i+ib-1, n )
249 CALL stplqt2( ib, nb, lb, a(i,i), lda, b( i, 1 ), ldb,
250 $ t(1, i ), ldt, iinfo )
255 CALL stprfb(
'R',
'N',
'F',
'R', m-i-ib+1, nb, ib, lb,
256 $ b( i, 1 ), ldb, t( 1, i ), ldt,
257 $ a( i+ib, i ), lda, b( i+ib, 1 ), ldb,
subroutine stplqt2(m, n, l, a, lda, b, ldb, t, ldt, info)
STPLQT2 computes a LQ factorization of a real or complex "triangular-pentagonal" matrix,...
subroutine stplqt(m, n, l, mb, a, lda, b, ldb, t, ldt, work, info)
STPLQT
subroutine stprfb(side, trans, direct, storev, m, n, k, l, v, ldv, t, ldt, a, lda, b, ldb, work, ldwork)
STPRFB applies a real "triangular-pentagonal" block reflector to a real matrix, which is composed of ...