From: wolf Date: Tue, 28 Mar 2006 23:18:21 +0000 (+0000) Subject: Use FiniteElementData::face_to_equivalent_cell_index to avoid an X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=2d63e9652652d7dd25d9d5c599203025c81ef126;p=dealii-svn.git Use FiniteElementData::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 &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 aux_quad (unit_support_points);