From 584beee94e5b981386774453ff07f9fb09bb725c Mon Sep 17 00:00:00 2001 From: Pratheek Shanthraj Date: Mon, 11 Aug 2014 08:34:44 +0000 Subject: [PATCH] dTdFe bug fix --- code/constitutive.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/code/constitutive.f90 b/code/constitutive.f90 index aefe0586e..bad830b82 100644 --- a/code/constitutive.f90 +++ b/code/constitutive.f90 @@ -496,7 +496,7 @@ subroutine constitutive_hooke_TandItsTangent(T, dT_dFe, Fe, ipc, ip, el) dT_dFe = 0.0_pReal forall (i=1_pInt:3_pInt, j=1_pInt:3_pInt, k=1_pInt:3_pInt, l=1_pInt:3_pInt) & - dT_dFe(i,j,k,l) = sum(C(i,j,l,1:3)*Fe(k,1:3)) ! dT*_ij/dFe_kl + dT_dFe(i,j,k,l) = damage*sum(C(i,j,l,1:3)*Fe(k,1:3)) ! dT*_ij/dFe_kl if (pressure < 0.0_pReal) then do i=1_pInt, 3_pInt; do k=1_pInt,3_pInt; do l=1_pInt,3_pInt; do j=1_pInt,3_pInt