diff --git a/dune/gfe/surfacecosseratstressassembler.hh b/dune/gfe/surfacecosseratstressassembler.hh
index b8246b8fa47a2826bbc63a915bd21924a4fdb15c..61e9b4434ef854702542ee55a448c325bb26337d 100644
--- a/dune/gfe/surfacecosseratstressassembler.hh
+++ b/dune/gfe/surfacecosseratstressassembler.hh
@@ -96,7 +96,6 @@ namespace Dune::GFE {
             // Compute gradients of Base functions
             for (size_t i=0; i<gradients.size(); ++i)
               gradients[i] = referenceGradients[i] * transpose(jacobianInverseTransposed);
-              //jacobianInverseTransposed.mv(referenceGradients[i][0], gradients[i]);
 
             // Deformation gradient in vector form
             size_t nDoubles = dim*dim;
@@ -249,8 +248,6 @@ namespace Dune::GFE {
               for (size_t i=0; i<gradients.size(); i++)
                 gradients[i] = referenceGradients[i] * transpose(jacobianInverseTransposed);
 
-                //jacobianInverseTransposed.mv(referenceGradients[i][0], gradients[i]);
-
               // Deformation gradient - call this U_es_minus_Id already
               FieldMatrix<double,dim,dim> U_es_minus_Id(0);
               for (size_t i=0; i<gradients.size(); i++)