From 11263ee920df88b6c75da24145b3026b293c760b Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Fri, 7 Oct 2022 10:59:23 +0200 Subject: [PATCH] MatrixFree face evaluation kernels: Fix check for same face type --- include/deal.II/matrix_free/evaluation_kernels.h | 12 ++++-------- 1 file changed, 4 insertions(+), 8 deletions(-) diff --git a/include/deal.II/matrix_free/evaluation_kernels.h b/include/deal.II/matrix_free/evaluation_kernels.h index 8acdac9d1e..a37aa58155 100644 --- a/include/deal.II/matrix_free/evaluation_kernels.h +++ b/include/deal.II/matrix_free/evaluation_kernels.h @@ -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) { -- 2.39.5