diff --git a/MathLib/KelvinVector-impl.h b/MathLib/KelvinVector-impl.h
index 7a860b8869dcf3a3d97db53701644e04f6f00816..884d279f048243a51d7187a180025b519a6094f8 100644
--- a/MathLib/KelvinVector-impl.h
+++ b/MathLib/KelvinVector-impl.h
@@ -19,7 +19,7 @@ double Invariants<KelvinVectorSize>::equivalentStress(
     Eigen::Matrix<double, KelvinVectorSize, 1> const& deviatoric_v)
 {
     assert(std::abs(trace(deviatoric_v)) <=
-           2e-13 * diagonal(deviatoric_v).norm());
+           2e-12 * diagonal(deviatoric_v).norm());
     return std::sqrt(3 * J2(deviatoric_v));
 }
 
@@ -35,7 +35,7 @@ double Invariants<KelvinVectorSize>::J2(
     Eigen::Matrix<double, KelvinVectorSize, 1> const& deviatoric_v)
 {
     assert(std::abs(trace(deviatoric_v)) <=
-           2e-13 * diagonal(deviatoric_v).norm());
+           2e-12 * diagonal(deviatoric_v).norm());
     return 0.5 * deviatoric_v.transpose() * deviatoric_v;
 }
 
@@ -46,7 +46,7 @@ double Invariants<KelvinVectorSize>::J3(
     Eigen::Matrix<double, KelvinVectorSize, 1> const& deviatoric_v)
 {
     assert(std::abs(trace(deviatoric_v)) <=
-           2e-13 * diagonal(deviatoric_v).norm());
+           2e-12 * diagonal(deviatoric_v).norm());
     return determinant(deviatoric_v);
 }