109 SUBROUTINE clarfy( UPLO, N, V, INCV, TAU, C, LDC, WORK )
122 COMPLEX c( ldc, * ), v( * ), work( * )
128 COMPLEX one, zero, half
129 parameter ( one = ( 1.0e+0, 0.0e+0 ),
130 $ zero = ( 0.0e+0, 0.0e+0 ),
131 $ half = ( 0.5e+0, 0.0e+0 ) )
150 CALL chemv( uplo, n, one, c, ldc, v, incv, zero, work, 1 )
152 alpha = -half*tau*
cdotc( n, work, 1, v, incv )
153 CALL caxpy( n, alpha, v, incv, work, 1 )
157 CALL cher2( uplo, n, -tau, v, incv, work, 1, c, ldc )
complex function cdotc(N, CX, INCX, CY, INCY)
CDOTC
subroutine chemv(UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
CHEMV
subroutine cher2(UPLO, N, ALPHA, X, INCX, Y, INCY, A, LDA)
CHER2
subroutine caxpy(N, CA, CX, INCX, CY, INCY)
CAXPY
subroutine clarfy(UPLO, N, V, INCV, TAU, C, LDC, WORK)
CLARFY