diff --git a/src/multi-body-problem.cc b/src/multi-body-problem.cc
index 5e4c505e45d55ca572de3a4532334194ab5360b5..f986549e8b66d38b6e2c619e706dfe1be2e16500 100644
--- a/src/multi-body-problem.cc
+++ b/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++) {
diff --git a/src/spatial-solving/tnnmg/tnnmgstep.hh b/src/spatial-solving/tnnmg/tnnmgstep.hh
index d09f3fd1650beae1d6f9dcdf30a754b676a4b7d8..41780c072d6c9a20af92a264f42589209c561c96 100644
--- a/src/spatial-solving/tnnmg/tnnmgstep.hh
+++ b/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