109 SUBROUTINE dlarfy( UPLO, N, V, INCV, TAU, C, LDC, WORK )
122 DOUBLE PRECISION c( ldc, * ), v( * ), work( * )
128 DOUBLE PRECISION one, zero, half
129 parameter( one = 1.0d+0, zero = 0.0d+0, half = 0.5d+0 )
132 DOUBLE PRECISION alpha
138 DOUBLE PRECISION ddot
148 CALL
dsymv( uplo, n, one, c, ldc, v, incv, zero, work, 1 )
150 alpha = -half*tau*
ddot( n, work, 1, v, incv )
151 CALL
daxpy( n, alpha, v, incv, work, 1 )
155 CALL
dsyr2( uplo, n, -tau, v, incv, work, 1, c, ldc )