geru: General matrix rank 1 update, unconjugated

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

Functions
template<typename TA , typename TX , typename TY >
void	blas::geru (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::geru (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, unconjugated x y^T version.
void	blas::geru (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, unconjugated x y^T version.
void	blas::geru (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>, unconjugated x y^T version.
void	blas::geru (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>, unconjugated x y^T version.

Detailed Description

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

Function Documentation

geru()

template<typename TA , typename TX , typename TY >

void blas::geru	(	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^T + 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]	layout	Matrix storage, Layout::ColMajor or Layout::RowMajor.
[in]	m	Number of rows of the matrix A. m >= 0.
[in]	n	Number of columns of the matrix A. n >= 0.
[in]	alpha	Scalar alpha. If alpha is zero, A is not updated.
[in]	x	The m-element vector x, in an array of length (m-1)*abs(incx) + 1.
[in]	incx	Stride 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]	y	The n-element vector y, in an array of length (n-1)*abs(incy) + 1.
[in]	incy	Stride 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]	A	The m-by-n matrix A, stored in an lda-by-n array [RowMajor: m-by-lda].
[in]	lda	Leading dimension of A. lda >= max(1, m) [RowMajor: lda >= max(1, n)].