Commit b3682271 authored by ESTERIE Pierre's avatar ESTERIE Pierre
Browse files

Temporary remove fma in m2l product

parent c8fd3a82
...@@ -229,10 +229,10 @@ namespace scalfmm::tensor ...@@ -229,10 +229,10 @@ namespace scalfmm::tensor
auto t_ = xsimd::load_aligned(&t.data()[i]); auto t_ = xsimd::load_aligned(&t.data()[i]);
auto k_ = xsimd::load_aligned(&k.data()[i]); auto k_ = xsimd::load_aligned(&k.data()[i]);
auto acc_ = xsimd::load_aligned(&accumulate.data()[i]); auto acc_ = xsimd::load_aligned(&accumulate.data()[i]);
auto times = k_*splat; //auto times = k_*splat;
//acc_ += t_ * k_ * splat; //auto tmp = xsimd::fma(t_, times, acc_);
auto tmp = xsimd::fma(t_, times, acc_); acc_ += t_ * k_ * splat;
xsimd::store_aligned(&accumulate.data()[i], tmp); xsimd::store_aligned(&accumulate.data()[i], acc_);
} }
for(std::size_t i{vec_size}; i<size; ++i) for(std::size_t i{vec_size}; i<size; ++i)
{ {
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment