Simplify some code.

diff --git a/include/deal.II/base/tensor_function.templates.h b/include/deal.II/base/tensor_function.templates.h
index d760d22070bf2fe47ffd432d6ca76178d8611c84..11c2f54cb73aec6dede6bc860e53963c5270943d 100644 (file)
--- a/include/deal.II/base/tensor_function.templates.h
+++ b/include/deal.II/base/tensor_function.templates.h
@@ -123,9 +123,8 @@ template <int rank, int dim, typename Number>
 typename TensorFunction<rank, dim, Number>::gradient_type
 ConstantTensorFunction<rank, dim, Number>::gradient(const Point<dim> &) const
 {
-  static const Tensor<rank + 1, dim, Number> zero{};
-
-  return zero;
+  // Return a zero (=default initialized) tensor
+  return {};
 }
 
 
@@ -140,10 +139,10 @@ ConstantTensorFunction<rank, dim, Number>::gradient_list(
   Assert(gradients.size() == points.size(),
          ExcDimensionMismatch(gradients.size(), points.size()));
 
-  static const Tensor<rank + 1, dim, Number> zero{};
-
-  for (unsigned int i = 0; i < gradients.size(); ++i)
-    gradients[i] = zero;
+  // Return an array of zero tensors.
+  std::fill(gradients.begin(),
+            gradients.end(),
+            typename TensorFunction<rank, dim, Number>::gradient_type());
 }