From bdaa703d37e25618faceb1e782f1bd073da07ef6 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Wed, 3 Apr 2019 16:13:30 +0200 Subject: [PATCH] forall is deprecated, do concurrent not yet established --- src/constitutive.f90 | 14 +++----------- 1 file changed, 3 insertions(+), 11 deletions(-) diff --git a/src/constitutive.f90 b/src/constitutive.f90 index 4df97b240..23ae3f07b 100644 --- a/src/constitutive.f90 +++ b/src/constitutive.f90 @@ -476,19 +476,11 @@ subroutine constitutive_LpAndItsTangents(Lp, dLp_dS, dLp_dFi, & end select plasticityType -#if defined(__INTEL_COMPILER) || defined(__PGI) - forall(i = 1_pInt:3_pInt, j = 1_pInt:3_pInt) -#else - do concurrent(i = 1_pInt:3_pInt, j = 1_pInt:3_pInt) -#endif + do i=1,3; do j=1,3 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) -#if defined(__INTEL_COMPILER) || defined(__PGI) - end forall -#else - enddo -#endif + 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 subroutine constitutive_LpAndItsTangents