From: Martin Kronbichler Date: Mon, 23 May 2022 18:35:47 +0000 (+0200) Subject: Fix bug in FEEvaluation without vectorization X-Git-Tag: v9.4.0-rc1~165^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F13797%2Fhead;p=dealii.git Fix bug in FEEvaluation without vectorization --- 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) {