diff --git a/dune/gfe/cosseratrodenergy.hh b/dune/gfe/cosseratrodenergy.hh
index a7c1386d4a40aebe08dfec1eba33ad9b8fed999e..6ccb14e2a701c89b72aeb6051826d88e75726fe4 100644
--- a/dune/gfe/cosseratrodenergy.hh
+++ b/dune/gfe/cosseratrodenergy.hh
@@ -327,7 +327,7 @@ getStrain(const std::vector<RigidBodyMotion<double,3> >& sol,
         int elementIdx = indexSet.index(element);
 
         // Extract local solution on this element
-        P1LocalFiniteElement<double,double,1> localFiniteElement;
+        Dune::LagrangeSimplexLocalFiniteElement<double, double, 1, 1> localFiniteElement;
         int numOfBaseFct = localFiniteElement.localCoefficients().size();
 
         std::vector<RigidBodyMotion<double,3> > localSolution(2);