diff --git a/dune/gfe/localgeodesicfestiffness.hh b/dune/gfe/localgeodesicfestiffness.hh
index a83533789739c9b035115054f4b477a5bade07e5..ae0786d009d51f719f6cc9cb478e33d67d7bfc92 100644
--- a/dune/gfe/localgeodesicfestiffness.hh
+++ b/dune/gfe/localgeodesicfestiffness.hh
@@ -516,9 +516,6 @@ assemble(const Entity& element,
 
     }
 
-    std::cout << "embeddedHessian" << std::endl;
-    printmatrix(std::cout, embeddedHessian, "embeddedHessian", "--");
-
     // transform to local tangent space bases
     std::vector<Dune::FieldMatrix<double,blocksize,embeddedBlocksize> > orthonormalFrames(nDofs);
     std::vector<Dune::FieldMatrix<double,embeddedBlocksize,blocksize> > orthonormalFramesTransposed(nDofs);
@@ -542,9 +539,6 @@ assemble(const Entity& element,
 
         }
     
-    std::cout << "localHessian" << std::endl;
-    printmatrix(std::cout, A_, "localHessian", "--");
-
 #endif
 }