From: Wolfgang Bangerth Date: Thu, 16 Aug 2001 11:20:42 +0000 (+0000) Subject: Reinsert an assertion. X-Git-Tag: v8.0.0~18909 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=bdbac370e490c6f919a06fc77119e49b71b8f385;p=dealii.git Reinsert an assertion. git-svn-id: https://svn.dealii.org/trunk@4889 0785d39b-7218-0410-832d-ea1e28bc413d --- 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);