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 )