109 SUBROUTINE slarfy( UPLO, N, V, INCV, TAU, C, LDC, WORK )
122 REAL c( ldc, * ), v( * ), work( * )
129 parameter( one = 1.0e+0, zero = 0.0e+0, half = 0.5e+0 )
148 CALL
ssymv( uplo, n, one, c, ldc, v, incv, zero, work, 1 )
150 alpha = -half*tau*
sdot( n, work, 1, v, incv )
151 CALL
saxpy( n, alpha, v, incv, work, 1 )
155 CALL
ssyr2( uplo, n, -tau, v, incv, work, 1, c, ldc )