her2

Computes a rank-2 update of a Hermitian matrix.

Description

The her2 routines compute two scalar-vector-vector products and adds them to a Hermitian matrix. The operation is defined as

\[A \leftarrow alpha*x*yH + conjg(alpha)*y*xH + A\]

where:

  • alpha is a scalar,

  • A is an n-by-n Hermitian matrix.

  • x and y are vectors or length n.

API

Syntax

namespace oneapi::mkl::blas::column_major {
    void her2(sycl::queue &queue,
              onemkl::uplo upper_lower,
              std::int64_t n,
              T alpha,
              sycl::buffer<T,1> &x,
              std::int64_t incx,
              sycl::buffer<T,1> &y,
              std::int64_t incy,
              sycl::buffer<T,1> &a,
              std::int64_t lda)
}
namespace oneapi::mkl::blas::row_major {
    void her2(sycl::queue &queue,
              onemkl::uplo upper_lower,
              std::int64_t n,
              T alpha,
              sycl::buffer<T,1> &x,
              std::int64_t incx,
              sycl::buffer<T,1> &y,
              std::int64_t incy,
              sycl::buffer<T,1> &a,
              std::int64_t lda)
}

her2 supports the following precisions and devices.

T

Devices Supported

std::complex<float>

Host, CPU, and GPU

std::complex<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 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 Storage for more details.

incx

Stride of vector x.

y

Buffer holding input/output vector y. The buffer must be of size at least (1 + (n - 1)*abs(incy)). See Matrix Storage for more details.

incy

Stride of vector y.

a

Buffer holding input matrix A. Must have size at least lda*n. See Matrix 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 elements are set to zero.