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

disabled a faulty ntrans loop (matmul was complaining)

This commit is contained in:
Martin Diehl 2015-01-23 11:23:35 +00:00
parent 4772efbc45
commit e9c3137a3e
2 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
code/constitutive.f90
View File

@ -1524,8 +1524,8 @@ function constitutive_getSlipDamage(nSlip, Tstar_v, ipc, ip, el)
el !< element number
real(pReal), intent(in), dimension(6) :: &
Tstar_v !< 2nd Piola-Kirchhoff stress
real(pReal) :: &
constitutive_getSlipDamage(nSlip)
real(pReal), dimension(nSlip) :: &
constitutive_getSlipDamage
select case (phase_damage(material_phase(ipc,ip,el)))
case (LOCAL_DAMAGE_isoDuctile_ID)

7
code/plastic_dislotwin.f90
View File

@ -1875,15 +1875,16 @@ subroutine plastic_dislotwin_dotState(Tstar_v,Temperature,nSlipDamage,slipDamage
shear_trans = 0.0_pReal
shearrate_trans = 0.0_pReal
probrate_trans = 0.0_pReal
#ifdef TRANS
j = 0_pInt
do f = 1_pInt,lattice_maxNtransFamily ! loop over all trans families
index_myFamily = sum(lattice_NtransSystem(1:f-1_pInt,ph)) ! at which index starts my family
!* Projection of shear and shear rate on fault band (twin) systems
shear_trans = matmul(plastic_dislotwin_projectionMatrix_Trans(:,:,instance), &
plasticState(ph)%state(2_pInt*ns+1_pInt:3_pInt*ns, of))
plasticState(ph)%state(2_pInt*ns+1_pInt:3_pInt*ns, of)) This is not safe
shearrate_trans = matmul(plastic_dislotwin_projectionMatrix_Trans(:,:,instance), &
plasticState(ph)%dotState(2_pInt*ns+1_pInt:3_pInt*ns, of))
plasticState(ph)%dotState(2_pInt*ns+1_pInt:3_pInt*ns, of)) This is not safe
do i = 1_pInt,plastic_dislotwin_Ntrans(f,instance) ! process each (active) trans system in family
j = j+1_pInt
@ -1933,7 +1934,7 @@ subroutine plastic_dislotwin_dotState(Tstar_v,Temperature,nSlipDamage,slipDamage
enddo
enddo
#endif
end subroutine plastic_dislotwin_dotState