! { dg-do compile } ! { dg-options "-O3" } implicit real*8 (a-h,o-z) common clop6(3),dps(6),aml6(6,6) dimension y1(3) dimension dclo(3) dimension dx(3),dy(3) save do 80 ii=1,itco y1(3)=dps(1) do 40 l=1,3 dy(l)=clop6(l)-y1(l) 40 continue dczp=abs(dy(3)) if(dcx.le.c1m10.and.dcz.le.c1m10.and.dcxp.le.c1m10.and.dczp + .le.c1m10.and.dcy.le.c1m10.and.dcyp.le.c1m10) goto 90 80 continue write(6) itco ii=itco 90 continue if(ii.ne.itco) then do 65 k=1,3 do 55 j=1,3 jj=2*j kk=2*k dclo(k)=aml6(kk-1,jj-1)*dx(j)+dclo(k) dclo(k)=aml6(kk-1,jj)*dy(j)+dclo(k) 55 continue 65 continue endif end