diff --git a/MathLib/KelvinVector-impl.h b/MathLib/KelvinVector-impl.h
index 5baca5d3eece4149d2043adb381e5153c51cf436..8424ca47e09bf67c94667e7e4980965638323b17 100644
--- a/MathLib/KelvinVector-impl.h
+++ b/MathLib/KelvinVector-impl.h
@@ -18,7 +18,7 @@ double Invariants<KelvinVectorSize>::equivalentStress(
     Eigen::Matrix<double, KelvinVectorSize, 1> const& deviatoric_v)
 {
     assert(std::abs(trace(deviatoric_v)) <=
-           5e-14 * diagonal(deviatoric_v).norm());
+           2e-13 * diagonal(deviatoric_v).norm());
     return std::sqrt(3 * J2(deviatoric_v));
 }
 
@@ -34,7 +34,7 @@ double Invariants<KelvinVectorSize>::J2(
     Eigen::Matrix<double, KelvinVectorSize, 1> const& deviatoric_v)
 {
     assert(std::abs(trace(deviatoric_v)) <=
-           5e-14 * diagonal(deviatoric_v).norm());
+           2e-13 * diagonal(deviatoric_v).norm());
     return 0.5 * deviatoric_v.transpose() * deviatoric_v;
 }
 
@@ -45,7 +45,7 @@ double Invariants<KelvinVectorSize>::J3(
     Eigen::Matrix<double, KelvinVectorSize, 1> const& deviatoric_v)
 {
     assert(std::abs(trace(deviatoric_v)) <=
-           5e-14 * diagonal(deviatoric_v).norm());
+           2e-13 * diagonal(deviatoric_v).norm());
     return determinant(deviatoric_v);
 }