diff --git a/dune/elasticity/materials/mooneyrivlindensity.hh b/dune/elasticity/materials/mooneyrivlindensity.hh
index e067923871b8e0be80c3ae078e55a734b2397f74..0e7a76658c6be1d2db1bb9b2ac19eef0d6a7389e 100644
--- a/dune/elasticity/materials/mooneyrivlindensity.hh
+++ b/dune/elasticity/materials/mooneyrivlindensity.hh
@@ -92,7 +92,10 @@ public:
     field_type strainEnergy = 0;
 
     if (mooneyrivlin_energy == "ciarlet")
-      return mooneyrivlin_a*normFSquared + mooneyrivlin_b*normFinvSquared*detF + mooneyrivlin_c*detF*detF - ((dim-1)*mooneyrivlin_a + mooneyrivlin_b + 2*mooneyrivlin_c)*std::log(detF);
+    {
+      using std::log;
+      return mooneyrivlin_a*normFSquared + mooneyrivlin_b*normFinvSquared*detF + mooneyrivlin_c*detF*detF - ((dim-1)*mooneyrivlin_a + mooneyrivlin_b + 2*mooneyrivlin_c)*log(detF);
+    }
     else {
       strainEnergy = mooneyrivlin_10 * trCTildeMinus3 +
                         mooneyrivlin_01 * c2TildeMinus3  +
@@ -104,7 +107,8 @@ public:
                         mooneyrivlin_12 * trCTildeMinus3 * c2TildeMinus3  * c2TildeMinus3  +
                         mooneyrivlin_03 * c2TildeMinus3  * c2TildeMinus3  * c2TildeMinus3;
       if (mooneyrivlin_energy == "log") {
-        field_type logDetF = std::log(detF);
+        using std::log;
+        field_type logDetF = log(detF);
         return strainEnergy + 0.5 * mooneyrivlin_k* logDetF * logDetF;
       } else if (mooneyrivlin_energy == "square") {
         field_type detFMinus1 = detF - 1;
diff --git a/dune/elasticity/materials/neohookedensity.hh b/dune/elasticity/materials/neohookedensity.hh
index 1541e3c25554e72542b62a2e1760df614ee070c2..f93acf04c5b0672a4af73b6dd5e8016e875a9017 100644
--- a/dune/elasticity/materials/neohookedensity.hh
+++ b/dune/elasticity/materials/neohookedensity.hh
@@ -50,15 +50,17 @@ public:
     Dune::FieldVector<field_type, dim> sigmaSquared;
     FMatrixHelp::eigenValues(C, sigmaSquared);
 
+    using std::sqrt;
+
     // singular values of F, i.e., eigenvalues of U
     std::array<field_type, dim> sigma;
     for (int i = 0; i < dim; i++)
-      sigma[i] = std::sqrt(sigmaSquared[i]);
+      sigma[i] = sqrt(sigmaSquared[i]);
 
     field_type detC = 1.0;
     for (int i = 0; i < dim; i++)
       detC *= sigmaSquared[i];
-    field_type detF = std::sqrt(detC);
+    field_type detF = sqrt(detC);
 
     // \tilde{C} = \tilde{F}^T\tilde{F} = \frac{1}{\det{F}^{2/3}}C
     field_type trCTilde = 0;