From 9e9a46d8c927816772f33ad3d06e720b8a6e4082 Mon Sep 17 00:00:00 2001 From: Matthias Maier Date: Thu, 7 Sep 2017 22:28:53 -0500 Subject: [PATCH] gracefully fall-back to Lagrange interpolation if conformity is unknown --- .../deal.II/numerics/vector_tools.templates.h | 28 ++++++++++--------- 1 file changed, 15 insertions(+), 13 deletions(-) diff --git a/include/deal.II/numerics/vector_tools.templates.h b/include/deal.II/numerics/vector_tools.templates.h index 27adfdd4fd..7ae278f396 100644 --- a/include/deal.II/numerics/vector_tools.templates.h +++ b/include/deal.II/numerics/vector_tools.templates.h @@ -242,16 +242,6 @@ namespace VectorTools // element for the correct transformation. switch (fe[fe_index].conforming_space) { - case FiniteElementData::H1: - DEAL_II_FALLTHROUGH; - case FiniteElementData::L2: - // See Monk, Finite Element Methods for Maxwell's Equations, - // p. 77ff, formula (3.74). - // For given mapping F_K: \hat K \to K, we have to transform - // \hat p = p\circ F_K - // i.e., do nothing. - break; - case FiniteElementData::Hcurl: // See Monk, Finite Element Methods for Maxwell's Equations, // p. 77ff, formula (3.76) and Corollary 3.58. @@ -295,10 +285,22 @@ namespace VectorTools } break; + case FiniteElementData::H1: + DEAL_II_FALLTHROUGH; + case FiniteElementData::L2: + // See Monk, Finite Element Methods for Maxwell's Equations, + // p. 77ff, formula (3.74). + // For given mapping F_K: \hat K \to K, we have to transform + // \hat p = p\circ F_K + // i.e., do nothing. + // + break; + default: - Assert(false, - ExcMessage( - "The supplied finite element has an unknown conformity.")); + // In case we deal with an unknown conformity, just assume we + // deal with a Lagrange element and do nothing. + break; + } /*switch*/ FETools::convert_generalized_support_point_values_to_dof_values( -- 2.39.5