188 INTEGER INFO, LDA, LDT, LWORK, M, N, NB1, NB2, MB1
191 DOUBLE PRECISION A( LDA, * ), T( LDT, * ), WORK( * )
198 PARAMETER ( ONE = 1.0d+0 )
202 INTEGER I, IINFO, J, LW1, LW2, LWT, LDWT, LWORKOPT,
203 $ nb1local, nb2local, num_all_row_blocks
211 INTRINSIC ceiling, dble, max, min
218 lquery = ( lwork.EQ.-1 )
221 ELSE IF( n.LT.0 .OR. m.LT.n )
THEN
223 ELSE IF( mb1.LE.n )
THEN
225 ELSE IF( nb1.LT.1 )
THEN
227 ELSE IF( nb2.LT.1 )
THEN
229 ELSE IF( lda.LT.max( 1, m ) )
THEN
231 ELSE IF( ldt.LT.max( 1, min( nb2, n ) ) )
THEN
242 IF( lwork.LT.n*n+1 .AND. .NOT.lquery )
THEN
248 nb1local = min( nb1, n )
250 num_all_row_blocks = max( 1,
251 $ ceiling( dble( m - n ) / dble( mb1 - n ) ) )
256 lwt = num_all_row_blocks * n * nb1local
266 lw2 = nb1local * max( nb1local, ( n - nb1local ) )
268 lworkopt = max( lwt + lw1, max( lwt+n*n+lw2, lwt+n*n+n ) )
269 lworkopt = max( 1, lworkopt )
271 IF( lwork.LT.lworkopt .AND. .NOT.lquery )
THEN
281 CALL xerbla(
'DGETSQRHRT', -info )
283 ELSE IF ( lquery )
THEN
284 work( 1 ) = dble( lworkopt )
290 IF( min( m, n ).EQ.0 )
THEN
291 work( 1 ) = dble( lworkopt )
295 nb2local = min( nb2, n )
300 CALL dlatsqr( m, n, mb1, nb1local, a, lda, work, ldwt,
301 $ work(lwt+1), lw1, iinfo )
308 CALL dcopy( j, a( 1, j ), 1, work( lwt + n*(j-1)+1 ), 1 )
315 CALL dorgtsqr_row( m, n, mb1, nb1local, a, lda, work, ldwt,
316 $ work( lwt+n*n+1 ), lw2, iinfo )
321 CALL dorhr_col( m, n, nb2local, a, lda, t, ldt,
322 $ work( lwt+n*n+1 ), iinfo )
339 IF( work( lwt+n*n+i ).EQ.-one )
THEN
341 a( i, j ) = -one * work( lwt+n*(j-1)+i )
344 CALL dcopy( n-i+1, work(lwt+n*(i-1)+i), n, a( i, i ),
349 work( 1 ) = dble( lworkopt )
subroutine dgetsqrhrt(m, n, mb1, nb1, nb2, a, lda, t, ldt, work, lwork, info)
DGETSQRHRT