Fix do_reinit of FEPointEvaluation for empty unit_points

dealii · May 4, 2023 · 702793e · 702793e
1 parent 72a4a92
commit 702793e
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/include/deal.II/matrix_free/fe_point_evaluation.h b/include/deal.II/matrix_free/fe_point_evaluation.h
@@ -1498,6 +1498,12 @@ FEPointEvaluation<n_components, dim, spacedim, Number>::do_reinit()
     n_q_points_scalar / n_lanes_user_interface +
     (n_q_points_scalar % n_lanes_user_interface > 0 ? 1 : 0);
 
+  if (n_q_points == 0)
+    {
+      is_reinitialized = true;
+      return;
+    }
+
   // set unit point pointer
   const unsigned int unit_point_offset =
     mapping_info->compute_unit_point_index_offset(current_cell_index,