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 dopgtr(UPLO, N, AP, TAU, Q, LDQ, WORK, INFO)
DOPGTR
subroutine xerbla(SRNAME, INFO)
XERBLA
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 ...