.

src/multi-body-problem.cc

@@ -199,7 +199,7 @@ int main(int argc, char *argv[]) {
      programState.setupInitialConditions(parset, contactNetwork);
 
 
-    DUNE_THROW(Dune::Exception, "Just need to stop here!");
+    //DUNE_THROW(Dune::Exception, "Just need to stop here!");
 
     auto& nBodyAssembler = contactNetwork.nBodyAssembler();
     for (size_t i=0; i<bodyCount; i++) {

src/spatial-solving/tnnmg/tnnmgstep.hh

@@ -151,7 +151,7 @@ class TNNMGStep :
 
     //print(f.linearPart(), "f.linearPart():");
 
-    //std::cout << "-- energy before smoothing: " << energy(x) << std::endl;
+    std::cout << "-- energy before smoothing: " << f(x) << std::endl;
 
     // Nonlinear presmoothing
     for (std::size_t i=0; i<preSmoothingSteps_; ++i)
@@ -160,7 +160,7 @@ class TNNMGStep :
     std::cout << "- nonlinear presmoothing: success" << std::endl;
     //print(x, "TNNMG iterate after smoothing:");
 
-    std::cout << "-- energy after presmoothing: " << energy(x) << std::endl;
+    std::cout << "-- energy after presmoothing: " << f(x) << std::endl;
 
     /*
     // Compute constraint/truncated linearization