From 3f732fc6df35455275f02b2e6d30f833208e8169 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Mon, 29 Jan 2001 14:54:22 +0000 Subject: [PATCH] Add comment. git-svn-id: https://svn.dealii.org/trunk@3822 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/deal.II/source/grid/grid_generator.cc | 23 +++++++++++++++++-- 1 file changed, 21 insertions(+), 2 deletions(-) diff --git a/deal.II/deal.II/source/grid/grid_generator.cc b/deal.II/deal.II/source/grid/grid_generator.cc index 03f4454bb4..d3bb2715e3 100644 --- a/deal.II/deal.II/source/grid/grid_generator.cc +++ b/deal.II/deal.II/source/grid/grid_generator.cc @@ -536,10 +536,29 @@ void GridGenerator::hyper_L<> (Triangulation<3> &tria, }; + template <> void GridGenerator::hyper_ball<> (Triangulation<3> &tria, const Point<3> &p, - const double radius) { + const double radius) +{ + // this function used to be + // implemented by the code below, + // but it turned out that it didn't + // work as expected: there were + // faces that we used in both + // orientations, leading the + // triangulation function to + // believe that the faces were + // external, even though they were + // in fact internal to the + // ball. this leads to strange + // results. + // + // since we haven't found a working + // enumeration of cell vertices and + // faces, this function is disabled + // altogether Assert(false, ExcNotImplemented()); const double a = 1./(1+sqrt(3)); // equilibrate cell sizes at transition @@ -550,7 +569,7 @@ void GridGenerator::hyper_ball<> (Triangulation<3> &tria, = { // first the vertices of the inner // cell - p+Point<3>(-1,-1,-1)*(radius/sqrt(3)*a), + p+Point<3>(-1,-1,-1)*(radius/sqrt(3)*a), p+Point<3>(+1,-1,-1)*(radius/sqrt(3)*a), p+Point<3>(+1,+1,-1)*(radius/sqrt(3)*a), p+Point<3>(-1,+1,-1)*(radius/sqrt(3)*a), -- 2.39.5