187 SUBROUTINE ctpqrt( M, N, L, NB, A, LDA, B, LDB, T, LDT, WORK,
195 INTEGER INFO, LDA, LDB, LDT, N, M, L, NB
198 COMPLEX A( LDA, * ), B( LDB, * ), T( LDT, * ), WORK( * )
205 INTEGER I, IB, LB, MB, IINFO
217 ELSE IF( n.LT.0 )
THEN
219 ELSE IF( l.LT.0 .OR. (l.GT.min(m,n) .AND. min(m,n).GE.0))
THEN
221 ELSE IF( nb.LT.1 .OR. (nb.GT.n .AND. n.GT.0))
THEN
223 ELSE IF( lda.LT.max( 1, n ) )
THEN
225 ELSE IF( ldb.LT.max( 1, m ) )
THEN
227 ELSE IF( ldt.LT.nb )
THEN
231 CALL xerbla(
'CTPQRT', -info )
237 IF( m.EQ.0 .OR. n.EQ.0 )
RETURN
243 ib = min( n-i+1, nb )
244 mb = min( m-l+i+ib-1, m )
251 CALL ctpqrt2( mb, ib, lb, a(i,i), lda, b( 1, i ), ldb,
252 $ t(1, i ), ldt, iinfo )
257 CALL ctprfb(
'L',
'C',
'F',
'C', mb, n-i-ib+1, ib, lb,
258 $ b( 1, i ), ldb, t( 1, i ), ldt,
259 $ a( i, i+ib ), lda, b( 1, i+ib ), ldb,
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine ctprfb(SIDE, TRANS, DIRECT, STOREV, M, N, K, L, V, LDV, T, LDT, A, LDA, B, LDB, WORK, LDWORK)
CTPRFB applies a complex "triangular-pentagonal" block reflector to a complex matrix,...
subroutine ctpqrt(M, N, L, NB, A, LDA, B, LDB, T, LDT, WORK, INFO)
CTPQRT
subroutine ctpqrt2(M, N, L, A, LDA, B, LDB, T, LDT, INFO)
CTPQRT2 computes a QR factorization of a real or complex "triangular-pentagonal" matrix,...