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 )
subroutine ssymv(UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
SSYMV
subroutine ssyr2(UPLO, N, ALPHA, X, INCX, Y, INCY, A, LDA)
SSYR2
subroutine saxpy(N, SA, SX, INCX, SY, INCY)
SAXPY
REAL function sdot(N, SX, INCX, SY, INCY)
SDOT
subroutine slarfy(UPLO, N, V, INCV, TAU, C, LDC, WORK)
SLARFY