115 SUBROUTINE sorg2r( M, N, K, A, LDA, TAU, WORK, INFO )
123 INTEGER INFO, K, LDA, M, N
126 REAL A( lda, * ), TAU( * ), WORK( * )
133 parameter ( one = 1.0e+0, zero = 0.0e+0 )
151 ELSE IF( n.LT.0 .OR. n.GT.m )
THEN
153 ELSE IF( k.LT.0 .OR. k.GT.n )
THEN
155 ELSE IF( lda.LT.max( 1, m ) )
THEN
159 CALL xerbla(
'SORG2R', -info )
183 CALL slarf(
'Left', m-i+1, n-i, a( i, i ), 1, tau( i ),
184 $ a( i, i+1 ), lda, work )
187 $
CALL sscal( m-i, -tau( i ), a( i+1, i ), 1 )
188 a( i, i ) = one - tau( i )
subroutine sorg2r(M, N, K, A, LDA, TAU, WORK, INFO)
SORG2R generates all or part of the orthogonal matrix Q from a QR factorization determined by sgeqrf ...
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine sscal(N, SA, SX, INCX)
SSCAL
subroutine slarf(SIDE, M, N, V, INCV, TAU, C, LDC, WORK)
SLARF applies an elementary reflector to a general rectangular matrix.