Remove a template parameter specifying the number type of a tensor.

This will prevent premature casting of one tensor type into another;
this operation is now taken up in the more general Tensor::operator*.

diff --git a/source/fe/fe_values.cc b/source/fe/fe_values.cc
index 686d6c35c98999aa8f745addd5bf550bef664292..8d73901763e6cd30ab56f6fc2394c16b9fb91a4d 100644 (file)
--- a/source/fe/fe_values.cc
+++ b/source/fe/fe_values.cc
@@ -2892,7 +2892,7 @@ namespace internal
           = &shape_derivatives[shape_func][0];
         for (unsigned int point=0; point<n_quadrature_points; ++point)
           derivatives[point] += value *
-                                dealii::Tensor<order,spacedim,Number>(*shape_derivative_ptr++);
+                                dealii::Tensor<order,spacedim>(*shape_derivative_ptr++);
       }
   }
 
@@ -2965,11 +2965,11 @@ namespace internal
               if (quadrature_points_fastest)
                 for (unsigned int point=0; point<n_quadrature_points; ++point)
                   derivatives[comp][point] += value *
-                                              dealii::Tensor<order,spacedim,Number>(*shape_derivative_ptr++);
+                                              dealii::Tensor<order,spacedim>(*shape_derivative_ptr++);
               else
                 for (unsigned int point=0; point<n_quadrature_points; ++point)
                   derivatives[point][comp] += value *
-                                              dealii::Tensor<order,spacedim,Number>(*shape_derivative_ptr++);
+                                              dealii::Tensor<order,spacedim>(*shape_derivative_ptr++);
             }
           else
             for (unsigned int c=0; c<n_components; ++c)
@@ -2987,11 +2987,11 @@ namespace internal
                 if (quadrature_points_fastest)
                   for (unsigned int point=0; point<n_quadrature_points; ++point)
                     derivatives[comp][point] += value *
-                                                dealii::Tensor<order,spacedim,Number>(*shape_derivative_ptr++);
+                                                dealii::Tensor<order,spacedim>(*shape_derivative_ptr++);
                 else
                   for (unsigned int point=0; point<n_quadrature_points; ++point)
                     derivatives[point][comp] += value *
-                                                dealii::Tensor<order,spacedim,Number>(*shape_derivative_ptr++);
+                                                dealii::Tensor<order,spacedim>(*shape_derivative_ptr++);
               }
         }
   }