diff --git a/dune/gfe/localgeodesicfestiffness.hh b/dune/gfe/localgeodesicfestiffness.hh
index 044bc52bb2e50493b9ce1c8c08027e195dbdfc5e..ab5ab42d0c89907bb20cc64de4927f4cec892956 100644
--- a/dune/gfe/localgeodesicfestiffness.hh
+++ b/dune/gfe/localgeodesicfestiffness.hh
@@ -359,12 +359,6 @@ public:
         DUNE_THROW(Dune::NotImplemented, "!");
     }
 
-    /** \todo Remove this once this methods is not in base class LocalStiffness anymore */
-    void assemble (const Entity& e, int k=1)
-    {
-        DUNE_THROW(Dune::NotImplemented, "!");
-    }
-
     virtual RT energy (const Entity& e,
                        const std::vector<TargetSpace>& localSolution) const = 0;
 
diff --git a/dune/gfe/localstiffness.hh b/dune/gfe/localstiffness.hh
index 9dc94e777af01f0c9c29040635606581169a7286..f7a9abf5874a650e1673e87310a771babcdbc4c7 100644
--- a/dune/gfe/localstiffness.hh
+++ b/dune/gfe/localstiffness.hh
@@ -76,26 +76,6 @@ namespace Dune
 	{
 	}
 
-	//! assemble local stiffness matrix including boundary conditions for given element and order
-	/*! On exit the following things have been done:
-	  - The stiffness matrix for the given entity and polynomial degree has been assembled and is
-        accessible with the mat() method.
-	  - The boundary conditions have been evaluated and are accessible with the bc() method. 
-        The boundary conditions are either neumann, process or dirichlet. Neumann indicates
-        that the corresponding node (assuming a nodal basis) is at the Neumann boundary, process
-        indicates that the node is at a process boundary (arising from the parallel decomposition of the mesh).
-        Process boundaries are treated as homogeneous Dirichlet conditions, i.e. the corresponding value
-        in the right hand side is set to 0. Finally, Dirichlet indicates that the node is at the Dirichlet
-        boundary.  
-	  - The right hand side has been assembled. It contains either the value of the essential boundary
-        condition or the assembled source term and neumann boundary condition. 
-		It is accessible via the rhs() method.
-
-	  @param[in]  e    a codim 0 entity reference
-	  @param[in]  k    order of Lagrange basis (default is 1)
-	 */
-      virtual void assemble (const Entity& e, int k=1) = 0;
-
       /** \brief assemble local stiffness matrix including boundary conditions for given element and order
           
       Unlike the method with only two arguments, this one additionally takes the local solution in order