From 87410e3216186b240f71e15fdc390394617640a7 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Wed, 29 Mar 2017 11:51:49 +0200 Subject: [PATCH] Disable scatter as long as the indices are not fixed --- include/deal.II/matrix_free/fe_evaluation.h | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h index 143002fc7c..906fe53911 100644 --- a/include/deal.II/matrix_free/fe_evaluation.h +++ b/include/deal.II/matrix_free/fe_evaluation.h @@ -2813,11 +2813,14 @@ namespace internal VectorizedArray &res, internal::bool2type) const { -#if DEAL_II_COMPILER_VECTORIZATION_LEVEL < 3 + // TODO: enable scatter path when indices are fixed + + //#if DEAL_II_COMPILER_VECTORIZATION_LEVEL < 3 +#if 1 for (unsigned int v=0; v::n_array_elements; ++v) vector_access(vec, indices[v]) += res[v]; #else - // only use gather in case there is also scatter + // only use gather in case there is also scatter. VectorizedArray tmp; tmp.gather(vec.begin(), indices); tmp += res; -- 2.39.5