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