diff --git a/MaterialLib/SolidModels/Ehlers-impl.h b/MaterialLib/SolidModels/Ehlers-impl.h
index f166ef19b6d7c69b3a832a123886f72d64ac3fae..191607bdaa2a58c6263d925dc3d7e799f839514a 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) {