From 209367df0e83e508e3a1e010ba3e41dc542a8af0 Mon Sep 17 00:00:00 2001 From: Franz Roters Date: Fri, 23 Oct 2020 15:41:51 +0200 Subject: [PATCH] save one inversion --- src/crystallite.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/crystallite.f90 b/src/crystallite.f90 index 05c409baa..82c485625 100644 --- a/src/crystallite.f90 +++ b/src/crystallite.f90 @@ -418,9 +418,9 @@ function crystallite_stress() crystallite_subF(1:3,1:3,c,i,e) = crystallite_subF0(1:3,1:3,c,i,e) & + crystallite_subStep(c,i,e) *( crystallite_partitionedF (1:3,1:3,c,i,e) & -crystallite_partitionedF0(1:3,1:3,c,i,e)) - crystallite_Fe(1:3,1:3,c,i,e) = matmul(matmul(crystallite_subF(1:3,1:3,c,i,e), & - math_inv33(crystallite_Fp(1:3,1:3,c,i,e))), & - math_inv33(crystallite_Fi(1:3,1:3,c,i,e))) + crystallite_Fe(1:3,1:3,c,i,e) = matmul(crystallite_subF(1:3,1:3,c,i,e), & + math_inv33(matmul(crystallite_Fi(1:3,1:3,c,i,e), & + crystallite_Fp(1:3,1:3,c,i,e)))) crystallite_subdt(c,i,e) = crystallite_subStep(c,i,e) * crystallite_dt(c,i,e) crystallite_converged(c,i,e) = .false. call integrateState(c,i,e)