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 )