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