From 4f1cb3a8af125e02babf549ad8112b3f17cb8fe6 Mon Sep 17 00:00:00 2001 From: Onur Guevenc Date: Wed, 3 Aug 2011 17:58:16 +0000 Subject: [PATCH] exchanged forall with nested do loops... Bye, bye smatrix warning, will miss you!! --- code/homogenization_RGC.f90 | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/code/homogenization_RGC.f90 b/code/homogenization_RGC.f90 index eb89c1d72..9d31dce00 100644 --- a/code/homogenization_RGC.f90 +++ b/code/homogenization_RGC.f90 @@ -594,8 +594,9 @@ function homogenization_RGC_updateState(& mornN = homogenization_RGC_interfaceNormal(intFaceN,ip,el) ! get normal of the interfaces iMun = homogenization_RGC_interface4to1(intFaceN,homID) ! translate the interfaces ID into local 4-dimensional index if (iMun .gt. 0) then ! get the corresponding tangent - forall(i=1:3,j=1:3,k=1:3,l=1:3) & - smatrix(3*(iNum-1)+i,3*(iMun-1)+j) = smatrix(3*(iNum-1)+i,3*(iMun-1)+j) + dPdF(i,k,j,l,iGrN)*normN(k)*mornN(l) + do i=1,3; do j=1,3; do k=1,3; do l=1,3 + smatrix(3*(iNum-1)+i,3*(iMun-1)+j) = smatrix(3*(iNum-1)+i,3*(iMun-1)+j) + dPdF(i,k,j,l,iGrN)*normN(k)*mornN(l) + enddo;enddo;enddo;enddo ! projecting the material tangent dPdF into the interface ! to obtain the Jacobian matrix contribution of dPdF endif @@ -612,8 +613,9 @@ function homogenization_RGC_updateState(& mornP = homogenization_RGC_interfaceNormal(intFaceP,ip,el) ! get normal of the interfaces iMun = homogenization_RGC_interface4to1(intFaceP,homID) ! translate the interfaces ID into local 4-dimensional index if (iMun .gt. 0) then ! get the corresponding tangent - forall(i=1:3,j=1:3,k=1:3,l=1:3) & - smatrix(3*(iNum-1)+i,3*(iMun-1)+j) = smatrix(3*(iNum-1)+i,3*(iMun-1)+j) + dPdF(i,k,j,l,iGrP)*normP(k)*mornP(l) + do i=1,3; do j=1,3; do k=1,3; do l=1,3 + smatrix(3*(iNum-1)+i,3*(iMun-1)+j) = smatrix(3*(iNum-1)+i,3*(iMun-1)+j) + dPdF(i,k,j,l,iGrP)*normP(k)*mornP(l) + enddo;enddo;enddo;enddo endif enddo enddo