diff --git a/dune/gfe/averagedistanceassembler.hh b/dune/gfe/averagedistanceassembler.hh
index b4def6ba1165b4aac40784f58edd8049e2c078f9..3c23d23285b334c81971459c86a6cbd50b3af045 100644
--- a/dune/gfe/averagedistanceassembler.hh
+++ b/dune/gfe/averagedistanceassembler.hh
@@ -69,21 +69,6 @@ public:
         orthonormalFrame.mv(embeddedGradient,gradient);
     }
 
-    void assembleEmbeddedHessian(const TargetSpace& x,
-                         Dune::FieldMatrix<ctype,embeddedSize,embeddedSize>& matrix) const
-    {
-        matrix = 0;
-        for (size_t i=0; i<coefficients_.size(); i++) {
-          Dune::SymmetricMatrix<ctype,embeddedSize> tmp = TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[i], x);
-          for (int j=0; j<embeddedSize; j++)
-            for (int k=0; k<=j; k++) {
-              matrix[j][k] += weights_[i] * tmp(j,k);
-              if (j!=k)
-                matrix[k][j] += weights_[i] * tmp(j,k);
-            }
-        }
-    }
-
     void assembleEmbeddedHessian(const TargetSpace& x,
                          Dune::SymmetricMatrix<ctype,embeddedSize>& matrix) const
     {
@@ -92,6 +77,7 @@ public:
             matrix.axpy(weights_[i], TargetSpace::secondDerivativeOfDistanceSquaredWRTSecondArgument(coefficients_[i], x));
     }
 
+    // TODO Use a Symmetric matrix for the result
     void assembleHessian(const TargetSpace& x,
                          Dune::FieldMatrix<ctype,size,size>& matrix) const
     {