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
COMPLEX *16 function zdotc(N, ZX, INCX, ZY, INCY)
ZDOTC
subroutine zher2(UPLO, N, ALPHA, X, INCX, Y, INCY, A, LDA)
ZHER2
subroutine zlarfy(UPLO, N, V, INCV, TAU, C, LDC, WORK)
ZLARFY