]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Generalised tensor l1_norm() on number type
authorJean-Paul Pelteret <jppelteret@gmail.com>
Wed, 16 Aug 2017 08:56:36 +0000 (02:56 -0600)
committerJean-Paul Pelteret <jppelteret@gmail.com>
Mon, 21 Aug 2017 16:53:30 +0000 (18:53 +0200)
include/deal.II/base/tensor.h

index d6982e69b33516f7d51b64a8250add0fc10077a7..35037bd3e9a11f43e90f64e03be52473b948f232 100644 (file)
@@ -2006,13 +2006,13 @@ cofactor (const Tensor<2,dim,Number> &t)
  */
 template <int dim, typename Number>
 inline
-double
+Number
 l1_norm (const Tensor<2,dim,Number> &t)
 {
-  double max = 0;
+  Number max = internal::NumberType<Number>::value(0.0);
   for (unsigned int j=0; j<dim; ++j)
     {
-      double sum = 0;
+      Number sum = internal::NumberType<Number>::value(0.0);
       for (unsigned int i=0; i<dim; ++i)
         sum += std::fabs(t[i][j]);
 

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.