diff --git a/dune/gfe/cosseratenergystiffness.hh b/dune/gfe/cosseratenergystiffness.hh
index dadf9615f880d3403a7d0f982cb2e841749d3e71..51ab17e0b762bb8696f848831bb6e8244cdddef4 100644
--- a/dune/gfe/cosseratenergystiffness.hh
+++ b/dune/gfe/cosseratenergystiffness.hh
@@ -135,7 +135,8 @@ public:
             
         return mu_ * sym(RT_DR3).frobenius_norm2()
                + mu_c_ * skew(RT_DR3).frobenius_norm2()
-               + mu_*lambda_/(2*mu_+lambda_) * traceSquared(RT_DR3);
+               /** \todo Is this sym correct?  It is in the paper, but not in the notes */
+               + mu_*lambda_/(2*mu_+lambda_) * traceSquared(sym(RT_DR3));
     }
 
     /** \brief The shell thickness */