From: wolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Thu, 16 Aug 2001 11:20:42 +0000 (+0000)
Subject: Reinsert an assertion.
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=c0fc5bc206147a8468d6f5608880b6ffa3844de8;p=dealii-svn.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<dim>::reinit (const typename DoFHandler<dim>::cell_iterator
   Assert (subface_no < GeometryInfo<dim>::subfaces_per_face,
 	  ExcIndexRange (subface_no, 0, GeometryInfo<dim>::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);