172 SUBROUTINE ctplqt( M, N, L, MB, A, LDA, B, LDB, T, LDT, WORK,
181 INTEGER INFO, LDA, LDB, LDT, N, M, L, MB
184 COMPLEX A( lda, * ), B( ldb, * ), T( ldt, * ), WORK( * )
191 INTEGER I, IB, LB, NB, IINFO
203 ELSE IF( n.LT.0 )
THEN
205 ELSE IF( l.LT.0 .OR. (l.GT.min(m,n) .AND. min(m,n).GE.0))
THEN
207 ELSE IF( mb.LT.1 .OR. (mb.GT.m .AND. m.GT.0))
THEN
209 ELSE IF( lda.LT.max( 1, m ) )
THEN
211 ELSE IF( ldb.LT.max( 1, m ) )
THEN
213 ELSE IF( ldt.LT.mb )
THEN
217 CALL xerbla(
'CTPLQT', -info )
223 IF( m.EQ.0 .OR. n.EQ.0 )
RETURN
229 ib = min( m-i+1, mb )
230 nb = min( n-l+i+ib-1, n )
237 CALL ctplqt2( ib, nb, lb, a(i,i), lda, b( i, 1 ), ldb,
238 $ t(1, i ), ldt, iinfo )
243 CALL ctprfb(
'R',
'N',
'F',
'R', m-i-ib+1, nb, ib, lb,
244 $ b( i, 1 ), ldb, t( 1, i ), ldt,
245 $ a( i+ib, i ), lda, b( i+ib, 1 ), ldb,
subroutine ctplqt(M, N, L, MB, A, LDA, B, LDB, T, LDT, WORK, INFO)
subroutine ctprfb(SIDE, TRANS, DIRECT, STOREV, M, N, K, L, V, LDV, T, LDT, A, LDA, B, LDB, WORK, LDWORK)
CTPRFB applies a real or complex "triangular-pentagonal" blocked reflector to a real or complex matri...
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine ctplqt2(M, N, L, A, LDA, B, LDB, T, LDT, INFO)