From 77f2ee51b1feabb57e735c09ba6a7ef831d36d5f Mon Sep 17 00:00:00 2001 From: heister Date: Fri, 15 Nov 2013 21:10:02 +0000 Subject: [PATCH] assert that we have a enough fe components in compute_no_normal_flux_constraints git-svn-id: https://svn.dealii.org/trunk@31675 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/include/deal.II/numerics/vector_tools.templates.h | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/deal.II/include/deal.II/numerics/vector_tools.templates.h b/deal.II/include/deal.II/numerics/vector_tools.templates.h index 92670d6ac5..258f83a2c9 100644 --- a/deal.II/include/deal.II/numerics/vector_tools.templates.h +++ b/deal.II/include/deal.II/numerics/vector_tools.templates.h @@ -4294,6 +4294,10 @@ namespace VectorTools internal::VectorDoFTuple vector_dofs; vector_dofs.dof_indices[0] = face_dofs[i]; + Assert(first_vector_component+dim<=fe.n_components(), + ExcMessage("Error: the finite element has not enough components " + "to define a normal direction.")); + for (unsigned int k=0; k