diff --git a/dune/gfe/averagedistanceassembler.hh b/dune/gfe/averagedistanceassembler.hh
index 3e072407e0d5bb42abd24d555fb1c8a4a915db3f..f6676be291f4694dfaef50f75a865dbf4cbb1080 100644
--- a/dune/gfe/averagedistanceassembler.hh
+++ b/dune/gfe/averagedistanceassembler.hh
@@ -65,7 +65,7 @@ public:
                 matrix[i][j] = (i==j);
     }
 
-    void assembleHessian(const TargetSpace& x,
+    void assembleEmbeddedHessian(const TargetSpace& x,
                          Dune::FieldMatrix<double,embeddedSize,embeddedSize>& matrix) const
     {
         matrix = 0;
diff --git a/dune/gfe/localgeodesicfefunction.hh b/dune/gfe/localgeodesicfefunction.hh
index b726f97735b24df87ed84ce9c14a0d7e51c71ede..5f2f4a20bcc748f582ad15b125ca7837dcfdaf1b 100644
--- a/dune/gfe/localgeodesicfefunction.hh
+++ b/dune/gfe/localgeodesicfefunction.hh
@@ -228,7 +228,7 @@ evaluateDerivative(const Dune::FieldVector<ctype, dim>& local) const
     AverageDistanceAssembler<TargetSpace> assembler(coefficients_, w);
     
     Dune::FieldMatrix<ctype,embeddedDim,embeddedDim> dFdq(0);
-    assembler.assembleHessian(q,dFdq);
+    assembler.assembleEmbeddedHessian(q,dFdq);
 
     // ////////////////////////////////////
     //   solve the system
@@ -305,7 +305,7 @@ evaluateDerivativeOfGradientWRTCoefficient(const Dune::FieldVector<ctype, dim>&
     AverageDistanceAssembler<TargetSpace> assembler(coefficients_, w);
     
     Dune::FieldMatrix<ctype,embeddedDim,embeddedDim> dFdq(0);
-    assembler.assembleHessian(q,dFdq);
+    assembler.assembleEmbeddedHessian(q,dFdq);
     
    
     Dune::FieldMatrix<ctype,embeddedDim,embeddedDim> mixedDerivative = TargetSpace::secondDerivativeOfDistanceSquaredWRTFirstAndSecondArgument(coefficients_[coefficient], q);