From a646b8b18f7bc92bb0af1ac259065a6e27fd31f8 Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <github@naumov.de>
Date: Sun, 14 May 2017 14:01:53 +0200
Subject: [PATCH] [MatL] E: calculatePlasticResidual returns result.
---
MaterialLib/SolidModels/Ehlers-impl.h | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/MaterialLib/SolidModels/Ehlers-impl.h b/MaterialLib/SolidModels/Ehlers-impl.h
index f166ef19b6d..191607bdaa2 100644
--- a/MaterialLib/SolidModels/Ehlers-impl.h
+++ b/MaterialLib/SolidModels/Ehlers-impl.h
@@ -613,12 +613,11 @@ newton(double const dt, MaterialProperties const& mp,
double const k_hardening =
calculateIsotropicHardening(mp.kappa, mp.hardening_coefficient,
solution[KelvinVectorSize * 2 + 1]);
- calculatePlasticResidual<DisplacementDim>(
+ residual = calculatePlasticResidual<DisplacementDim>(
eps_D, eps_V, s,
solution.template segment<KelvinVectorSize>(KelvinVectorSize),
eps_p_D_dot, solution[KelvinVectorSize * 2], eps_p_V_dot,
- eps_p_eff_dot, solution[KelvinVectorSize * 2 + 2], k_hardening, mp,
- residual);
+ eps_p_eff_dot, solution[KelvinVectorSize * 2 + 2], k_hardening, mp);
};
auto const update_jacobian = [&](JacobianMatrix& jacobian) {
--
GitLab