11#ifndef TLAPACK_BLAS_GER_HH
12#define TLAPACK_BLAS_GER_HH
37 class T = type_t<matrixA_t>,
38 disable_if_allow_optblas_t<pair<alpha_t, T>,
41 pair<vectorY_t, T> > = 0>
52 const idx_t m = nrows(
A);
53 const idx_t n = ncols(
A);
59 for (idx_t j = 0; j < n; ++j) {
61 for (idx_t i = 0; i < m; ++i)
62 A(i, j) +=
x[i] *
tmp;
66#ifdef TLAPACK_USE_LAPACKPP
72 class T = type_t<matrixA_t>,
73 enable_if_allow_optblas_t<pair<alpha_t, T>,
76 pair<vectorY_t, T> > = 0>
77void ger(
const alpha_t alpha,
83 auto A_ = legacy_matrix(A);
84 auto x_ = legacy_vector(x);
85 auto y_ = legacy_vector(y);
88 constexpr Layout L = layout<matrixA_t>;
92 return ::blas::ger((::blas::Layout)L, m, n, alpha, x_.ptr, x_.inc, y_.ptr,
93 y_.inc, A_.ptr, A_.ldim);
constexpr T conj(const T &x) noexcept
Extends std::conj() to real datatypes.
Definition utils.hpp:100
#define TLAPACK_SCALAR
Macro for tlapack::concepts::Scalar compatible with C++17.
Definition concepts.hpp:915
#define TLAPACK_LEGACY_VECTOR
Macro for tlapack::concepts::LegacyVector compatible with C++17.
Definition concepts.hpp:954
#define TLAPACK_LEGACY_MATRIX
Macro for tlapack::concepts::LegacyMatrix compatible with C++17.
Definition concepts.hpp:951
#define TLAPACK_VECTOR
Macro for tlapack::concepts::Vector compatible with C++17.
Definition concepts.hpp:906
#define TLAPACK_MATRIX
Macro for tlapack::concepts::Matrix compatible with C++17.
Definition concepts.hpp:896
void ger(const alpha_t &alpha, const vectorX_t &x, const vectorY_t &y, matrixA_t &A)
General matrix rank-1 update:
Definition ger.hpp:42
#define tlapack_check_false(cond)
Throw an error if cond is true.
Definition exceptionHandling.hpp:113
typename traits::real_type_traits< Types..., int >::type real_type
The common real type of the list of types.
Definition scalar_type_traits.hpp:113