107 SUBROUTINE slarfy( UPLO, N, V, INCV, TAU, C, LDC, WORK )
119 REAL C( LDC, * ), V( * ), WORK( * )
126 parameter( one = 1.0e+0, zero = 0.0e+0, half = 0.5e+0 )
145 CALL ssymv( uplo, n, one, c, ldc, v, incv, zero, work, 1 )
147 alpha = -half*tau*sdot( n, work, 1, v, incv )
148 CALL saxpy( n, alpha, v, incv, work, 1 )
152 CALL ssyr2( uplo, n, -tau, v, incv, work, 1, c, ldc )
subroutine slarfy(UPLO, N, V, INCV, TAU, C, LDC, WORK)
SLARFY
subroutine saxpy(N, SA, SX, INCX, SY, INCY)
SAXPY
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