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

◆ zspr2()

void zspr2 ( char  uplo,
int  n,
doublecomplex  alpha,
doublecomplex  x[],
int  incx,
doublecomplex  y[],
int  incy,
doublecomplex  ap[] 
)

Rank 2 operation: A <- αxyT + αyxT + A (complex symmetric matrices in packed form) (BLAS 2)

Purpose
This routine performs the rank 2 operation
A <- alpha*x*y^T + alpha*y*x^T + A
where alpha is a scalar, x and y are n element vectors and A is an n x n symmetric matrix.
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]y[]Array y[ly] (ly >= 1 + (n - 1)*abs(incy))
Vector y.
[in]incyStorage spacing between elements of y. (incy != 0)
[in,out]ap[]Array ap[lap] (lap >= n(n + 1)/2)
[in] n x n symmetric matrix A in packed form. (Upper or lower triangle)
[out] Output matrix in packed form. (= alpha*x*y^T + alpha*y*x^T + A) (Upper or lower triangle)
Reference
BLAS