diff --git a/src/rodassembler.cc b/src/rodassembler.cc
index 47c1d20b563efc5587d9dee07707c95b57c0e622..c89266cafff69258eb47368ede61af158d344fb8 100644
--- a/src/rodassembler.cc
+++ b/src/rodassembler.cc
@@ -5,7 +5,6 @@
 
 #include <dune/grid/common/quadraturerules.hh>
 
-#include <dune/disc/shapefunctions/lagrangeshapefunctions.hh>
 #include <dune/localfunctions/lagrange/p1.hh>
 
 #include "src/rodlocalstiffness.hh"
@@ -118,9 +117,8 @@ getStrain(const std::vector<RigidBodyMotion<3> >& sol,
         int elementIdx = indexSet.index(*it);
 
         // Extract local solution on this element
-        const LagrangeShapeFunctionSet<double, double, gridDim> & baseSet 
-            = Dune::LagrangeShapeFunctions<double, double, gridDim>::general(it->type(), elementOrder);
-        int numOfBaseFct = baseSet.size();
+        Dune::P1LocalFiniteElement<double,double,gridDim> localFiniteElement;
+        int numOfBaseFct = localFiniteElement.localCoefficients().size();
 
         std::vector<RigidBodyMotion<3> > localSolution(2);