From 5744cd1e76cbd49f29a61b0e743f04344dbe33b0 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Wed, 12 Dec 2018 16:01:37 +0100
Subject: [PATCH] [MaL] KelvinVector; Relax deviatoric assertion.
---
MathLib/KelvinVector-impl.h | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/MathLib/KelvinVector-impl.h b/MathLib/KelvinVector-impl.h
index 5baca5d3eec..8424ca47e09 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);
}
--
GitLab