From 824cec5ca1c688c64d2be92cef724610729eca36 Mon Sep 17 00:00:00 2001 From: Franz Roters Date: Tue, 11 Nov 2008 14:57:20 +0000 Subject: [PATCH] correct error in case of outdated FFN1: element number is cp_en not e --- trunk/CPFEM_Taylor.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/trunk/CPFEM_Taylor.f90 b/trunk/CPFEM_Taylor.f90 index b87fb30ca..1af6a5326 100644 --- a/trunk/CPFEM_Taylor.f90 +++ b/trunk/CPFEM_Taylor.f90 @@ -156,7 +156,7 @@ if (any(abs(ffn1 - CPFEM_ffn1_bar(:,:,CPFEM_in,cp_en)) > relevantStrain)) then CPFEM_stress_bar(1:CPFEM_ngens,:,:) = CPFEM_odd_stress odd_jaco = CPFEM_odd_jacobian*math_identity2nd(CPFEM_ngens) - forall (i=1:mesh_NcpElems, j=1:FE_Nips(mesh_element(2,e))) & + forall (i=1:mesh_NcpElems, j=1:FE_Nips(mesh_element(2,cp_en))) & CPFEM_jaco_bar(1:CPFEM_ngens,1:CPFEM_ngens,j,i) = odd_jaco outdatedFFN1 = .true. if (.not. CPFEM_calc_done .AND.CPFEM_mode == 1) then