11#ifndef TLAPACK_BLAS_GER_HH
12#define TLAPACK_BLAS_GER_HH
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
83 auto A_ = legacy_matrix(
A);
84 auto x_ = legacy_vector(
x);
85 auto y_ = legacy_vector(
y);
92 return ::blas::ger((::blas::Layout)L, m, n,
alpha,
x_.ptr,
x_.inc,
y_.ptr,
#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
Sort the numbers in D in increasing order (if ID = 'I') or in decreasing order (if ID = 'D' ).
Definition arrayTraits.hpp:15
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
constexpr T conj(const T &x) noexcept
Extends std::conj() to real datatypes.
Definition utils.hpp:100
Layout
Definition types.hpp:29