1 SUBROUTINE dlarf( SIDE, M, N, V, INCV, TAU, C, LDC, WORK )
10 INTEGER INCV, LDC, M, N
14 DOUBLE PRECISION C( LDC, * ), V( * ), WORK( * )
69 DOUBLE PRECISION ONE, ZERO
70 parameter( one = 1.0d+0, zero = 0.0d+0 )
81 IF( lsame( side,
'L' ) )
THEN
85 IF( tau.NE.zero )
THEN
89 CALL dgemv(
'Transpose', m, n, one, c, ldc, v, incv, zero,
94 CALL dger( m, n, -tau, v, incv, work, 1, c, ldc )
100 IF( tau.NE.zero )
THEN
104 CALL dgemv(
'No transpose', m, n, one, c, ldc, v, incv,
109 CALL dger( m, n, -tau, work, 1, v, incv, c, ldc )
subroutine dgemv(TRANS, M, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
subroutine dger(M, N, ALPHA, X, INCX, Y, INCY, A, LDA)
subroutine dlarf(SIDE, M, N, V, INCV, TAU, C, LDC, WORK)