]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Generalised tensor linfty_norm() on number type
authorJean-Paul Pelteret <jppelteret@gmail.com>
Wed, 16 Aug 2017 08:57:03 +0000 (02:57 -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 35037bd3e9a11f43e90f64e03be52473b948f232..a52f0720b5dd76df6de1a2bb632459488b6d4d9c 100644 (file)
@@ -2033,13 +2033,13 @@ l1_norm (const Tensor<2,dim,Number> &t)
  */
 template <int dim, typename Number>
 inline
-double
+Number
 linfty_norm (const Tensor<2,dim,Number> &t)
 {
-  double max = 0;
+  Number max = internal::NumberType<Number>::value(0.0);
   for (unsigned int i=0; i<dim; ++i)
     {
-      double sum = 0;
+      Number sum = internal::NumberType<Number>::value(0.0);
       for (unsigned int j=0; j<dim; ++j)
         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.