diff --git a/dune/elasticity/common/elasticityhelpers.hh b/dune/elasticity/common/elasticityhelpers.hh
index 279fdafab14b11b408f98b9d5ba2c9f1bb63872a..4decd0675a924d2de557bf9e290b42c80e69813b 100644
--- a/dune/elasticity/common/elasticityhelpers.hh
+++ b/dune/elasticity/common/elasticityhelpers.hh
@@ -5,6 +5,8 @@
 
 #include <dune/common/fmatrix.hh>
 
+#include <dune/fufem/symmetrictensor.hh>
+
 namespace Dune {
 
     namespace Elasticity {
@@ -133,7 +135,7 @@ namespace Dune {
             linDet = 0;
             for (int i=0; i<2; i++)                       
                 for (int j=(i+1)%3; j<3; j++) {
-                    int k=(j+1)%3 + (j+1)/3;          
+                    int k=(-(i+j))%3;
                     linDet[i][j] = u[j][k]*u[k][i] - u[j][i]*(1+u[k][k]);
                     linDet[j][i] = u[k][j]*u[i][k] - u[i][j]*(1+u[k][k]);
             }    
@@ -213,7 +215,7 @@ namespace Dune {
             // 1 1 1 0   -->  0  - 0
             // 1 1 1 1   -->  0  - 0
         }
-        
+
         double det_val(const Dune::FieldMatrix<double,3,3>& u) {
             return (1+u[0][0])*(1+u[1][1])*(1+u[2][2])
                 + u[0][1]*u[1][2]*u[2][0]