From c67d7b9e086ba484a814000dc20fa288d3593885 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Wed, 5 May 2021 09:50:54 -0600 Subject: [PATCH] Add an assertion. --- include/deal.II/numerics/vector_tools_boundary.templates.h | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/include/deal.II/numerics/vector_tools_boundary.templates.h b/include/deal.II/numerics/vector_tools_boundary.templates.h index 5f52784853..b82c6f7fb2 100644 --- a/include/deal.II/numerics/vector_tools_boundary.templates.h +++ b/include/deal.II/numerics/vector_tools_boundary.templates.h @@ -1027,6 +1027,7 @@ namespace VectorTools // Get boundary function values // at quadrature points. + AssertDimension(boundary_function.n_components, fe.n_components()); boundary_function.vector_value_list(quadrature_points, values); const std::vector> &reference_quadrature_points = @@ -1219,6 +1220,7 @@ namespace VectorTools // Get boundary function // values at quadrature // points. + AssertDimension(boundary_function.n_components, fe.n_components()); boundary_function.vector_value_list(quadrature_points, values); switch (dim) @@ -2133,6 +2135,7 @@ namespace VectorTools // Get boundary function values // at quadrature points. + AssertDimension(boundary_function.n_components, fe.n_components()); boundary_function.vector_value_list(quadrature_points, values); // Find the group of vector components we want to project onto @@ -2440,6 +2443,7 @@ namespace VectorTools Vector(fe.n_components())); // Get boundary function values at quadrature points. + AssertDimension(boundary_function.n_components, fe.n_components()); boundary_function.vector_value_list(quadrature_points, values); // Find where the group of vector components (dim of them, -- 2.39.5