11#ifndef TLAPACK_BLAS_HER_HH
12#define TLAPACK_BLAS_HER_HH
60 const idx_t n = nrows(
A);
68 for (idx_t j = 0; j < n; ++j) {
70 for (idx_t i = 0; i < j; ++i)
71 A(i, j) +=
x[i] *
tmp;
77 for (idx_t j = 0; j < n; ++j) {
81 for (idx_t i = j + 1; i < n; ++i)
82 A(i, j) +=
x[i] *
tmp;
87#ifdef TLAPACK_USE_LAPACKPP
99 auto A_ = legacy_matrix(
A);
100 auto x_ = legacy_vector(
x);
104 const auto& n =
A_.n;
106 return ::blas::her((::blas::Layout)L, (::blas::Uplo)
uplo, n,
alpha,
x_.ptr,
#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
#define TLAPACK_REAL
Macro for tlapack::concepts::Real compatible with C++17.
Definition concepts.hpp:918
void her(Uplo uplo, const alpha_t &alpha, const vectorX_t &x, matrixA_t &A)
Hermitian matrix rank-1 update:
Definition her.hpp:53
#define tlapack_check_false(cond)
Throw an error if cond is true.
Definition exceptionHandling.hpp:113
Concept for types that represent tlapack::Uplo.
Sort the numbers in D in increasing order (if ID = 'I') or in decreasing order (if ID = 'D' ).
Definition arrayTraits.hpp:15
enable_if_t<(!allow_optblas< T1, Ts... >), int > disable_if_allow_optblas_t
Disable if the list of types allows optimized BLAS library.
Definition utils.hpp:303
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 real_type< T > real(const T &x) noexcept
Extends std::real() to real datatypes.
Definition utils.hpp:71
constexpr T conj(const T &x) noexcept
Extends std::conj() to real datatypes.
Definition utils.hpp:100
constexpr real_type< T > imag(const T &x) noexcept
Extends std::imag() to real datatypes.
Definition utils.hpp:86
Uplo
Definition types.hpp:50
@ Upper
0 <= i <= j, 0 <= j <= n.
@ Lower
0 <= i <= m, 0 <= j <= i.
Layout
Definition types.hpp:29