Revert "Dirichlet constraints can be distributed here"

This reverts commit b0f8edb.

src/matrixfree/computeLHS.cc

@@ -35,8 +35,15 @@ MatrixFreePDE<dim, degree>::vmult(vectorType &dst, const vectorType &src) const
                                  true);
     }
 
-  // Account for Dirichlet BC's
-  constraintsDirichletSet[currentFieldIndex]->distribute(dst);
+  // Account for Dirichlet BC's (essentially copy dirichlet DOF values present in src to
+  // dst, although it is unclear why the constraints can't just be distributed here)
+  for (auto &it : *valuesDirichletSet[currentFieldIndex])
+    {
+      if (dst.in_local_range(it.first))
+        {
+          dst(it.first) = src(it.first);
+        }
+    }
 
   // end log
   computing_timer.leave_subsection("matrixFreePDE: computeLHS");