From 95a0bb3fec1eec5a4ef7270d7546574d198ae06d Mon Sep 17 00:00:00 2001 From: Oreste Marquis Date: Tue, 13 Aug 2024 18:52:14 -0400 Subject: [PATCH] Assert correction for issue #17450 --- .../deal.II/numerics/vector_tools_interpolate.templates.h | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/include/deal.II/numerics/vector_tools_interpolate.templates.h b/include/deal.II/numerics/vector_tools_interpolate.templates.h index f5fafb1488..24ac01d045 100644 --- 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 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(points)); if (fe_i.n_base_elements() == 1 && -- 2.39.5