diff --git a/dune/gfe/rodlocalstiffness.hh b/dune/gfe/rodlocalstiffness.hh
index c215606277aa95e56c815b1be69c44286a90d70e..cc43e0e000cd55d4ad28f9cc5c6c4de1331c2e73 100644
--- a/dune/gfe/rodlocalstiffness.hh
+++ b/dune/gfe/rodlocalstiffness.hh
@@ -465,7 +465,7 @@ getStrain(const std::vector<RigidBodyMotion<RT,3> >& localSolution,
         // multiply with jacobian inverse
         Dune::FieldVector<double,1> tmp(0);
         inv.umv(shapeGrad[dof][0], tmp);
-        shapeGrad[dof] = tmp;
+        shapeGrad[dof][0] = tmp;
 
     }