223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

correct error in case of outdated FFN1: element number is cp_en not e

This commit is contained in:
Franz Roters 2008-11-11 14:57:20 +00:00
parent db245e8935
commit 824cec5ca1
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
trunk/CPFEM_Taylor.f90
View File

@ -156,7 +156,7 @@
if (any(abs(ffn1 - CPFEM_ffn1_bar(:,:,CPFEM_in,cp_en)) > relevantStrain)) then if (any(abs(ffn1 - CPFEM_ffn1_bar(:,:,CPFEM_in,cp_en)) > relevantStrain)) then
CPFEM_stress_bar(1:CPFEM_ngens,:,:) = CPFEM_odd_stress CPFEM_stress_bar(1:CPFEM_ngens,:,:) = CPFEM_odd_stress
odd_jaco = CPFEM_odd_jacobian*math_identity2nd(CPFEM_ngens) 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 CPFEM_jaco_bar(1:CPFEM_ngens,1:CPFEM_ngens,j,i) = odd_jaco
outdatedFFN1 = .true. outdatedFFN1 = .true.
if (.not. CPFEM_calc_done .AND.CPFEM_mode == 1) then if (.not. CPFEM_calc_done .AND.CPFEM_mode == 1) then