11#ifndef TLAPACK_BLAS_HERK_HH
12#define TLAPACK_BLAS_HERK_HH
94 for (idx_t j = 0; j < n; ++j) {
95 for (idx_t i = 0; i < j; ++i)
99 for (idx_t
l = 0;
l <
k; ++
l) {
103 for (idx_t i = 0; i < j; ++i)
111 for (idx_t j = 0; j < n; ++j) {
113 for (idx_t i = j + 1; i < n; ++i)
116 for (idx_t
l = 0;
l <
k; ++
l) {
122 for (idx_t i = j + 1; i < n; ++i)
131 for (idx_t j = 0; j < n; ++j) {
132 for (idx_t i = 0; i < j; ++i) {
134 for (idx_t
l = 0;
l <
k; ++
l)
139 for (idx_t
l = 0;
l <
k; ++
l)
146 for (idx_t j = 0; j < n; ++j) {
147 for (idx_t i = j + 1; i < n; ++i) {
149 for (idx_t
l = 0;
l <
k; ++
l)
154 for (idx_t
l = 0;
l <
k; ++
l)
163 for (idx_t j = 0; j < n; ++j) {
164 for (idx_t i = j + 1; i < n; ++i)
170#ifdef TLAPACK_USE_LAPACKPP
202 auto A_ = legacy_matrix(
A);
203 auto C_ = legacy_matrix(
C);
207 const auto& n =
C_.n;
213 "Infs and NaNs in A will not propagate to C on output");
217 "Infs and NaNs in C on input will not propagate to C on output");
219 return ::blas::herk((::blas::Layout)L, (::blas::Uplo)
uplo,
#define TLAPACK_SCALAR
Macro for tlapack::concepts::Scalar compatible with C++17.
Definition concepts.hpp:915
#define TLAPACK_LEGACY_MATRIX
Macro for tlapack::concepts::LegacyMatrix compatible with C++17.
Definition concepts.hpp:951
#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 herk(Uplo uplo, Op trans, const alpha_t &alpha, const matrixA_t &A, const beta_t &beta, matrixC_t &C)
Hermitian rank-k update:
Definition herk.hpp:68
#define tlapack_check_false(cond)
Throw an error if cond is true.
Definition exceptionHandling.hpp:113
#define tlapack_warning(info, detailedInfo)
Warning handler.
Definition exceptionHandling.hpp:156
Concept for types that represent tlapack::Op.
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
Op
Definition types.hpp:227
@ ConjTrans
conjugate transpose
constexpr real_type< T > imag(const T &x) noexcept
Extends std::imag() to real datatypes.
Definition utils.hpp:86
Uplo
Definition types.hpp:50
@ General
0 <= i <= m, 0 <= j <= n.
@ Upper
0 <= i <= j, 0 <= j <= n.
@ Lower
0 <= i <= m, 0 <= j <= i.
Layout
Definition types.hpp:29
Strong zero type.
Definition StrongZero.hpp:43