]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Bugfix: Fix misplaced brackets in condition 907/head
authorMatthias Maier <matthias.maier@iwr.uni-heidelberg.de>
Thu, 7 May 2015 08:03:57 +0000 (10:03 +0200)
committerMatthias Maier <matthias.maier@iwr.uni-heidelberg.de>
Thu, 7 May 2015 08:03:57 +0000 (10:03 +0200)
source/grid/grid_generator.cc

index 7e7fa1641c795c462e6bb6738f9827e9ab27b63e..2c16bbd2dc68190187d489301ae75bfce7bdebb5 100644 (file)
@@ -3100,9 +3100,10 @@ namespace GridGenerator
             if (cell->at_boundary(i))
               cell->face(i)->set_all_boundary_ids(2);
 
-        // Next look for the curved boundaries. If the first component is
-        // not near to zero, the boundary is curved. Then decide whether the center
-        // is nearer to the inner or outer boundary to set the correct boundary id.
+        // Next look for the curved boundaries. If the x value of the
+        // center of the face is not equal to center(0), we're on a curved
+        // boundary. Then decide whether the center is nearer to the inner
+        // or outer boundary to set the correct boundary id.
         for (cell=tria.begin(); cell!=tria.end(); ++cell)
           for (unsigned int i=0; i<GeometryInfo<3>::faces_per_cell; ++i)
             if (cell->at_boundary(i))
@@ -3111,7 +3112,7 @@ namespace GridGenerator
                   = cell->face(i);
 
                 const Point<3> face_center (face->center());
-                if (std::abs((face_center(0)-center(0)) > 1.e-3))
+                if (std::abs(face_center(0)-center(0)) > 1.e-6)
                   {
                     if (std::abs((face_center-center).norm()-inner_radius) <
                         std::abs((face_center-center).norm()-outer_radius))

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.