115 SUBROUTINE dopgtr( UPLO, N, AP, TAU, Q, LDQ, WORK, INFO )
127 DOUBLE PRECISION AP( * ), Q( ldq, * ), TAU( * ), WORK( * )
133 DOUBLE PRECISION ZERO, ONE
134 parameter ( zero = 0.0d+0, one = 1.0d+0 )
138 INTEGER I, IINFO, IJ, J
155 upper = lsame( uplo,
'U' )
156 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
158 ELSE IF( n.LT.0 )
THEN
160 ELSE IF( ldq.LT.max( 1, n ) )
THEN
164 CALL xerbla(
'DOPGTR', -info )
197 CALL dorg2l( n-1, n-1, n-1, q, ldq, tau, work, iinfo )
224 CALL dorg2r( n-1, n-1, n-1, q( 2, 2 ), ldq, tau, work,
subroutine dorg2l(M, N, K, A, LDA, TAU, WORK, INFO)
DORG2L generates all or part of the orthogonal matrix Q from a QL factorization determined by sgeqlf ...
subroutine dorg2r(M, N, K, A, LDA, TAU, WORK, INFO)
DORG2R generates all or part of the orthogonal matrix Q from a QR factorization determined by sgeqrf ...
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine dopgtr(UPLO, N, AP, TAU, Q, LDQ, WORK, INFO)
DOPGTR