1 SUBROUTINE dorgl2( M, N, K, A, LDA, TAU, WORK, INFO )
9 INTEGER INFO, K, LDA, M, N
12 DOUBLE PRECISION A( LDA, * ), TAU( * ), WORK( * )
61 DOUBLE PRECISION ONE, ZERO
62 parameter( one = 1.0d+0, zero = 0.0d+0 )
80 ELSE IF( n.LT.m )
THEN
82 ELSE IF( k.LT.0 .OR. k.GT.m )
THEN
84 ELSE IF( lda.LT.max( 1, m ) )
THEN
88 CALL xerbla(
'DORGL2', -info )
105 IF( j.GT.k .AND. j.LE.m )
117 CALL dlarf(
'Right', m-i, n-i+1, a( i, i ), lda,
118 $ tau( i ), a( i+1, i ), lda, work )
120 CALL dscal( n-i, -tau( i ), a( i, i+1 ), lda )
122 a( i, i ) = one - tau( i )
subroutine dlarf(SIDE, M, N, V, INCV, TAU, C, LDC, WORK)
subroutine dorgl2(M, N, K, A, LDA, TAU, WORK, INFO)
subroutine dscal(n, da, dx, incx)
subroutine xerbla(SRNAME, INFO)