BLAS++ 2024.05.31
BLAS C++ API
Loading...
Searching...
No Matches
ger: General matrix rank 1 update

\(A = \alpha xy^H + A\) More...

Functions

template<typename TA , typename TX , typename TY >
void blas::ger (blas::Layout layout, int64_t m, int64_t n, blas::scalar_type< TA, TX, TY > alpha, TX const *x, int64_t incx, TY const *y, int64_t incy, TA *A, int64_t lda)
 General matrix rank-1 update:
 
void blas::ger (blas::Layout layout, int64_t m, int64_t n, float alpha, float const *x, int64_t incx, float const *y, int64_t incy, float *A, int64_t lda)
 CPU, float version.
 
void blas::ger (blas::Layout layout, int64_t m, int64_t n, double alpha, double const *x, int64_t incx, double const *y, int64_t incy, double *A, int64_t lda)
 CPU, double version.
 
void blas::ger (blas::Layout layout, int64_t m, int64_t n, std::complex< float > alpha, std::complex< float > const *x, int64_t incx, std::complex< float > const *y, int64_t incy, std::complex< float > *A, int64_t lda)
 CPU, complex<float> version.
 
void blas::ger (blas::Layout layout, int64_t m, int64_t n, std::complex< double > alpha, std::complex< double > const *x, int64_t incx, std::complex< double > const *y, int64_t incy, std::complex< double > *A, int64_t lda)
 CPU, complex<double> version.
 

Detailed Description

\(A = \alpha xy^H + A\)

Function Documentation

◆ ger()

template<typename TA , typename TX , typename TY >
void blas::ger ( blas::Layout  layout,
int64_t  m,
int64_t  n,
blas::scalar_type< TA, TX, TY >  alpha,
TX const *  x,
int64_t  incx,
TY const *  y,
int64_t  incy,
TA *  A,
int64_t  lda 
)

General matrix rank-1 update:

\[ A = \alpha x y^H + A, \]

where alpha is a scalar, x and y are vectors, and A is an m-by-n matrix.

Generic implementation for arbitrary data types.

Parameters
[in]layoutMatrix storage, Layout::ColMajor or Layout::RowMajor.
[in]mNumber of rows of the matrix A. m >= 0.
[in]nNumber of columns of the matrix A. n >= 0.
[in]alphaScalar alpha. If alpha is zero, A is not updated.
[in]xThe m-element vector x, in an array of length (m-1)*abs(incx) + 1.
[in]incxStride between elements of x. incx must not be zero. If incx < 0, uses elements of x in reverse order: x(n-1), ..., x(0).
[in]yThe n-element vector y, in an array of length (n-1)*abs(incy) + 1.
[in]incyStride between elements of y. incy must not be zero. If incy < 0, uses elements of y in reverse order: y(n-1), ..., y(0).
[in,out]AThe m-by-n matrix A, stored in an lda-by-n array [RowMajor: m-by-lda].
[in]ldaLeading dimension of A. lda >= max(1, m) [RowMajor: lda >= max(1, n)].