Assert correction for issue #17450

diff --git a/include/deal.II/numerics/vector_tools_interpolate.templates.h b/include/deal.II/numerics/vector_tools_interpolate.templates.h
index f5fafb1488a592a442a442aef404f9fe134a4a65..24ac01d0459ac9401845cde710beda1b094ea36f 100644 (file)
--- a/include/deal.II/numerics/vector_tools_interpolate.templates.h
+++ b/include/deal.II/numerics/vector_tools_interpolate.templates.h
@@ -320,7 +320,11 @@ namespace VectorTools
       hp::QCollection<dim> support_quadrature;
       for (unsigned int fe_index = 0; fe_index < fe.size(); ++fe_index)
         {
-          const auto &fe_i   = fe[fe_index];
+          const auto &fe_i = fe[fe_index];
+          Assert(fe_i.has_generalized_support_points(),
+                 ExcMessage(
+                   "The finite element does not have generalized support "
+                   "points. This is required for interpolation."));
           const auto &points = fe_i.get_generalized_support_points();
           support_quadrature.push_back(Quadrature<dim>(points));
           if (fe_i.n_base_elements() == 1 &&