11#ifndef TLAPACK_BLAS_SYR2K_HH
12#define TLAPACK_BLAS_SYR2K_HH
92 for (idx_t j = 0; j < n; ++j) {
93 for (idx_t i = 0; i <= j; ++i)
96 for (idx_t
l = 0;
l <
k; ++
l) {
99 for (idx_t i = 0; i <= j; ++i)
105 for (idx_t j = 0; j < n; ++j) {
106 for (idx_t i = j; i < n; ++i)
109 for (idx_t
l = 0;
l <
k; ++
l) {
112 for (idx_t i = j; i < n; ++i)
123 for (idx_t j = 0; j < n; ++j) {
124 for (idx_t i = 0; i <= j; ++i) {
127 for (idx_t
l = 0;
l <
k; ++
l) {
136 for (idx_t j = 0; j < n; ++j) {
137 for (idx_t i = j; i < n; ++i) {
140 for (idx_t
l = 0;
l <
k; ++
l) {
151 for (idx_t j = 0; j < n; ++j) {
152 for (idx_t i = j + 1; i < n; ++i)
158#ifdef TLAPACK_USE_LAPACKPP
193 auto A_ = legacy_matrix(
A);
194 auto B_ = legacy_matrix(
B);
195 auto C_ = legacy_matrix(
C);
199 const auto& n =
C_.n;
205 -3,
"Infs and NaNs in A or B will not propagate to C on output");
209 "Infs and NaNs in C on input will not propagate to C on output");
211 return ::blas::syr2k((::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
void syr2k(Uplo uplo, Op trans, const alpha_t &alpha, const matrixA_t &A, const matrixB_t &B, const beta_t &beta, matrixC_t &C)
Symmetric rank-k update:
Definition syr2k.hpp:64
#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
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
Op
Definition types.hpp:227
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