194 INTEGER INFO, LDA, LDT, LWORK, M, N, MB, NB
197 REAL A( LDA, * ), T( LDT, * ), WORK( * )
204 parameter( one = 1.0e+0, zero = 0.0e+0 )
208 INTEGER NBLOCAL, MB2, M_PLUS_ONE, ITMP, IB_BOTTOM,
209 $ lworkopt, num_all_row_blocks, jb_t, ib, imb,
210 $ kb, kb_last, knb, mb1
217 EXTERNAL sroundup_lwork
233 ELSE IF( n.LT.0 .OR. m.LT.n )
THEN
235 ELSE IF( mb.LE.n )
THEN
237 ELSE IF( nb.LT.1 )
THEN
239 ELSE IF( lda.LT.max( 1, m ) )
THEN
241 ELSE IF( ldt.LT.max( 1, min( nb, n ) ) )
THEN
243 ELSE IF( lwork.LT.1 .AND. .NOT.lquery )
THEN
247 nblocal = min( nb, n )
252 lworkopt = nblocal * max( nblocal, ( n - nblocal ) )
258 CALL xerbla(
'SORGTSQR_ROW', -info )
260 ELSE IF ( lquery )
THEN
261 work( 1 ) = sroundup_lwork( lworkopt )
267 IF( min( m, n ).EQ.0 )
THEN
268 work( 1 ) = sroundup_lwork( lworkopt )
275 CALL slaset(
'U', m, n, zero, one, a, lda )
280 kb_last = ( ( n-1 ) / nblocal ) * nblocal + 1
302 itmp = ( m - mb - 1 ) / mb2
303 ib_bottom = itmp * mb2 + mb + 1
304 num_all_row_blocks = itmp + 2
305 jb_t = num_all_row_blocks * n + 1
307 DO ib = ib_bottom, mb+1, -mb2
312 imb = min( m_plus_one - ib, mb2 )
324 DO kb = kb_last, 1, -nblocal
329 knb = min( nblocal, n - kb + 1 )
332 $ t( 1, jb_t+kb-1 ), ldt, a( kb, kb ), lda,
333 $ a( ib, kb ), lda, work, knb )
352 DO kb = kb_last, 1, -nblocal
357 knb = min( nblocal, n - kb + 1 )
359 IF( mb1-kb-knb+1.EQ.0 )
THEN
366 $ t( 1, kb ), ldt, a( kb, kb ), lda,
367 $ dummy( 1, 1 ), 1, work, knb )
370 $ t( 1, kb ), ldt, a( kb, kb ), lda,
371 $ a( kb+knb, kb), lda, work, knb )
377 work( 1 ) = sroundup_lwork( lworkopt )