XLPack 7.0
XLPack Numerical Library (C API) Reference Manual
Loading...
Searching...
No Matches

◆ zhpr()

void zhpr ( char  uplo,
int  n,
double  alpha,
doublecomplex  x[],
int  incx,
doublecomplex  ap[] 
)

Rank 1 operation: A <- αxxH + A (Hermitian matrices in packed form) (BLAS 2)

Purpose
This routine performs the Hermitian rank 1 operation
A <- alpha*x*x^H + A
where alpha is a real scalar, x is an n element vector and A is an n x n Hermitian matrix, supplied in packed form.
Parameters
[in]uploSpecifies whether the upper or lower triangular part of the matrix A is supplied in the packed array ap[] as follows:
= 'U': The upper triangular part of A is supplied in ap[].
= 'L': The lower triangular part of A is supplied in ap[].
[in]nOrder of the matrix A. (n >= 0) (If n = 0, returns without computation)
[in]alphaScalar alpha.
[in]x[]Array x[lx] (lx >= 1 + (n - 1)*abs(incx))
Vector x.
[in]incxStorage spacing between elements of x. (incx != 0)
[in,out]ap[]Array ap[lap] (lap >= n(n + 1)/2)
[in] n x n Hermitian matrix A in packed form. According to uplo, upper or lower triangular part is to be supplied. The imaginary parts of the diagonal elements need not be set and are assumed to be zero.
[out] Output matrix in packed form (= alpha*x*x^H + A). According to uplo, upper or lower triangular part is supplied. The imaginary parts of the diagonal elements are set to zero.
Reference
BLAS