]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Reinsert an assertion.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 16 Aug 2001 11:20:42 +0000 (11:20 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 16 Aug 2001 11:20:42 +0000 (11:20 +0000)
git-svn-id: https://svn.dealii.org/trunk@4889 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe_values.cc

index f3653a108c6379d32970ccf16c8cba82a755f1ee..860f9c3f5cf36434079bce801a3befa892fad0c2 100644 (file)
@@ -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);

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.