From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Tue, 28 Mar 2006 23:18:21 +0000 (+0000)
Subject: Use FiniteElementData<dim>::face_to_equivalent_cell_index to avoid an
X-Git-Tag: v8.0.0~11976
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=4b870bb7ce1ac960b7379f31c396cb3c141959c9;p=dealii.git

Use FiniteElementData<dim>::face_to_equivalent_cell_index to avoid an
unjustified exception.


git-svn-id: https://svn.dealii.org/trunk@12736 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/deal.II/include/numerics/vectors.templates.h b/deal.II/deal.II/include/numerics/vectors.templates.h
index 961563d4de..c1e57726a5 100644
--- a/deal.II/deal.II/include/numerics/vectors.templates.h
+++ b/deal.II/deal.II/include/numerics/vectors.templates.h
@@ -946,13 +946,30 @@ interpolate_boundary_values (const Mapping<dim>            &mapping,
 				             // initialized by the default
 				             // constructor), since we don't
 				             // need them anyway.
+					     //
+					     // As a detour, we must
+					     // make sure we only
+					     // query
+					     // face_system_to_component_index
+					     // if the index
+					     // corresponds to a
+					     // primitive shape
+					     // function. since we
+					     // know that all the
+					     // components we are
+					     // interested in are
+					     // primitive (by the
+					     // above check), we can
+					     // safely put such a
+					     // check in front
 	    if (unit_support_points.size() == 0)
 	      {
 	        unit_support_points.resize (fe.dofs_per_face);
                 for (unsigned int i=0; i<fe.dofs_per_face; ++i)
-                  if (component_mask[fe.face_system_to_component_index(i).first]
-                      == true)
-		    unit_support_points[i] = fe.unit_face_support_point(i);
+		  if (fe.is_primitive (fe.face_to_equivalent_cell_index(i)))
+		    if (component_mask[fe.face_system_to_component_index(i).first]
+			== true)
+		      unit_support_points[i] = fe.unit_face_support_point(i);
               }
 
             Quadrature<dim-1> aux_quad (unit_support_points);