diff --git a/dune/gfe/nonplanarcosseratshellenergy.hh b/dune/gfe/nonplanarcosseratshellenergy.hh
index 85e97321ea5719bec0602bcb7507fb9c45628922..48506937e535f63b93066fa090d40d7d444c42b2 100644
--- a/dune/gfe/nonplanarcosseratshellenergy.hh
+++ b/dune/gfe/nonplanarcosseratshellenergy.hh
@@ -230,7 +230,7 @@ energy(const typename Basis::LocalView& localView,
 
     for (int alpha=0; alpha<2; alpha++)
       for (int beta=0; beta<2; beta++)
-        c += sqrt(aScalar) * eps[alpha][beta] * Dune::GFE::dyadicProduct(aContravariant[alpha], aContravariant[beta]);
+        c += aScalar * eps[alpha][beta] * Dune::GFE::dyadicProduct(aContravariant[alpha], aContravariant[beta]);
 
     // Second fundamental form
     // The derivative of the normal field
diff --git a/dune/gfe/surfacecosseratenergy.hh b/dune/gfe/surfacecosseratenergy.hh
index 1aad21afc16e0769232786189914aa4772a8305a..b7b5d1226c2e8279c7c44de3bc62be586cd8a63e 100644
--- a/dune/gfe/surfacecosseratenergy.hh
+++ b/dune/gfe/surfacecosseratenergy.hh
@@ -245,7 +245,7 @@ RT energy(const typename Basis::LocalView& localView,
 
       for (int alpha=0; alpha<2; alpha++)
         for (int beta=0; beta<2; beta++)
-          c += sqrt(aScalar) * eps[alpha][beta] * Dune::GFE::dyadicProduct(aContravariant[alpha], aContravariant[beta]);
+          c += aScalar * eps[alpha][beta] * Dune::GFE::dyadicProduct(aContravariant[alpha], aContravariant[beta]);
 
       // Second fundamental form
       // The derivative of the normal field