114 SUBROUTINE sorgl2( M, N, K, A, LDA, TAU, WORK, INFO )
122 INTEGER INFO, K, LDA, M, N
125 REAL A( lda, * ), TAU( * ), WORK( * )
132 parameter ( one = 1.0e+0, zero = 0.0e+0 )
150 ELSE IF( n.LT.m )
THEN
152 ELSE IF( k.LT.0 .OR. k.GT.m )
THEN
154 ELSE IF( lda.LT.max( 1, m ) )
THEN
158 CALL xerbla(
'SORGL2', -info )
175 IF( j.GT.k .AND. j.LE.m )
187 CALL slarf(
'Right', m-i, n-i+1, a( i, i ), lda,
188 $ tau( i ), a( i+1, i ), lda, work )
190 CALL sscal( n-i, -tau( i ), a( i, i+1 ), lda )
192 a( i, i ) = one - tau( i )
subroutine sorgl2(M, N, K, A, LDA, TAU, WORK, INFO)
SORGL2
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.