172 SUBROUTINE dorgtsqr( M, N, MB, NB, A, LDA, T, LDT, WORK, LWORK,
181 INTEGER INFO, LDA, LDT, LWORK, M, N, MB, NB
184 DOUBLE PRECISION A( LDA, * ), T( LDT, * ), WORK( * )
190 DOUBLE PRECISION ONE, ZERO
191 parameter( one = 1.0d+0, zero = 0.0d+0 )
195 INTEGER IINFO, LDC, LWORKOPT, LC, LW, NBLOCAL, J
201 INTRINSIC dble, max, min
211 ELSE IF( n.LT.0 .OR. m.LT.n )
THEN
213 ELSE IF( mb.LE.n )
THEN
215 ELSE IF( nb.LT.1 )
THEN
217 ELSE IF( lda.LT.max( 1, m ) )
THEN
219 ELSE IF( ldt.LT.max( 1, min( nb, n ) ) )
THEN
227 IF( lwork.LT.2 .AND. (.NOT.lquery) )
THEN
233 nblocal = min( nb, n )
245 IF( ( lwork.LT.max( 1, lworkopt ) ).AND.(.NOT.lquery) )
THEN
255 CALL xerbla(
'DORGTSQR', -info )
257 ELSE IF ( lquery )
THEN
258 work( 1 ) = dble( lworkopt )
264 IF( min( m, n ).EQ.0 )
THEN
265 work( 1 ) = dble( lworkopt )
280 CALL dlaset(
'F', m, n, zero, one, work, ldc )
287 CALL dlamtsqr(
'L',
'N', m, n, n, mb, nblocal, a, lda, t, ldt,
288 $ work, ldc, work( lc+1 ), lw, iinfo )
295 CALL dcopy( m, work( (j-1)*ldc + 1 ), 1, a( 1, j ), 1 )
298 work( 1 ) = dble( lworkopt )
subroutine dlamtsqr(side, trans, m, n, k, mb, nb, a, lda, t, ldt, c, ldc, work, lwork, info)
DLAMTSQR
subroutine dlaset(uplo, m, n, alpha, beta, a, lda)
DLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values.
subroutine dorgtsqr(m, n, mb, nb, a, lda, t, ldt, work, lwork, info)
DORGTSQR