]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Fix bug in FEEvaluation without vectorization 13797/head
authorMartin Kronbichler <martin.kronbichler@uni-a.de>
Mon, 23 May 2022 18:35:47 +0000 (20:35 +0200)
committerMartin Kronbichler <martin.kronbichler@uni-a.de>
Mon, 23 May 2022 18:35:47 +0000 (20:35 +0200)
include/deal.II/matrix_free/evaluation_kernels.h

index 7958f5bf11141e185d273b76e82efc99a9ccd81f..886a8d7b62c2e822b9f4126fcc9255e48522cc93 100644 (file)
@@ -4769,8 +4769,12 @@ namespace internal
                orientation[v][i];
     };
 
+    const unsigned int cell_index =
+      dof_access_index == MatrixFreeFunctions::DoFInfo::dof_access_cell ?
+        fe_eval.get_cell_ids()[0] :
+        cell * n_lanes;
     const unsigned int *dof_indices =
-      &dof_info.dof_indices_contiguous[dof_access_index][cell * n_lanes];
+      &dof_info.dof_indices_contiguous[dof_access_index][cell_index];
 
     for (unsigned int comp = 0; comp < n_components; ++comp)
       {

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.