diff --git a/dune/gfe/geodesicfeassembler.hh b/dune/gfe/geodesicfeassembler.hh
index bf1593c94c53608f9b850f05e103a78a9e627e07..2a13170ad9641e9c25ee4f8fe624b2637fa0d38d 100644
--- a/dune/gfe/geodesicfeassembler.hh
+++ b/dune/gfe/geodesicfeassembler.hh
@@ -126,7 +126,7 @@ assembleMatrix(const std::vector<TargetSpace>& sol,
             localSolution[i] = sol[indexSet.subIndex(*it,i,gridDim)];
 
         // setup matrix 
-        localStiffness_->assemble(*it, localSolution);
+        localStiffness_->assembleHessian(*it, localSolution);
 
         // Add element matrix to global stiffness matrix
         for(int i=0; i<numOfBaseFct; i++) { 
diff --git a/dune/gfe/localgeodesicfestiffness.hh b/dune/gfe/localgeodesicfestiffness.hh
index ae0786d009d51f719f6cc9cb478e33d67d7bfc92..f2a2239cb8a5897dc5f3d080ca09dc5c8792debc 100644
--- a/dune/gfe/localgeodesicfestiffness.hh
+++ b/dune/gfe/localgeodesicfestiffness.hh
@@ -102,7 +102,7 @@ public:
 
     This default implementation used finite-difference approximations to compute the second derivatives
     */
-    virtual void assemble(const Entity& e,
+    virtual void assembleHessian(const Entity& e,
                   const std::vector<TargetSpace>& localSolution);
    
     virtual RT energy (const Entity& e,
@@ -158,7 +158,7 @@ assembleGradient(const Entity& element,
 
 template <class GridView, class TargetSpace>
 void LocalGeodesicFEStiffness<GridView,TargetSpace>::
-assemble(const Entity& element,
+assembleHessian(const Entity& element,
          const std::vector<TargetSpace>& localSolution)
 {
     // 1 degree of freedom per element vertex
@@ -338,7 +338,7 @@ public:
 
     This default implementation used finite-difference approximations to compute the second derivatives
     */
-    virtual void assemble(const Entity& e,
+    virtual void assembleHessian(const Entity& e,
                   const std::vector<TargetSpace>& localSolution);
     
     virtual RT energy (const Entity& e,
@@ -484,7 +484,7 @@ assembleGradient(const Entity& element,
 // ///////////////////////////////////////////////////////////
 template <class GridType, int dim>
 void LocalGeodesicFEStiffness<GridType,UnitVector<dim> >::
-assemble(const Entity& element,
+assembleHessian(const Entity& element,
          const std::vector<TargetSpace>& localSolution)
 {
     // 1 degree of freedom per element vertex