From b59b6523b4a539b61134a4cdee344de9304cdbd7 Mon Sep 17 00:00:00 2001 From: Test User Date: Tue, 2 Oct 2018 02:53:01 +0200 Subject: [PATCH 1/4] [skip ci] updated version information after successful test of v2.0.2-638-gd5cf7d21 --- VERSION | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/VERSION b/VERSION index 26c337c96..79ee48b68 100644 --- a/VERSION +++ b/VERSION @@ -1 +1 @@ -v2.0.2-591-ga00d15b8 +v2.0.2-638-gd5cf7d21 From 7ac6a5e739031647c06b6983d396d4559a227e96 Mon Sep 17 00:00:00 2001 From: Test User Date: Tue, 2 Oct 2018 18:14:32 +0200 Subject: [PATCH 2/4] [skip ci] updated version information after successful test of v2.0.2-700-g27bde055 --- VERSION | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/VERSION b/VERSION index 26c337c96..684700ac8 100644 --- a/VERSION +++ b/VERSION @@ -1 +1 @@ -v2.0.2-591-ga00d15b8 +v2.0.2-700-g27bde055 From e74a852ef7eb6034471ffcc800834fb3ae9ff603 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Fri, 5 Oct 2018 06:38:20 +0200 Subject: [PATCH 3/4] new style is noticable slower runtime (polycrystal example) increased from 70s to 90s. Need to figure out how to write readable AND fast code --- src/plastic_phenopowerlaw.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/plastic_phenopowerlaw.f90 b/src/plastic_phenopowerlaw.f90 index f7c723521..76e7dd4b4 100644 --- a/src/plastic_phenopowerlaw.f90 +++ b/src/plastic_phenopowerlaw.f90 @@ -637,7 +637,7 @@ end subroutine plastic_phenopowerlaw_dotState !> @details: Shear rates are calculated only optionally. NOTE: Agains the common convention, the !> result (i.e. intent(out)) variables are the last to have the optional arguments at the end !-------------------------------------------------------------------------------------------------- -subroutine kinetics_slip(prm,stt,of,Mp,gdot_slip_pos,gdot_slip_neg, & +pure subroutine kinetics_slip(prm,stt,of,Mp,gdot_slip_pos,gdot_slip_neg, & dgdot_dtau_slip_pos,dgdot_dtau_slip_neg) use prec, only: & dNeq0 @@ -703,7 +703,7 @@ end subroutine kinetics_slip !> @details: Shear rates are calculated only optionally. NOTE: Agains the common convention, the !> result (i.e. intent(out)) variables are the last to have the optional arguments at the end !-------------------------------------------------------------------------------------------------- -subroutine kinetics_twin(prm,stt,of,Mp,gdot_twin,dgdot_dtau_twin) +pure subroutine kinetics_twin(prm,stt,of,Mp,gdot_twin,dgdot_dtau_twin) use prec, only: & dNeq0 use math, only: & From 47a748831fad8397cb96a4cc38d5602f3ca41509 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Fri, 5 Oct 2018 07:26:13 +0200 Subject: [PATCH 4/4] Cannot write out resolved stress including non schmid Would need to write out tau_pos and tau_neg separately, so limit to resolved stress in the sense of Schmid --- src/plastic_phenopowerlaw.f90 | 10 +--------- 1 file changed, 1 insertion(+), 9 deletions(-) diff --git a/src/plastic_phenopowerlaw.f90 b/src/plastic_phenopowerlaw.f90 index 76e7dd4b4..edba9a71a 100644 --- a/src/plastic_phenopowerlaw.f90 +++ b/src/plastic_phenopowerlaw.f90 @@ -769,8 +769,6 @@ function plastic_phenopowerlaw_postResults(Mp,instance,of) result(postResults) integer(pInt) :: & o,c,i,j - real(pReal) :: & - tau_slip_pos, tau_slip_neg real(pReal), dimension(param(instance)%totalNslip) :: & gdot_slip_pos,gdot_slip_neg @@ -797,13 +795,7 @@ function plastic_phenopowerlaw_postResults(Mp,instance,of) result(postResults) c = c + prm%totalNslip case (resolvedstress_slip_ID) do i = 1_pInt, prm%totalNslip - tau_slip_pos = math_mul33xx33(Mp,prm%Schmid_slip(1:3,1:3,i)) - tau_slip_neg = tau_slip_pos - !do j = 1,size(prm%nonSchmidCoeff) - ! tau_slip_pos = tau_slip_pos + math_mul33xx33(S,prm%nonSchmid_pos(1:3,1:3,j,i)) - ! tau_slip_neg = tau_slip_neg + math_mul33xx33(S,prm%nonSchmid_neg(1:3,1:3,j,i)) - !enddo - postResults(c+i) = 0.5_pReal*(tau_slip_pos+tau_slip_neg) + postResults(c+i) = math_mul33xx33(Mp,prm%Schmid_slip(1:3,1:3,i)) enddo c = c + prm%totalNslip