diff --git a/dune/tectonic/ellipticenergy.hh b/dune/tectonic/ellipticenergy.hh
index 39aa3f0faa5b72d5206eb74472ee21b5dc61c351..d64b8ab123d864462ac1c2c8ae318334e27c8922 100644
--- a/dune/tectonic/ellipticenergy.hh
+++ b/dune/tectonic/ellipticenergy.hh
@@ -68,16 +68,10 @@ template <int dim> class EllipticEnergy {
     double const mgx = mg * x;
     if (pgx >= 0 && mgx >= 0) {
       ret = pg;
-      dverb << "## Directional derivative (as per scalar product w/ "
-               "semigradient): " << -(ret * mg)
-            << " (coordinates of the restriction)" << std::endl;
       ret *= -1;
       return true;
     } else if (pgx <= 0 && mgx <= 0) {
       ret = mg;
-      dverb << "## Directional derivative (as per scalar product w/ "
-               "semigradient): " << -(ret * pg)
-            << " (coordinates of the restriction)" << std::endl;
       ret *= -1;
       return true;
     } else {
@@ -86,9 +80,6 @@ template <int dim> class EllipticEnergy {
       SmallVector d;
       smoothGradient(x, d);
       complementaryProjection(d, x, ret);
-      dverb << "## Directional derivative (as per scalar product w/ "
-               "semigradient): " << -(ret * ret)
-            << " (coordinates of the restriction)" << std::endl;
       ret *= -1;
       return false;
     }