diff --git a/MaterialLib/FractureModels/MohrCoulomb.cpp b/MaterialLib/FractureModels/MohrCoulomb.cpp
index 8cfddaaf8ca41c3e9c60d8c9e027622e3e681a54..48d24bb312bfa6bafc0baa38ee8c202a2ab12752 100644
--- a/MaterialLib/FractureModels/MohrCoulomb.cpp
+++ b/MaterialLib/FractureModels/MohrCoulomb.cpp
@@ -78,7 +78,7 @@ void MohrCoulomb<DisplacementDim>::computeConstitutiveRelation(
     }
 
     // check shear yield function (Fs)
-    double const Fs = std::abs(sigma[0]) + sigma[1] * tan(mat.phi) - mat.c;
+    double const Fs = std::abs(sigma[0]) + sigma[1] * std::tan(mat.phi) - mat.c;
     if (Fs < .0)
     {
         Kep = Ke;
@@ -87,12 +87,12 @@ void MohrCoulomb<DisplacementDim>::computeConstitutiveRelation(
 
     Eigen::VectorXd dFs_dS(2);
     dFs_dS[0] = MathLib::sgn(sigma[0]);
-    dFs_dS[1] = tan(mat.phi);
+    dFs_dS[1] = std::tan(mat.phi);
 
     // plastic potential function: Qs = |tau| + Sn * tan da
     Eigen::VectorXd dQs_dS(2);
     dQs_dS[0] = MathLib::sgn(sigma[0]);
-    dQs_dS[1] = tan(mat.psi);
+    dQs_dS[1] = std::tan(mat.psi);
 
     // plastic multiplier
     Eigen::RowVectorXd const A = dFs_dS.transpose() * Ke / (dFs_dS.transpose() * Ke * dQs_dS);