diff options
author | Gael Guennebaud <g.gael@free.fr> | 2019-02-07 15:21:35 +0100 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2019-02-07 15:21:35 +0100 |
commit | b3c4344a6852e55c849976dd46ec4e861399bf16 (patch) | |
tree | 54656766450e12fba109e0535f55d7b35f2fecfa /bench | |
parent | 3091c0389808216a1892444e16209ed411c4093a (diff) |
bug #1676: workaround GCC's bug in c++17 mode.
Diffstat (limited to 'bench')
-rw-r--r-- | bench/bench_gemm.cpp | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/bench/bench_gemm.cpp b/bench/bench_gemm.cpp index 688d99c4a..7c6dbea61 100644 --- a/bench/bench_gemm.cpp +++ b/bench/bench_gemm.cpp @@ -112,6 +112,7 @@ void matlab_cplx_cplx(const M& ar, const M& ai, const M& br, const M& bi, M& cr, cr.noalias() -= ai * bi; ci.noalias() += ar * bi; ci.noalias() += ai * br; + // [cr ci] += [ar ai] * br + [-ai ar] * bi } void matlab_real_cplx(const M& a, const M& br, const M& bi, M& cr, M& ci) @@ -240,7 +241,7 @@ int main(int argc, char ** argv) blas_gemm(a,b,r); c.noalias() += a * b; if(!r.isApprox(c)) { - std::cout << r - c << "\n"; + std::cout << (r - c).norm() << "\n"; std::cerr << "Warning, your product is crap!\n\n"; } #else @@ -249,7 +250,7 @@ int main(int argc, char ** argv) gemm(a,b,c); r.noalias() += a.cast<Scalar>() .lazyProduct( b.cast<Scalar>() ); if(!r.isApprox(c)) { - std::cout << r - c << "\n"; + std::cout << (r - c).norm() << "\n"; std::cerr << "Warning, your product is crap!\n\n"; } } |