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

preventing NaN in dotState

This commit is contained in:
Martin Diehl 2018-09-12 11:47:35 +02:00
parent 2337dde525
commit 253f318e55
1 changed files with 2 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/plastic_phenopowerlaw.f90
View File

@ -440,6 +440,7 @@ subroutine plastic_phenopowerlaw_init
startIndex = 1_pInt startIndex = 1_pInt
endIndex = prm%totalNslip endIndex = prm%totalNslip
state (instance)%s_slip => plasticState(p)%state (startIndex:endIndex,:) state (instance)%s_slip => plasticState(p)%state (startIndex:endIndex,:)
state (instance)%s_slip = spread(math_expand(prm%tau0_slip, prm%Nslip), 2, NipcMyPhase)
dotState(instance)%s_slip => plasticState(p)%dotState(startIndex:endIndex,:) dotState(instance)%s_slip => plasticState(p)%dotState(startIndex:endIndex,:)
plasticState(p)%state0(startIndex:endIndex,:) = & plasticState(p)%state0(startIndex:endIndex,:) = &
spread(math_expand(prm%tau0_slip, prm%Nslip), 2, NipcMyPhase) spread(math_expand(prm%tau0_slip, prm%Nslip), 2, NipcMyPhase)
@ -448,6 +449,7 @@ subroutine plastic_phenopowerlaw_init
startIndex = endIndex + 1_pInt startIndex = endIndex + 1_pInt
endIndex = endIndex + prm%totalNtwin endIndex = endIndex + prm%totalNtwin
state (instance)%s_twin => plasticState(p)%state (startIndex:endIndex,:) state (instance)%s_twin => plasticState(p)%state (startIndex:endIndex,:)
state (instance)%s_twin = spread(math_expand(prm%tau0_twin, prm%Ntwin), 2, NipcMyPhase)
dotState(instance)%s_twin => plasticState(p)%dotState(startIndex:endIndex,:) dotState(instance)%s_twin => plasticState(p)%dotState(startIndex:endIndex,:)
plasticState(p)%state0(startIndex:endIndex,:) = & plasticState(p)%state0(startIndex:endIndex,:) = &
spread(math_expand(prm%tau0_twin, prm%Ntwin), 2, NipcMyPhase) spread(math_expand(prm%tau0_twin, prm%Ntwin), 2, NipcMyPhase)