diff --git a/src/rodassembler.cc b/src/rodassembler.cc
index b30a80940c6465912397b98ed00f9696edc3a127..b1a86a94594c35a81bd3f180d8d5954edaf67201 100644
--- a/src/rodassembler.cc
+++ b/src/rodassembler.cc
@@ -214,39 +214,6 @@ getResultantForce(const BoundaryPatchBase<PatchGridView>& boundary,
 }
 
 
-template <class GridView>
-void PlanarRodAssembler<GridView>::getNeighborsPerVertex(Dune::MatrixIndexSet& nb) const
-{
-    const int gridDim = GridView::dimension;
-    const typename GridView::IndexSet& indexSet = this->gridView_.indexSet();
-    
-    int i, j;
-    int n = this->gridView_.size(gridDim);
-    
-    nb.resize(n, n);
-    
-    ElementIterator it    = this->gridView_.template begin<0>();
-    ElementIterator endit = this->gridView_.template end<0>  ();
-    
-    for (; it!=endit; ++it) {
-        
-        for (i=0; i<it->template count<gridDim>(); i++) {
-            
-            for (j=0; j<it->template count<gridDim>(); j++) {
-                
-                int iIdx = indexSet.subIndex(*it,i,gridDim);
-                int jIdx = indexSet.subIndex(*it,j,gridDim);
-                
-                nb.add(iIdx, jIdx);
-                
-            }
-            
-        }
-        
-    }
-    
-}
-
 template <class GridView>
 void PlanarRodAssembler<GridView>::
 assembleMatrix(const std::vector<RigidBodyMotion<2> >& sol,
@@ -255,7 +222,7 @@ assembleMatrix(const std::vector<RigidBodyMotion<2> >& sol,
     const typename GridView::IndexSet& indexSet = this->gridView_.indexSet();
 
     Dune::MatrixIndexSet neighborsPerVertex;
-    getNeighborsPerVertex(neighborsPerVertex);
+    this->getNeighborsPerVertex(neighborsPerVertex);
     
     matrix = 0;
     
diff --git a/src/rodassembler.hh b/src/rodassembler.hh
index af2b47d6e1b381c50525b3d8a11d2631e4140e40..f5c67235f2d4450c52e36573e7a6faf977ee5d16 100644
--- a/src/rodassembler.hh
+++ b/src/rodassembler.hh
@@ -132,8 +132,6 @@ public:
     /** \brief Compute the energy of a deformation state */
     double computeEnergy(const std::vector<RigidBodyMotion<2> >& sol) const;
     
-    void getNeighborsPerVertex(Dune::MatrixIndexSet& nb) const;
-    
 protected:
     
     /** \brief Compute the element tangent stiffness matrix  */