From: heister Date: Fri, 15 Nov 2013 21:10:02 +0000 (+0000) Subject: assert that we have a enough fe components in compute_no_normal_flux_constraints X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=77f2ee51b1feabb57e735c09ba6a7ef831d36d5f;p=dealii-svn.git 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 --- 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