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

bugfixes 

flaws introduced in 42-xx branch and detected by new test
- outputsize was wrong
- nconstituents is product along the 3 directions, not the sum
This commit is contained in:
Martin Diehl 2018-11-03 15:13:11 +01:00
parent c5b98e544a
commit a678e9b94f
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/homogenization_RGC.f90
View File

@ -197,14 +197,13 @@ subroutine homogenization_RGC_init(fileUnit)
outputID = magnitudemismatch_ID outputID = magnitudemismatch_ID
outputSize = 3_pInt outputSize = 3_pInt
case default case default
end select
if (outputID /= undefined_ID) then if (outputID /= undefined_ID) then
homogenization_RGC_output(i,instance) = outputs(i) homogenization_RGC_output(i,instance) = outputs(i)
homogenization_RGC_sizePostResult(i,instance) = outputSize homogenization_RGC_sizePostResult(i,instance) = outputSize
prm%outputID = [prm%outputID , outputID] prm%outputID = [prm%outputID , outputID]
endif endif
end select
enddo enddo
!-------------------------------------------------------------------------------------------------- !--------------------------------------------------------------------------------------------------
! * assigning cluster orientations ! * assigning cluster orientations
elementLooping: do e = 1_pInt,mesh_NcpElems elementLooping: do e = 1_pInt,mesh_NcpElems
@ -243,6 +242,7 @@ subroutine homogenization_RGC_init(fileUnit)
! (6) Volume discrepancy, (7) Avg relaxation rate component, (8) Max relaxation rate component ! (6) Volume discrepancy, (7) Avg relaxation rate component, (8) Max relaxation rate component
homogState(h)%sizeState = sizeHState homogState(h)%sizeState = sizeHState
homogenization_RGC_sizePostResults(instance) = sum(homogenization_RGC_sizePostResult(:,instance))
homogState(h)%sizePostResults = homogenization_RGC_sizePostResults(instance) homogState(h)%sizePostResults = homogenization_RGC_sizePostResults(instance)
allocate(homogState(h)%state0 (sizeHState,NofMyHomog), source=0.0_pReal) allocate(homogState(h)%state0 (sizeHState,NofMyHomog), source=0.0_pReal)
allocate(homogState(h)%subState0(sizeHState,NofMyHomog), source=0.0_pReal) allocate(homogState(h)%subState0(sizeHState,NofMyHomog), source=0.0_pReal)
@ -1456,7 +1456,7 @@ subroutine grainDeformation(F, avgF, ip, el)
! compute the deformation gradient of individual grains due to relaxations ! compute the deformation gradient of individual grains due to relaxations
instance = homogenization_typeInstance(mesh_homogenizationAt(el)) instance = homogenization_typeInstance(mesh_homogenizationAt(el))
F = 0.0_pReal F = 0.0_pReal
do iGrain = 1_pInt,sum(param(instance)%Nconstituents) do iGrain = 1_pInt,product(param(instance)%Nconstituents)
iGrain3 = grain1to3(iGrain,instance) iGrain3 = grain1to3(iGrain,instance)
do iFace = 1_pInt,6_pInt do iFace = 1_pInt,6_pInt
intFace = getInterface(iFace,iGrain3) intFace = getInterface(iFace,iGrain3)