From 4dc63903bcffdd82e61aa790f3b55c953dd92654 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Mon, 23 May 2022 20:35:47 +0200 Subject: [PATCH] Fix bug in FEEvaluation without vectorization --- include/deal.II/matrix_free/evaluation_kernels.h | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/include/deal.II/matrix_free/evaluation_kernels.h b/include/deal.II/matrix_free/evaluation_kernels.h index 7958f5bf11..886a8d7b62 100644 --- a/include/deal.II/matrix_free/evaluation_kernels.h +++ b/include/deal.II/matrix_free/evaluation_kernels.h @@ -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) { -- 2.39.5