]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Use the Tensor member function norm_square() instead of the Point function square... 501/head
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 5 Feb 2015 12:50:14 +0000 (06:50 -0600)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Thu, 5 Feb 2015 12:50:14 +0000 (06:50 -0600)
examples/step-8/step-8.cc
examples/step-9/step-9.cc

index 8625efab2cae84337ccb551f32deb1820d7a1bf9..98e3def95f211e0c7bc1c85794463ed93e1a4ed0 100644 (file)
@@ -1,6 +1,6 @@
 /* ---------------------------------------------------------------------
  *
- * Copyright (C) 2000 - 2014 by the deal.II authors
+ * Copyright (C) 2000 - 2015 by the deal.II authors
  *
  * This file is part of the deal.II library.
  *
@@ -210,15 +210,15 @@ namespace Step8
     // If now the point <code>p</code> is in a circle (sphere) of radius 0.2
     // around one of these points, then set the force in x-direction to one,
     // otherwise to zero:
-    if (((p-point_1).square() < 0.2*0.2) ||
-        ((p-point_2).square() < 0.2*0.2))
+    if (((p-point_1).norm_square() < 0.2*0.2) ||
+        ((p-point_2).norm_square() < 0.2*0.2))
       values(0) = 1;
     else
       values(0) = 0;
 
     // Likewise, if <code>p</code> is in the vicinity of the origin, then set
     // the y-force to 1, otherwise to zero:
-    if (p.square() < 0.2*0.2)
+    if (p.norm_square() < 0.2*0.2)
       values(1) = 1;
     else
       values(1) = 0;
index 24939609bb04355739313dcca2908d202a078612..76dcaad9e99cc71472b2acc7db8cb6aa1146c4e7 100644 (file)
@@ -1,6 +1,6 @@
 /* ---------------------------------------------------------------------
  *
- * Copyright (C) 2000 - 2014 by the deal.II authors
+ * Copyright (C) 2000 - 2015 by the deal.II authors
  *
  * This file is part of the deal.II library.
  *
@@ -329,7 +329,7 @@ namespace Step9
   {
     Assert (component == 0, ExcIndexRange (component, 0, 1));
     const double diameter = 0.1;
-    return ( (p-center_point).square() < diameter*diameter ?
+    return ( (p-center_point).norm_square() < diameter*diameter ?
              .1/std::pow(diameter,dim) :
              0);
   }
@@ -376,8 +376,8 @@ namespace Step9
   {
     Assert (component == 0, ExcIndexRange (component, 0, 1));
 
-    const double sine_term = std::sin(16*numbers::PI*std::sqrt(p.square()));
-    const double weight    = std::exp(-5*p.square()) / std::exp(-5.);
+    const double sine_term = std::sin(16*numbers::PI*std::sqrt(p.norm_square()));
+    const double weight    = std::exp(-5*p.norm_square()) / std::exp(-5.);
     return sine_term * weight;
   }
 
@@ -1254,8 +1254,8 @@ namespace Step9
         // two cells. Note that as opposed to the introduction, we denote
         // by <code>y</code> the normalized difference vector, as this is
         // the quantity used everywhere in the computations.
-        Point<dim>   y        = neighbor_center - this_center;
-        const double distance = std::sqrt(y.square());
+        Tensor<1,dim> y        = neighbor_center - this_center;
+        const double  distance = y.norm();
         y /= distance;
 
         // Then add up the contribution of this cell to the Y matrix...
@@ -1300,7 +1300,7 @@ namespace Step9
     // using this quantity and the right powers of the mesh width:
     const Tensor<2,dim> Y_inverse = invert(Y);
 
-    Point<dim> gradient;
+    Tensor<1,dim> gradient;
     contract (gradient, Y_inverse, projected_gradient);
 
     // The last part of this function is the one where we
@@ -1311,7 +1311,7 @@ namespace Step9
     // difficult:
     *(std_cxx11::get<1>(cell.iterators)) = (std::pow(std_cxx11::get<0>(cell.iterators)->diameter(),
                                                      1+1.0*dim/2) *
-                                            std::sqrt(gradient.square()));
+                                            std::sqrt(gradient.norm_square()));
 
   }
 }

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.