Try to fix valence

diff --git a/source/distributed/solution_transfer.cc b/source/distributed/solution_transfer.cc
index eeb35c3954b463475f893b910805b869f3787a6c..1307bc5117f71b837787efab5dda148d6d4072e9 100644 (file)
--- a/source/distributed/solution_transfer.cc
+++ b/source/distributed/solution_transfer.cc
@@ -267,9 +267,12 @@ namespace parallel
       if (average_values)
         {
           // finalize valence: compress and invert
+          using Number = typename VectorType::value_type;
           valence.compress(::dealii::VectorOperation::add);
           for (const auto i : valence.locally_owned_elements())
-            valence[i] = (valence[i] == 0.0 ? 0.0 : (1.0 / valence[i]));
+            valence[i] = (static_cast<Number>(valence[i]) == Number() ?
+                            Number() :
+                            (Number(1.0) / static_cast<Number>(valence[i])));
           valence.compress(::dealii::VectorOperation::insert);
 
           for (const auto vec : all_out)