Add AssertIndexRange for FEPointEvaluation

diff --git a/doc/news/changes/minor/20210605Pauen b/doc/news/changes/minor/20210605Pauen
new file mode 100644 (file)
index 0000000..62252c8
--- /dev/null
+++ b/doc/news/changes/minor/20210605Pauen
@@ -0,0 +1,4 @@
+Improved: An assert has been added in the constructor of
+FEPointEvaluation to prevent accessing non-existing components.
+<br>
+(Judith Pauen, 2021/06/05)

diff --git a/include/deal.II/matrix_free/fe_point_evaluation.h b/include/deal.II/matrix_free/fe_point_evaluation.h
index 8acf9d16c20f4a768cecbbfa451386be8ca7fd12..f059e9b097e61a867c28a8c994c17b4495cf6f14 100644 (file)
--- a/include/deal.II/matrix_free/fe_point_evaluation.h
+++ b/include/deal.II/matrix_free/fe_point_evaluation.h
@@ -691,6 +691,9 @@ FEPointEvaluation<n_components, dim, spacedim, Number>::FEPointEvaluation(
   , update_flags(update_flags)
   , update_flags_mapping(update_default)
 {
+  AssertIndexRange(first_selected_component + n_components,
+                   fe.n_components() + 1);
+
   bool         same_base_element   = true;
   unsigned int base_element_number = 0;
   unsigned int component           = 0;