141 SUBROUTINE dlatrz( M, N, L, A, LDA, TAU, WORK )
152 DOUBLE PRECISION A( lda, * ), TAU( * ), WORK( * )
158 DOUBLE PRECISION ZERO
159 parameter ( zero = 0.0d+0 )
175 ELSE IF( m.EQ.n )
THEN
187 CALL dlarfg( l+1, a( i, i ), a( i, n-l+1 ), lda, tau( i ) )
191 CALL dlarz(
'Right', i-1, n-i+1, l, a( i, n-l+1 ), lda,
192 $ tau( i ), a( 1, i ), lda, work )
subroutine dlarz(SIDE, M, N, L, V, INCV, TAU, C, LDC, WORK)
DLARZ applies an elementary reflector (as returned by stzrzf) to a general matrix.
subroutine dlarfg(N, ALPHA, X, INCX, TAU)
DLARFG generates an elementary reflector (Householder matrix).
subroutine dlatrz(M, N, L, A, LDA, TAU, WORK)
DLATRZ factors an upper trapezoidal matrix by means of orthogonal transformations.