Fix a bug in GridGenerator::cylinder<3>

git-svn-id: https://svn.dealii.org/trunk@9727 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/deal.II/source/grid/grid_generator.cc b/deal.II/deal.II/source/grid/grid_generator.cc
index a11ad36fe3d77c1c584ee2765804fd98cd4cc606..d87aa11418a390c962b2874f92b730475116140a 100644 (file)
--- a/deal.II/deal.II/source/grid/grid_generator.cc
+++ b/deal.II/deal.II/source/grid/grid_generator.cc
@@ -1151,9 +1151,9 @@ GridGenerator::cylinder (Triangulation<3> &tria,
          if (cell->face(i)->boundary_indicator() == 255)
            continue;
          
-         if (cell->face(i)->center()(0) > 1.-1.e-5)
+         if (cell->face(i)->center()(0) > half_length-1.e-5)
            cell->face(i)->set_boundary_indicator(2);
-         else if (cell->face(i)->center()(0) < -1.+1.e-5)
+         else if (cell->face(i)->center()(0) < -half_length+1.e-5)
            cell->face(i)->set_boundary_indicator(1);
        }
       ++cell;

diff --git a/deal.II/doc/news/c-5.0.html b/deal.II/doc/news/c-5.0.html
index 11eb90cd8ce114b3286e6c81c3313235cf7ea888..95def95d42c3f47f693b02fb19dbe91d0f8bd76a 100644 (file)
--- a/deal.II/doc/news/c-5.0.html
+++ b/deal.II/doc/news/c-5.0.html
@@ -360,8 +360,18 @@ inconvenience this causes.
 
 <ol>
   <li> <p>
-       Improved: Now the FE_Q class supports elements of arbitrary
-       polynomial degree also in 3D.
+       Fixed: The <code>GridGenerator::cylinder</code> function in 3d
+       assigned the wrong boundary value to the top and bottom part of
+       the cylinder if the half length of the cylinder was not equal
+       to 1. This is now fixed.
+       <br>
+       (Ralf Schulz 2004/10/27)
+       </p>
+
+  <li> <p>
+       Improved: Now the <code>FE_Q</code> class supports hanging node
+       constraints for elements of arbitrary polynomial degree also in
+       3D.
        <br>
        (Oliver Kayser-Herold 2004/10/21)
        </p>