Comment

src/one-body-sample.cc

@@ -283,6 +283,7 @@ int main(int argc, char *argv[]) {
               false); // absolute error
           overallSolver.solve();
 
+          // Update the state
           for (size_t i = 0; i < frictionalNodes.size(); ++i) {
             if (frictionalNodes[i][0]) {
               double const unorm = u_diff[i].two_norm();