her

Computes a rank-1 update of a Hermitian matrix.

Syntax

void her(queue &exec_queue, uplo upper_lower, std::int64_t n, T alpha, buffer<T, 1> &x, std::int64_t incx, buffer<T, 1> &a, std::int64_t lda)

her supports the following precisions and devices.

T

Devices Supported

std::complex<float>

Host, CPU, and GPU

std::complex<double>

Host, CPU, and GPU

Description

The her routines compute a scalar-vector-vector product and add the result to a Hermitian matrix. The operation is defined as

A <- alpha*x*xH + A

where:

alpha is scalar,

A is an n-by-n Hermitian matrix,

x is a vector of length n.

Input Parameters

exec_queue

The queue where the routine should be executed.

upper_lower

Specifies whether A is upper or lower triangular. See Data Types for more details.

n

Number of rows and columns of A. Must be at least zero.

alpha

Scaling factor for the matrix-vector product.

x

Buffer holding input vector x. The buffer must be of size at least (1 + (n - 1)*abs(incx)). See Matrix and Vector Storage for more details.

incx

Stride of vector x.

a

Buffer holding input matrix A. Must have size at least lda*n. See Matrix and Vector Storage for more details.

lda

Leading dimension of matrix A. Must be at least n, and positive.

Output Parameters

a

Buffer holding the updated upper triangular part of theHermitian matrix A if upper_lower = upper or the updated lowertriangular part of the Hermitian matrix A if upper_lower = lower.

The imaginary parts of the diagonal elementsare set to zero.