diff --git a/MaterialLib/SolidModels/Ehlers-impl.h b/MaterialLib/SolidModels/Ehlers-impl.h
index 191607bdaa2a58c6263d925dc3d7e799f839514a..2b0367c5db5e4a3f820902f7e90765b00d29b86d 100644
--- a/MaterialLib/SolidModels/Ehlers-impl.h
+++ b/MaterialLib/SolidModels/Ehlers-impl.h
@@ -621,8 +621,8 @@ newton(double const dt, MaterialProperties const& mp,
     };
 
     auto const update_jacobian = [&](JacobianMatrix& jacobian) {
-        calculatePlasticJacobian<DisplacementDim>(
-            dt, jacobian, s, solution[KelvinVectorSize * 2 + 2], mp);
+        jacobian = calculatePlasticJacobian<DisplacementDim>(
+            dt, s, solution[KelvinVectorSize * 2 + 2], mp);
     };
 
     auto const update_solution = [&](ResidualVectorType const& increment) {