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 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
subroutine zaxpy(N, ZA, ZX, INCX, ZY, INCY)
ZAXPY