[MaL] KelvinVector; Relax deviatoric assertion.

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);
 }