diff --git a/src/localgeodesicfestiffness.hh b/src/localgeodesicfestiffness.hh
index 42aa5744657cf3345058ec517ee3367f28195100..9b084c4a5ce6a0cec8e59693fc8300e3cd7a4688 100644
--- a/src/localgeodesicfestiffness.hh
+++ b/src/localgeodesicfestiffness.hh
@@ -5,9 +5,8 @@
 #include <dune/common/fmatrix.hh>
 #include <dune/istl/matrixindexset.hh>
 #include <dune/istl/matrix.hh>
-#include <dune/disc/operators/localstiffness.hh>
-#include<dune/disc/operators/boundaryconditions.hh>
 
+#include "localstiffness.hh"
 #include "rigidbodymotion.hh"
 
 template<class GridView, class TargetSpace>
@@ -145,12 +144,6 @@ assemble(const Entity& element,
 
     this->A = 0;
 
-    for (int i=0; i<nDofs; i++) {
-        this->b[i] = 0;
-        for (int j=0; j<this->bctype[i].size(); j++)
-            this->bctype[i][j] = Dune::BoundaryConditions::neumann;
-    }
-
     double eps = 1e-4;
 
     typedef typename Dune::Matrix<Dune::FieldMatrix<double,blocksize,blocksize> >::row_type::iterator ColumnIterator;