Change type of P1LocalFiniteElement to PQkLocalFiniteElementCache<..,..,..,..>::FiniteElementType

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::PQkLocalFiniteElementCache<double, double, 1, 1>::FiniteElementType localFiniteElement;
         int numOfBaseFct = localFiniteElement.localCoefficients().size();
 
         std::vector<RigidBodyMotion<double,3> > localSolution(2);