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 )
double precision function ddot(N, DX, INCX, DY, INCY)
DDOT
subroutine daxpy(N, DA, DX, INCX, DY, INCY)
DAXPY
subroutine dsyr2(UPLO, N, ALPHA, X, INCX, Y, INCY, A, LDA)
DSYR2
subroutine dlarfy(UPLO, N, V, INCV, TAU, C, LDC, WORK)
DLARFY
subroutine dsymv(UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
DSYMV