MatrixFree: Use separate data arrays for two separate evaluators

dealii · Jul 20, 2022 · 9688af2 · 9688af2
1 parent b8f23f1
commit 9688af2
Showing 1 changed file with 4 additions and 5 deletions.
diff --git a/include/deal.II/matrix_free/mapping_info.templates.h b/include/deal.II/matrix_free/mapping_info.templates.h
@@ -2101,11 +2101,10 @@ namespace internal
         FEEvaluationData<dim, VectorizedDouble, true> eval_ext(shape_info,
                                                                false);
 
-        // Let both evaluators use the same array as their use will not
-        // overlap
-        AlignedVector<VectorizedDouble> evaluation_data;
-        eval_int.set_data_pointers(&evaluation_data, dim);
-        eval_ext.set_data_pointers(&evaluation_data, dim);
+        AlignedVector<VectorizedDouble> evaluation_data_int,
+          evaluation_data_ext;
+        eval_int.set_data_pointers(&evaluation_data_int, dim);
+        eval_ext.set_data_pointers(&evaluation_data_ext, dim);
 
         for (unsigned int face = begin_face; face < end_face; ++face)
           for (unsigned vv = 0; vv < n_lanes; vv += n_lanes_d)