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 )
subroutine clarfy(UPLO, N, V, INCV, TAU, C, LDC, WORK)
CLARFY
subroutine caxpy(N, CA, CX, INCX, CY, INCY)
CAXPY
subroutine chemv(UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
CHEMV
COMPLEX function cdotc(N, CX, INCX, CY, INCY)
CDOTC
subroutine cher2(UPLO, N, ALPHA, X, INCX, Y, INCY, A, LDA)
CHER2