1 SUBROUTINE dorg2l( M, N, K, A, LDA, TAU, WORK, INFO )
9 INTEGER INFO, K, LDA, M, N
12 DOUBLE PRECISION A( LDA, * ), TAU( * ), WORK( * )
62 DOUBLE PRECISION ONE, ZERO
63 parameter( one = 1.0d+0, zero = 0.0d+0 )
81 ELSE IF( n.LT.0 .OR. n.GT.m )
THEN
83 ELSE IF( k.LT.0 .OR. k.GT.n )
THEN
85 ELSE IF( lda.LT.max( 1, m ) )
THEN
89 CALL xerbla(
'DORG2L', -info )
112 a( m-n+ii, ii ) = one
113 CALL dlarf(
'Left', m-n+ii, ii-1, a( 1, ii ), 1, tau( i ), a,
115 CALL dscal( m-n+ii-1, -tau( i ), a( 1, ii ), 1 )
116 a( m-n+ii, ii ) = one - tau( i )
120 DO 30 l = m - n + ii + 1, m
subroutine dlarf(SIDE, M, N, V, INCV, TAU, C, LDC, WORK)
subroutine dorg2l(M, N, K, A, LDA, TAU, WORK, INFO)
subroutine dscal(n, da, dx, incx)
subroutine xerbla(SRNAME, INFO)