]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Add comment.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Mon, 29 Jan 2001 14:54:22 +0000 (14:54 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Mon, 29 Jan 2001 14:54:22 +0000 (14:54 +0000)
git-svn-id: https://svn.dealii.org/trunk@3822 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/grid/grid_generator.cc

index 03f4454bb4a9fa7bb2a2cac1ba43ad0aaad9df0b..d3bb2715e3932160f0e9bfa4411b5023302d40d6 100644 (file)
@@ -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),

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.