From c0fc5bc206147a8468d6f5608880b6ffa3844de8 Mon Sep 17 00:00:00 2001 From: wolf Date: Thu, 16 Aug 2001 11:20:42 +0000 Subject: [PATCH] Reinsert an assertion. git-svn-id: https://svn.dealii.org/trunk@4889 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/fe/fe_values.cc | 13 ++++++++++--- 1 file changed, 10 insertions(+), 3 deletions(-) diff --git a/deal.II/deal.II/source/fe/fe_values.cc b/deal.II/deal.II/source/fe/fe_values.cc index f3653a108c..860f9c3f5c 100644 --- a/deal.II/deal.II/source/fe/fe_values.cc +++ b/deal.II/deal.II/source/fe/fe_values.cc @@ -804,9 +804,16 @@ void FESubfaceValues::reinit (const typename DoFHandler::cell_iterator Assert (subface_no < GeometryInfo::subfaces_per_face, ExcIndexRange (subface_no, 0, GeometryInfo::subfaces_per_face)); -//TODO:[GK] Reinsert this assertion? It tests a necessary, not a sufficient condition - // Assert (cell->face(face_no)->at_boundary() == false, - // ExcReinitCalledWithBoundaryFace()); + // check that the face we are + // presently working on is not at + // the boundary, since boundary + // faces cannot, by definition, be + // further refined as there is no + // other cell behind the face, and + // then there is no point in using + // subfacevalues on these faces + Assert (cell->face(face_no)->at_boundary() == false, + ExcReinitCalledWithBoundaryFace()); present_cell = cell; present_face = cell->face(face_no); -- 2.39.5