From: Martin Kronbichler <kronbichler@lnm.mw.tum.de>
Date: Fri, 3 Sep 2021 17:46:11 +0000 (+0200)
Subject: FEFaceEvaluation: Do not run into reduced reads with Hessians
X-Git-Tag: v9.4.0-rc1~1003^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=c48b05c9c8b50dac0d0391214713ae307280bd84;p=dealii.git

FEFaceEvaluation: Do not run into reduced reads with Hessians
---

diff --git a/include/deal.II/matrix_free/evaluation_kernels.h b/include/deal.II/matrix_free/evaluation_kernels.h
index d4dd77eafa..8ca4b63eb9 100644
--- a/include/deal.II/matrix_free/evaluation_kernels.h
+++ b/include/deal.II/matrix_free/evaluation_kernels.h
@@ -3499,12 +3499,13 @@ namespace internal
         // to use the functions for at least linear functions for values on
         // the faces and quadratic functions for gradients on the faces, so
         // include the switch here
-        if ((do_gradients == false &&
-             data.data.front().nodal_at_cell_boundaries == true &&
-             fe_degree > 0) ||
-            (data.element_type ==
-               MatrixFreeFunctions::tensor_symmetric_hermite &&
-             fe_degree > 1))
+        if (((do_gradients == false &&
+              data.data.front().nodal_at_cell_boundaries == true &&
+              fe_degree > 0) ||
+             (data.element_type ==
+                MatrixFreeFunctions::tensor_symmetric_hermite &&
+              fe_degree > 1)) &&
+            do_hessians == false)
           {
             // case 1: contiguous and interleaved indices
             if (n_face_orientations == 1 &&
@@ -3792,8 +3793,7 @@ namespace internal
             else if ((n_face_orientations > 1 ||
                       dof_info.index_storage_variants[dof_access_index][cell] ==
                         MatrixFreeFunctions::DoFInfo::IndexStorageVariants::
-                          contiguous) &&
-                     !do_hessians)
+                          contiguous))
               {
                 const unsigned int *indices =
                   &dof_info.dof_indices_contiguous[dof_access_index]