]> https://gitweb.dealii.org/ - dealii.git/commitdiff
FEFaceEvaluation: Do not run into reduced reads with Hessians 12740/head
authorMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Fri, 3 Sep 2021 17:46:11 +0000 (19:46 +0200)
committerMartin Kronbichler <kronbichler@lnm.mw.tum.de>
Fri, 3 Sep 2021 17:48:48 +0000 (19:48 +0200)
include/deal.II/matrix_free/evaluation_kernels.h

index d4dd77eafaa7fc95cf5462b82e96a1dbaeda5e84..8ca4b63eb9afe77f22ca97133f8043466b1e29c7 100644 (file)
@@ -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]

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.