From 38b6e7f410f59d47b4d43f139a714e098f9ea1d2 Mon Sep 17 00:00:00 2001 From: guido Date: Tue, 1 Jul 2003 07:45:32 +0000 Subject: [PATCH] has_support_on_face returns true by default git-svn-id: https://svn.dealii.org/trunk@7825 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/include/fe/fe_raviart_thomas.h | 17 +++++++---------- deal.II/deal.II/source/fe/fe_raviart_thomas.cc | 10 +++++++--- 2 files changed, 14 insertions(+), 13 deletions(-) diff --git a/deal.II/deal.II/include/fe/fe_raviart_thomas.h b/deal.II/deal.II/include/fe/fe_raviart_thomas.h index ba2a09c954..9d88b11fb0 100644 --- a/deal.II/deal.II/include/fe/fe_raviart_thomas.h +++ b/deal.II/deal.II/include/fe/fe_raviart_thomas.h @@ -255,16 +255,13 @@ class FE_RaviartThomas : public FiniteElement virtual unsigned int element_multiplicity (const unsigned int index) const; /** - * This function returns - * @p{true}, if the shape - * function @p{shape_index} has - * non-zero values on the face - * @p{face_index}. For the lowest - * order Nedelec elements, this - * is actually the case for the - * one on which the shape - * function is defined and all - * neighboring ones. + * Check whether a shape function + * is non-zero on a face. + * + * Right now, this is only + * implemented for RT0 in + * 1D. Otherwise, returns always + * @p{true}. * * Implementation of the * interface in diff --git a/deal.II/deal.II/source/fe/fe_raviart_thomas.cc b/deal.II/deal.II/source/fe/fe_raviart_thomas.cc index ccaafa22d3..51fe053653 100644 --- a/deal.II/deal.II/source/fe/fe_raviart_thomas.cc +++ b/deal.II/deal.II/source/fe/fe_raviart_thomas.cc @@ -1928,6 +1928,9 @@ FE_RaviartThomas::has_support_on_face (const unsigned int shape_index, Assert (face_index < GeometryInfo::faces_per_cell, ExcIndexRange (face_index, 0, GeometryInfo::faces_per_cell)); + // Return computed values if we + // know them easily. Otherwise, it + // is always safe to return true. switch (degree) { case 0: @@ -1946,14 +1949,15 @@ FE_RaviartThomas::has_support_on_face (const unsigned int shape_index, = { 2, 3, 0, 1}; return (face_index != opposite_faces[shape_index]); - }; + } - default: Assert (false, ExcNotImplemented()); + default: + return true; }; }; default: // other degree - Assert (false, ExcNotImplemented()); + return true; }; return true; -- 2.39.5