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
subroutine cher2(UPLO, N, ALPHA, X, INCX, Y, INCY, A, LDA)
CHER2