]> https://gitweb.dealii.org/ - dealii.git/commitdiff
MatrixFree face evaluation kernels: Fix check for same face type 14341/head
authorMartin Kronbichler <martin.kronbichler@uni-a.de>
Fri, 7 Oct 2022 08:59:23 +0000 (10:59 +0200)
committerMartin Kronbichler <martin.kronbichler@uni-a.de>
Fri, 7 Oct 2022 09:07:01 +0000 (11:07 +0200)
include/deal.II/matrix_free/evaluation_kernels.h

index 8acdac9d1e86d189fd3757c1e76bd99294a2eb22..a37aa58155c10914efd64d95465cfcfeaa55be2c 100644 (file)
@@ -4207,14 +4207,10 @@ namespace internal
       if (fe_eval.get_dof_access_index() ==
             MatrixFreeFunctions::DoFInfo::dof_access_cell &&
           fe_eval.is_interior_face() == false) // exterior faces in the ECL loop
-        use_vectorization =
-          fe_eval.get_cell_ids()[0] != numbers::invalid_unsigned_int &&
-          std::all_of(fe_eval.get_cell_ids().begin() + 1,
-                      fe_eval.get_cell_ids().end(),
-                      [&](const auto &v) {
-                        return v == fe_eval.get_cell_ids()[0] ||
-                               v == numbers::invalid_unsigned_int;
-                      });
+        for (unsigned int v = 0; v < Number::size(); ++v)
+          if (fe_eval.get_cell_ids()[v] != numbers::invalid_unsigned_int &&
+              fe_eval.get_face_no(v) != fe_eval.get_face_no(0))
+            use_vectorization = false;
 
       if (use_vectorization == false)
         {

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.