From 766601375873ce573664f5e0881ff8fbeb950fd6 Mon Sep 17 00:00:00 2001 From: Pratheek Shanthraj Date: Fri, 15 Aug 2014 09:04:00 +0000 Subject: [PATCH] use subF instead of F for analytic jacobian calculation --- code/crystallite.f90 | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/code/crystallite.f90 b/code/crystallite.f90 index 4a9f4ca95..fa65c3796 100644 --- a/code/crystallite.f90 +++ b/code/crystallite.f90 @@ -1133,8 +1133,7 @@ subroutine crystallite_stressAndItsTangent(updateJaco) dFedF = 0.0_pReal forall(p=1_pInt:3_pInt, o=1_pInt:3_pInt) & dFedF(1:3,1:3,p,o) = temp_3333(1:3,1:3,p,o) + & - math_mul33x33(math_mul33x33(crystallite_Fe(1:3,1:3,g,i,e), & - crystallite_Fp(1:3,1:3,g,i,e)), & + math_mul33x33(crystallite_subF(1:3,1:3,g,i,e), & dFpinvdF(1:3,1:3,p,o)) forall(p=1_pInt:3_pInt, o=1_pInt:3_pInt) &