zhpr2 - perform the Hermitian rank 2 operation A := alpha*x*conjg( y' ) + conjg( alpha )*y*conjg( x' ) + A
SUBROUTINE ZHPR2( UPLO, N, ALPHA, X, INCX, Y, INCY, A) CHARACTER * 1 UPLO DOUBLE COMPLEX ALPHA DOUBLE COMPLEX X(*), Y(*), A(*) INTEGER N, INCX, INCY
SUBROUTINE ZHPR2_64( UPLO, N, ALPHA, X, INCX, Y, INCY, A) CHARACTER * 1 UPLO DOUBLE COMPLEX ALPHA DOUBLE COMPLEX X(*), Y(*), A(*) INTEGER*8 N, INCX, INCY
SUBROUTINE HPR2( UPLO, [N], ALPHA, X, [INCX], Y, [INCY], A) CHARACTER(LEN=1) :: UPLO COMPLEX(8) :: ALPHA COMPLEX(8), DIMENSION(:) :: X, Y, A INTEGER :: N, INCX, INCY
SUBROUTINE HPR2_64( UPLO, [N], ALPHA, X, [INCX], Y, [INCY], A) CHARACTER(LEN=1) :: UPLO COMPLEX(8) :: ALPHA COMPLEX(8), DIMENSION(:) :: X, Y, A INTEGER(8) :: N, INCX, INCY
#include <sunperf.h>
void zhpr2(char uplo, int n, doublecomplex alpha, doublecomplex *x, int incx, doublecomplex *y, int incy, doublecomplex *a);
void zhpr2_64(char uplo, long n, doublecomplex alpha, doublecomplex *x, long incx, doublecomplex *y, long incy, doublecomplex *a);
zhpr2 performs the Hermitian rank 2 operation A := alpha*x*conjg( y' ) + conjg( alpha )*y*conjg( x' ) + A where alpha is a scalar, x and y are n element vectors and A is an n by n hermitian matrix, supplied in packed form.
UPLO = 'U' or 'u' The upper triangular part of A is supplied in A.
UPLO = 'L' or 'l' The lower triangular part of A is supplied in A.
Unchanged on exit.