// $Id$
// Version: $Name$
//
-// Copyright (C) 2006, 2007 by the deal.II authors
+// Copyright (C) 2006, 2007, 2008 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
//---------------------------- kelly_crash_01.cc ---------------------------
-// a modified version of step-27 that crashes because of circular constraints.
+// a modified version of step-27 that crashes due to circular constraints
char logname[] = "kelly_crash_01/output";
for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
if ((cell->face(f)->center()[2] != -4)
&&
- (cell->face(f)->center()[2] != 7))
+ (cell->face(f)->center()[2] != 7)
+ &&
+ (cell->face(f)->at_boundary()))
cell->face(f)->set_boundary_indicator (1);
triangulation.refine_global (1);
//---------------------------- kelly_crash_02.cc ---------------------------
-// like kelly_crash_01 but much reduced. it turned out that
-// cell->face(f)->at_boundary() and cell->at_boundary(f) did not always return
-// the same thing, although they of course should. as a result, the
-// KellyErrorEstimator forgot to work on certain faces
+// like kelly_crash_01 (which has since been renamed circular_01 because it
+// continued to fail, though for other reasons) but much reduced. it turned
+// out that cell->face(f)->at_boundary() and cell->at_boundary(f) did not
+// always return the same thing, although they of course should. as a result,
+// the KellyErrorEstimator forgot to work on certain faces
//
// the problem turned out to be that we were setting a boundary indicator for
// an interior face. while cell->face(f)->at_boundary() checks for a boundary