diff --git a/dune/tectonic/spatial-solving/tnnmg/linesearchsolver.hh b/dune/tectonic/spatial-solving/tnnmg/linesearchsolver.hh
index 2b9e86dc06ea38c304d38d67c58679c87a485b37..695c3e9a5d44d3bbc95dcae7d868fec2aa2c54a1 100644
--- a/dune/tectonic/spatial-solving/tnnmg/linesearchsolver.hh
+++ b/dune/tectonic/spatial-solving/tnnmg/linesearchsolver.hh
@@ -24,10 +24,10 @@ class LineSearchSolver
     Dune::Solvers::Interval<double> D;
     D = f.subDifferential(0);
 
-    std::cout << "f.A " << f.quadraticPart() << " f.b " << f.linearPart() << std::endl;
+    //std::cout << "f.A " << f.quadraticPart() << " f.b " << f.linearPart() << std::endl;
 
-    std::cout << D[0] << " " << D[1] << std::endl;
-    std::cout << "domain: " << f.domain()[0] << " " << f.domain()[1] << std::endl;
+    //std::cout << D[0] << " " << D[1] << std::endl;
+    //std::cout << "domain: " << f.domain()[0] << " " << f.domain()[1] << std::endl;
 
 
     if (D[1] > 0) // NOTE: Numerical instability can actually get us here
@@ -44,9 +44,9 @@ class LineSearchSolver
     const Bisection globalBisection; //(0.0, 1.0, 0.0, 0.0);
 
     x = globalBisection.minimize(f, f.scaling(), 0.0, bisectionsteps);
-    std::cout << "x: " << x << "scaling: " << f.scaling();
+    //std::cout << "x: " << x << "scaling: " << f.scaling();
     x /= f.scaling();
-    std::cout << "final x: "  << x << std::endl;
+    //std::cout << "final x: "  << x << std::endl;
     //x = f.domain().projectIn(x);
   }
 };