Added dXdF as input to update_extras! (#11)

Co-authored-by: Jacob Gunnarsson <jacob.gunnarsson@chalmers.se>

Author: GunnarssonJ

src/FiniteStrainPlastic.jl

@@ -122,7 +122,6 @@ function MMB.material_response(m::FiniteStrainPlastic, F::Tensor{2,3}, old::Fini
         if converged
             x_sol = fromvector(x_vector, x0)
             check_solution(x_sol)
-            update_extras!(extras, x_sol, ∂R∂X) # Use dRdx before calling inv!
             # dPdF = ∂P∂F + ∂P∂X dXdF
             # dRdF = 0 = ∂R∂F + ∂R∂X dXdF => dXdF = - ∂R∂X\∂R∂F
 
@@ -131,6 +130,8 @@ function MMB.material_response(m::FiniteStrainPlastic, F::Tensor{2,3}, old::Fini
 
             dXdF = - ∂R∂X \ ∂R∂F
 
+            update_extras!(extras, x_sol, ∂R∂X, dXdF)
+
             P = calculate_PKstress(m, x_sol, old, F)
             P_X_fun = Tensor2VectorFun(x_sol, x -> calculate_PKstress(m, x, old, F))
             ∂P∂X = ForwardDiff.jacobian(P_X_fun, tovector(SVector, x_sol))