109 SUBROUTINE zlarfy( UPLO, N, V, INCV, TAU, C, LDC, WORK )
122 COMPLEX*16 C( ldc, * ), V( * ), WORK( * )
128 COMPLEX*16 ONE, ZERO, HALF
129 parameter( one = ( 1.0d+0, 0.0d+0 ),
130 $ zero = ( 0.0d+0, 0.0d+0 ),
131 $ half = ( 0.5d+0, 0.0d+0 ) )
150 CALL
zhemv( uplo, n, one, c, ldc, v, incv, zero, work, 1 )
152 alpha = -half*tau*zdotc( n, work, 1, v, incv )
153 CALL
zaxpy( n, alpha, v, incv, work, 1 )
157 CALL
zher2( uplo, n, -tau, v, incv, work, 1, c, ldc )
subroutine zaxpy(N, ZA, ZX, INCX, ZY, INCY)
ZAXPY
subroutine zhemv(UPLO, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
ZHEMV
subroutine zher2(UPLO, N, ALPHA, X, INCX, Y, INCY, A, LDA)
ZHER2
subroutine zlarfy(UPLO, N, V, INCV, TAU, C, LDC, WORK)
ZLARFY