spr

cspr

Functions

void cspr(
    const char*          uplo,
    const INT            n,
    const c64            alpha,
    const c64*  restrict X,
    const INT            incx,
          c64*  restrict AP
);

void cspr(const char *uplo, const INT n, const c64 alpha, const c64 *restrict X, const INT incx, c64 *restrict AP)

CSPR performs the symmetric rank 1 operation.

A := alpha*x*x**T + A,

where alpha is a complex scalar, x is an n element vector and A is an n by n symmetric matrix, supplied in packed form.

Parameters

in

uplo

On entry, UPLO specifies whether the upper or lower triangular part of the matrix A is supplied in the packed array AP as follows:

• = ‘U’ or ‘u’: The upper triangular part of A is supplied in AP.

• = ‘L’ or ‘l’: The lower triangular part of A is supplied in AP.

in

n

The order of the matrix A. n >= 0.

in

alpha

The scalar alpha.

in

X

Complex*16 array, dimension at least (1 + (n-1)*abs(incx)). The n-element vector x.

in

incx

The increment for the elements of X. incx must not be zero.

inout

AP

Complex*16 array, dimension at least (n*(n+1)/2). On entry, the upper or lower triangular part of the symmetric matrix packed sequentially, column by column. On exit, overwritten by the updated matrix.

zspr

Functions

void zspr(
    const char*          uplo,
    const INT            n,
    const c128           alpha,
    const c128* restrict X,
    const INT            incx,
          c128* restrict AP
);

void zspr(const char *uplo, const INT n, const c128 alpha, const c128 *restrict X, const INT incx, c128 *restrict AP)

ZSPR performs the symmetric rank 1 operation.

A := alpha*x*x**T + A,

where alpha is a complex scalar, x is an n element vector and A is an n by n symmetric matrix, supplied in packed form.

Parameters

in

uplo

On entry, UPLO specifies whether the upper or lower triangular part of the matrix A is supplied in the packed array AP as follows:

• = ‘U’ or ‘u’: The upper triangular part of A is supplied in AP.

• = ‘L’ or ‘l’: The lower triangular part of A is supplied in AP.

in

n

The order of the matrix A. n >= 0.

in

alpha

The scalar alpha.

in

X

Complex*16 array, dimension at least (1 + (n-1)*abs(incx)). The n-element vector x.

in

incx

The increment for the elements of X. incx must not be zero.

inout

AP

Complex*16 array, dimension at least (n*(n+1)/2). On entry, the upper or lower triangular part of the symmetric matrix packed sequentially, column by column. On exit, overwritten by the updated matrix.