]> https://gitweb.dealii.org/ - dealii.git/commitdiff
MatrixFree: Use separate data arrays for two separate evaluators 14153/head
authorMartin Kronbichler <martin.kronbichler@uni-a.de>
Wed, 20 Jul 2022 16:11:38 +0000 (18:11 +0200)
committerMartin Kronbichler <martin.kronbichler@uni-a.de>
Wed, 20 Jul 2022 16:24:49 +0000 (18:24 +0200)
include/deal.II/matrix_free/mapping_info.templates.h

index 81b1849743219930c9aca4ba79e11536dbe746b8..f17fbb1949a9c1ab50b3878c6aa49539bc9f5a8f 100644 (file)
@@ -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)

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.