From f7dbda31e8f4c13079e8fdcac68bcc67a33abb96 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Wed, 10 Nov 2021 06:51:57 +0100 Subject: [PATCH] style adjustments - two empty lines after variables - enddo, endif should have space --- src/phase_mechanical_plastic.f90 | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/src/phase_mechanical_plastic.f90 b/src/phase_mechanical_plastic.f90 index fc69a57da..671af5009 100644 --- a/src/phase_mechanical_plastic.f90 +++ b/src/phase_mechanical_plastic.f90 @@ -261,6 +261,7 @@ module subroutine plastic_LpAndItsTangents(Lp, dLp_dS, dLp_dFi, & integer :: & i, j + if (phase_plasticity(ph) == PLASTICITY_NONE_ID) then Lp = 0.0_pReal dLp_dFi = 0.0_pReal @@ -295,7 +296,7 @@ module subroutine plastic_LpAndItsTangents(Lp, dLp_dS, dLp_dFi, & dLp_dFi(i,j,1:3,1:3) = matmul(matmul(Fi,S),transpose(dLp_dMp(i,j,1:3,1:3))) + & matmul(matmul(Fi,dLp_dMp(i,j,1:3,1:3)),S) dLp_dS(i,j,1:3,1:3) = matmul(matmul(transpose(Fi),Fi),dLp_dMp(i,j,1:3,1:3)) ! ToDo: @PS: why not: dLp_dMp:(FiT Fi) - enddo; enddo + end do; end do end if @@ -319,6 +320,7 @@ module function plastic_dotState(subdt,co,ip,el,ph,en) result(broken) Mp logical :: broken + if (phase_plasticity(ph) /= PLASTICITY_NONE_ID) then Mp = matmul(matmul(transpose(phase_mechanical_Fi(ph)%data(1:3,1:3,en)),& phase_mechanical_Fi(ph)%data(1:3,1:3,en)),phase_mechanical_S(ph)%data(1:3,1:3,en))