From 1d78b3eb222ffd9ccb728ae1dda914c97ea4888f Mon Sep 17 00:00:00 2001 From: wolf Date: Wed, 27 Oct 2004 18:46:11 +0000 Subject: [PATCH] Fix a bug in GridGenerator::cylinder<3> git-svn-id: https://svn.dealii.org/trunk@9727 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/grid/grid_generator.cc | 4 ++-- deal.II/doc/news/c-5.0.html | 14 ++++++++++++-- 2 files changed, 14 insertions(+), 4 deletions(-) diff --git a/deal.II/deal.II/source/grid/grid_generator.cc b/deal.II/deal.II/source/grid/grid_generator.cc index a11ad36fe3..d87aa11418 100644 --- 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 11eb90cd8c..95def95d42 100644 --- 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.
  1. - Improved: Now the FE_Q class supports elements of arbitrary - polynomial degree also in 3D. + Fixed: The GridGenerator::cylinder 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. +
    + (Ralf Schulz 2004/10/27) +

    + +
  2. + Improved: Now the FE_Q class supports hanging node + constraints for elements of arbitrary polynomial degree also in + 3D.
    (Oliver Kayser-Herold 2004/10/21)

    -- 2.39.5