From: Martin Kronbichler <kronbichler@lnm.mw.tum.de>
Date: Wed, 29 Mar 2017 09:51:49 +0000 (+0200)
Subject: Disable scatter as long as the indices are not fixed
X-Git-Tag: v9.0.0-rc1~1768^2~1
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=87410e3216186b240f71e15fdc390394617640a7;p=dealii.git

Disable scatter as long as the indices are not fixed
---

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<Number> &res,
                              internal::bool2type<true>) 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<VectorizedArray<Number>::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<Number> tmp;
       tmp.gather(vec.begin(), indices);
       tmp += res;