Example —Loop Nest

subroutine fem3d(…)

10 call addmon(…)

if(numelh.ne0) call solide

subroutine solide

do 20 i=1,nelt

do 20 j=1,nelg

call unpki

call strain

call force

20 continue

if(…) return

goto 10

subroutine force(…)

do 10 i=lft,llt

sgv(i) = sig1(i)-qp(i)*vol(i)

10 continue

do 50 n=1,nnc

i0=ia(n)

i1=ia(n+1)-1

do 50 i=i0,i1

e(1,ix(i))=e(1,ix(i))+ep11(i)

50 continue

Previous slide

Next slide

Back to first slide

View graphic version