]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Disable scatter as long as the indices are not fixed
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Wed, 29 Mar 2017 09:51:49 +0000 (11:51 +0200)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Wed, 29 Mar 2017 09:51:49 +0000 (11:51 +0200)
include/deal.II/matrix_free/fe_evaluation.h

index 143002fc7c200ad9967a27d3dab82049ed23a443..906fe539116b597441a6428caf245aca2cc656d4 100644 (file)
@@ -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;

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.