chpr - perform the hermitian rank 1 operation A := alpha*x*conjg( x' ) + A
SUBROUTINE CHPR( UPLO, N, ALPHA, X, INCX, A) CHARACTER * 1 UPLO COMPLEX X(*), A(*) INTEGER N, INCX REAL ALPHA
SUBROUTINE CHPR_64( UPLO, N, ALPHA, X, INCX, A) CHARACTER * 1 UPLO COMPLEX X(*), A(*) INTEGER*8 N, INCX REAL ALPHA
SUBROUTINE HPR( UPLO, [N], ALPHA, X, [INCX], A) CHARACTER(LEN=1) :: UPLO COMPLEX, DIMENSION(:) :: X, A INTEGER :: N, INCX REAL :: ALPHA
SUBROUTINE HPR_64( UPLO, [N], ALPHA, X, [INCX], A) CHARACTER(LEN=1) :: UPLO COMPLEX, DIMENSION(:) :: X, A INTEGER(8) :: N, INCX REAL :: ALPHA
#include <sunperf.h>
void chpr(char uplo, int n, float alpha, complex *x, int incx, complex *a);
void chpr_64(char uplo, long n, float alpha, complex *x, long incx, complex *a);
chpr performs the hermitian rank 1 operation A := alpha*x*conjg( x' ) + A where alpha is a real scalar, x is an n element vector 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.