]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix missing std::s and disambiguate call to sqrt(int).
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 3 Jul 2001 11:42:42 +0000 (11:42 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 3 Jul 2001 11:42:42 +0000 (11:42 +0000)
git-svn-id: https://svn.dealii.org/trunk@4803 0785d39b-7218-0410-832d-ea1e28bc413d

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

index ba4e69ee28093bb6257dd8e97929c229a5a73ac8..9f9f0124514805445636c17eb54aa7619706b850 100644 (file)
@@ -656,32 +656,32 @@ GridGenerator::hyper_ball (Triangulation<3> &tria,
                                   // altogether
   Assert(false, ExcNotImplemented());
   
-  const double a = 1./(1+std::sqrt(3)); // equilibrate cell sizes at transition
-                                       // from the inner part to the radial
-                                       // cells
+  const double a = 1./(1+std::sqrt(3.0)); // equilibrate cell sizes at transition
+                                         // from the inner part to the radial
+                                         // cells
   const unsigned int n_vertices = 16;
   const Point<3> vertices[n_vertices]
     = {
                                     // 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),
-      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/std::sqrt(3.0)*a),
+      p+Point<3>(+1,-1,-1)*(radius/std::sqrt(3.0)*a),
+      p+Point<3>(+1,+1,-1)*(radius/std::sqrt(3.0)*a),
+      p+Point<3>(-1,+1,-1)*(radius/std::sqrt(3.0)*a),
+      p+Point<3>(-1,-1,+1)*(radius/std::sqrt(3.0)*a),
+      p+Point<3>(+1,-1,+1)*(radius/std::sqrt(3.0)*a),
+      p+Point<3>(+1,+1,+1)*(radius/std::sqrt(3.0)*a),
+      p+Point<3>(-1,+1,+1)*(radius/std::sqrt(3.0)*a),
                                       // now the eight vertices at
                                       // the outer sphere
-      p+Point<3>(-1,-1,-1)*(radius/sqrt(3)),
-      p+Point<3>(+1,-1,-1)*(radius/sqrt(3)),
-      p+Point<3>(+1,+1,-1)*(radius/sqrt(3)),
-      p+Point<3>(-1,+1,-1)*(radius/sqrt(3)),
-      p+Point<3>(-1,-1,+1)*(radius/sqrt(3)),
-      p+Point<3>(+1,-1,+1)*(radius/sqrt(3)),
-      p+Point<3>(+1,+1,+1)*(radius/sqrt(3)),
-      p+Point<3>(-1,+1,+1)*(radius/sqrt(3))
+      p+Point<3>(-1,-1,-1)*(radius/std::sqrt(3.0)),
+      p+Point<3>(+1,-1,-1)*(radius/std::sqrt(3.0)),
+      p+Point<3>(+1,+1,-1)*(radius/std::sqrt(3.0)),
+      p+Point<3>(-1,+1,-1)*(radius/std::sqrt(3.0)),
+      p+Point<3>(-1,-1,+1)*(radius/std::sqrt(3.0)),
+      p+Point<3>(+1,-1,+1)*(radius/std::sqrt(3.0)),
+      p+Point<3>(+1,+1,+1)*(radius/std::sqrt(3.0)),
+      p+Point<3>(-1,+1,+1)*(radius/std::sqrt(3.0))
       };
 
                                   // one needs to draw the seven cubes to

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.