PGI does not no "do concurrent" (yet)
This commit is contained in:
parent
388e661433
commit
714c9b1ecb
|
@ -478,7 +478,7 @@ subroutine constitutive_LpAndItsTangents(Lp, dLp_dS, dLp_dFi, &
|
||||||
|
|
||||||
end select plasticityType
|
end select plasticityType
|
||||||
|
|
||||||
#ifdef __INTEL_COMPILER
|
#if defined(__INTEL_COMPILER) || defined(__PGI)
|
||||||
forall(i = 1_pInt:3_pInt, j = 1_pInt:3_pInt)
|
forall(i = 1_pInt:3_pInt, j = 1_pInt:3_pInt)
|
||||||
#else
|
#else
|
||||||
do concurrent(i = 1_pInt:3_pInt, j = 1_pInt:3_pInt)
|
do concurrent(i = 1_pInt:3_pInt, j = 1_pInt:3_pInt)
|
||||||
|
@ -486,7 +486,7 @@ subroutine constitutive_LpAndItsTangents(Lp, dLp_dS, dLp_dFi, &
|
||||||
dLp_dFi(i,j,1:3,1:3) = math_mul33x33(math_mul33x33(Fi,S),transpose(dLp_dMp(i,j,1:3,1:3))) + &
|
dLp_dFi(i,j,1:3,1:3) = math_mul33x33(math_mul33x33(Fi,S),transpose(dLp_dMp(i,j,1:3,1:3))) + &
|
||||||
math_mul33x33(math_mul33x33(Fi,dLp_dMp(i,j,1:3,1:3)),S)
|
math_mul33x33(math_mul33x33(Fi,dLp_dMp(i,j,1:3,1:3)),S)
|
||||||
dLp_dS(i,j,1:3,1:3) = math_mul33x33(math_mul33x33(transpose(Fi),Fi),dLp_dMp(i,j,1:3,1:3)) ! ToDo: @PS: why not: dLp_dMp:(FiT Fi)
|
dLp_dS(i,j,1:3,1:3) = math_mul33x33(math_mul33x33(transpose(Fi),Fi),dLp_dMp(i,j,1:3,1:3)) ! ToDo: @PS: why not: dLp_dMp:(FiT Fi)
|
||||||
#ifdef __INTEL_COMPILER
|
#if defined(__INTEL_COMPILER) || defined(__PGI)
|
||||||
end forall
|
end forall
|
||||||
#else
|
#else
|
||||||
enddo
|
enddo
|
||||||
|
@ -1167,4 +1167,5 @@ subroutine constitutive_results()
|
||||||
|
|
||||||
end subroutine constitutive_results
|
end subroutine constitutive_results
|
||||||
|
|
||||||
|
|
||||||
end module constitutive
|
end module constitutive
|
||||||
|
|
Loading…
Reference in New Issue