Merge pull request #13920 from peterrum/valence_fix

Try to fix valence
dealii · Jun 7, 2022 · f262a26 · f262a26
2 parents 35bcb43 + 95fbdfd
commit f262a26
Showing 1 changed file with 4 additions and 1 deletion.
diff --git 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)