141 SUBROUTINE slatrz( M, N, L, A, LDA, TAU, WORK )
152 REAL A( lda, * ), TAU( * ), WORK( * )
159 parameter ( zero = 0.0e+0 )
175 ELSE IF( m.EQ.n )
THEN
187 CALL slarfg( l+1, a( i, i ), a( i, n-l+1 ), lda, tau( i ) )
191 CALL slarz(
'Right', i-1, n-i+1, l, a( i, n-l+1 ), lda,
192 $ tau( i ), a( 1, i ), lda, work )
subroutine slatrz(M, N, L, A, LDA, TAU, WORK)
SLATRZ factors an upper trapezoidal matrix by means of orthogonal transformations.
subroutine slarz(SIDE, M, N, L, V, INCV, TAU, C, LDC, WORK)
SLARZ applies an elementary reflector (as returned by stzrzf) to a general matrix.
subroutine slarfg(N, ALPHA, X, INCX, TAU)
SLARFG generates an elementary reflector (Householder matrix).