.. _onemkl_blas_spr:

spr
===

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;

.. math::

      A \leftarrow alpha*x*x^T + A

where:

- ``alpha`` is scalar

- ``A`` is ``n`` x ``n`` symmetric matrix, supplied in packed form

- ``x`` is a vector of length ``n``

``spr`` supports the following precisions:

.. list-table::
   :header-rows: 1

   * -  T
   * -  ``float``
   * -  ``double``


spr (Buffer Version)
********************

Syntax
------

.. code-block:: cpp

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

.. code-block:: cpp

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

Input Parameters
----------------

queue
   The queue where the routine should be executed.

upper_lower
   Specifies whether matrix ``A`` is upper or lower triangular. See :ref:`data-types` for more details.

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

alpha
   Scaling factor for the matrix-vector product.

x
   Buffer holding input vector ``x``. Size of the buffer must be at
   least (1 + (``n`` - 1)\*abs(``incx``)). See :ref:`matrix-storage` for more details.

incx
   Stride of vector ``x``.

a
   Buffer holding input matrix ``A``. Size of the buffer must be at least
   (``n``\ \*(``n``-n))/2. See :ref:`matrix-storage` for more details.

Output Parameters
-----------------

a
   Buffer holding updated upper triangular part of the symmetric
   matrix ``A`` if ``upper_lower=upper``, or updated lower
   triangular part of the symmetric matrix ``A`` if
   ``upper_lower=lower``.


spr (USM Version)
*****************

Syntax
------

.. code-block:: cpp

   namespace oneapi::mkl::blas::column_major {
       sycl::event spr(sycl::queue &queue,
                       oneapi::mkl::uplo upper_lower,
                       std::int64_t n,
                       T alpha,
                       const T *x,
                       std::int64_t incx,
                       T *a,
                       const std::vector<sycl::event> &dependencies = {})
   }

.. code-block:: cpp

   namespace oneapi::mkl::blas::row_major {
       sycl::event spr(sycl::queue &queue,
                       oneapi::mkl::uplo upper_lower,
                       std::int64_t n,
                       T alpha,
                       const T *x,
                       std::int64_t incx,
                       T *a,
                       const std::vector<sycl::event> &dependencies = {})
   }

Input Parameters
----------------

queue
   The queue where the routine should be executed.

upper_lower
   Specifies whether matrix ``A`` is upper or lower triangular. See :ref:`data-types` for more details.

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

alpha
   Scaling factor for the matrix-vector product.

x
   Pointer to input vector ``x``. Size of the array holding input vector
   ``x`` must be at least (1 + (``n`` - 1)\*abs(``incx``)).
   See :ref:`matrix-storage` for more details.

incx
   Stride of vector ``x``.

a
   Pointer to input matrix ``A``. Size of the array holding input matrix
   ``A`` must be at least (``n``\ \*(``n``-n))/2. See :ref:`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 updated upper triangular part of the symmetric
   matrix ``A`` if ``upper_lower=upper``, or 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.