set dirichlet data for dirichlet body correctly

dune/tectonic/data-structures/program_state.hh

@@ -159,8 +159,9 @@ template <class VectorTEMPLATE, class ScalarVectorTEMPLATE> class ProgramState {
         const auto& body = *contactNetwork.body(i);
 
         if (body.data()->getYoungModulus() == 0.0) {
-            for (; dof<body.nVertices(); dof++) {
+            for (size_t j=0; j<body.nVertices(); j++) {
                 dirichletNodes[dof] = true;
+                dof++;
             }
         } else {
             dof += body.nVertices();