Ensure that DoFs not touched by VectorTools::interpolate don't get set.

diff --git a/include/deal.II/numerics/vector_tools.templates.h b/include/deal.II/numerics/vector_tools.templates.h
index 6eac0117ad39d56cc6772e04217db019ee431b54..b9f6b2f47c031c418e644e44b130cc9fd6c9e3e5 100644 (file)
--- a/include/deal.II/numerics/vector_tools.templates.h
+++ b/include/deal.II/numerics/vector_tools.templates.h
@@ -515,14 +515,21 @@ namespace VectorTools
 
       for (const auto i : interpolation.locally_owned_elements())
         {
-          const auto value =
-            ::dealii::internal::ElementAccess<VectorType>::get(interpolation,
-                                                               i);
           const auto weight =
             ::dealii::internal::ElementAccess<VectorType>::get(weights, i);
-          ::dealii::internal::ElementAccess<VectorType>::set(value / weight,
-                                                             i,
-                                                             vec);
+
+          // See if we touched this DoF at all. If so, set the average
+          // of the value we computed in the output vector. Otherwise,
+          // don't touch the value at all.
+          if (weight != number(0))
+            {
+              const auto value =
+                ::dealii::internal::ElementAccess<VectorType>::get(
+                  interpolation, i);
+              ::dealii::internal::ElementAccess<VectorType>::set(value / weight,
+                                                                 i,
+                                                                 vec);
+            }
         }
       vec.compress(VectorOperation::insert);
     }