diff --git a/dune/tectonic/mydirectionalconvexfunction.hh b/dune/tectonic/mydirectionalconvexfunction.hh
index 9f147ae43b0e87f2b11af9cecc4576473530ef91..389eb4104b2b76ed07b996560dccc42919af5775 100644
--- a/dune/tectonic/mydirectionalconvexfunction.hh
+++ b/dune/tectonic/mydirectionalconvexfunction.hh
@@ -4,9 +4,6 @@
 #ifndef MY_DIRECTIONAL_CONVEX_FUNCTION_HH
 #define MY_DIRECTIONAL_CONVEX_FUNCTION_HH
 
-// just for debugging
-#include <dune/solvers/computeenergy.hh>
-
 #include <dune/fufem/interval.hh>
 
 /*
@@ -43,13 +40,6 @@ template <class NonlinearityType> class MyDirectionalConvexFunction {
     phi.directionalDomain(u, v, dom);
   }
 
-  /* Just for debugging */
-  double operator()(double x) const {
-    VectorType tmp = v;
-    tmp *= x;
-    return computeEnergy(A, x, b) + phi(tmp);
-  }
-
   double quadraticPart() const { return A; }
 
   double linearPart() const { return b; }