Fix bug in FEEvaluation for constraints in face integrals

diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h
index 6ecc583b301394f97ae1f540e72660b90a2d175f..c27bfacd719643fc561568c7270efd37bedfa8d3 100644 (file)
--- a/include/deal.II/matrix_free/fe_evaluation.h
+++ b/include/deal.II/matrix_free/fe_evaluation.h
@@ -3899,7 +3899,7 @@ FEEvaluationBase<dim, n_components_, Number, is_face, VectorizedArrayType>::
           Assert(cells[v] < dof_info->row_starts.size() - 1,
                  ExcInternalError());
           has_constraints =
-            has_constraints &&
+            has_constraints ||
             dof_info
                 ->row_starts[cells[v] * n_fe_components +
                              first_selected_component + n_components]