spr (USM Version)

Performs a rank-1 update of a symmetric packed matrix.

Description

The spr routines compute a scalar-vector-vector product and add the result to a symmetric packed matrix. The operation is defined as

A <- alpha*x*xT + A

where:

  • alpha is scalar,

  • A is an n-by-n symmetric matrix, supplied in packed form,

  • x is a vector of length n.

API

Syntax

event spr(queue &exec_queue,
  uplo upper_lower,
  std::int64_t n,
  T alpha,
  const T *x,
  std::int64_t incx,
  T *a,
  const vector_class<event> &dependencies = {})

The USM version of spr supports the following precisions and devices.

T

Devices Supported

float

Host, CPU, and GPU

double

Host, CPU, and GPU

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

Pointer to input vector x. The array holding input vector x must be of size at least (1 + (n - 1)*abs(incx)). See Matrix Storage for more details.

incx

Stride of vector x.

a

Pointer to input matrix A. The array holding input matrix A must have size at least (n*(n-n))/2. See Matrix Storage for more details.

dependencies

List of events to wait for before starting computation, if any. If omitted, defaults to no dependencies.

Output Parameters

a

Pointer to the updated upper triangular part of the symmetric matrix A if upper_lower =upper, or the updated lower triangular part of the symmetric matrix A if upper_lower =lower.

Return Values

Output event to wait on to ensure computation is complete.