From 4b870bb7ce1ac960b7379f31c396cb3c141959c9 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Tue, 28 Mar 2006 23:18:21 +0000 Subject: [PATCH] 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 --- .../include/numerics/vectors.templates.h | 23 ++++++++++++++++--- 1 file changed, 20 insertions(+), 3 deletions(-) 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); -- 2.39.5