Explicitly initialise tensor contraction results to zero.

diff --git a/include/deal.II/base/tensor_accessors.h b/include/deal.II/base/tensor_accessors.h
index 52c7ca7bf1ab576d94d57122ad90dec485b33ce8..2ca2674e47f7068fcd4168abb428ded36b3c09d0 100644 (file)
--- a/include/deal.II/base/tensor_accessors.h
+++ b/include/deal.II/base/tensor_accessors.h
@@ -665,7 +665,9 @@ namespace TensorAccessors
       inline DEAL_II_ALWAYS_INLINE static
       T1 contract2(const T2 &left, const T3 &right)
       {
-        T1 result = T1();
+        // Some auto-differentiable numbers need explicit
+        // zero initialization.
+        T1 result = dealii::internal::NumberType<T1>::value(0.0);
         for (unsigned int i = 0; i < dim; ++i)
           result += Contract2<no_contr - 1, dim>::template contract2<T1>(left[i], right[i]);
         return result;
@@ -710,7 +712,9 @@ namespace TensorAccessors
       static inline
       T1 contract3(const T2 &left, const T3 &middle, const T4 &right)
       {
-        T1 result = T1();
+        // Some auto-differentiable numbers need explicit
+        // zero initialization.
+        T1 result = dealii::internal::NumberType<T1>::value(0.0);
         for (unsigned int i = 0; i < dim; ++i)
           result += Contract3<rank_1 - 1, rank_2, dim>::template contract3<T1>(left[i], middle[i], right);
         return result;
@@ -736,7 +740,9 @@ namespace TensorAccessors
       static inline
       T1 contract3(const T2 &left, const T3 &middle, const T4 &right)
       {
-        T1 result = T1();
+        // Some auto-differentiable numbers need explicit
+        // zero initialization.
+        T1 result = dealii::internal::NumberType<T1>::value(0.0);
         for (unsigned int i = 0; i < dim; ++i)
           result += Contract3<0, rank_2 - 1, dim>::template contract3<T1>(left, middle[i], right[i]);
         return result;