1! PR middle-end/49721 2! { dg-do compile } 3! { dg-options "-O3 -funroll-loops -std=legacy" } 4 5 subroutine midbloc6(c,a2,a2i,q) 6 parameter (ndim2=6) 7 parameter (ndim=3) 8 dimension ri(ndim2),cr(ndim2,ndim2),xj(ndim2,ndim2),q(*) 9 @,sai(ndim2,ndim2),cm(ndim2,ndim2),w(ndim2,ndim2) 10 dimension vr(ndim2,ndim2),vi(ndim2,ndim2),s1(ndim2,ndim2),p(ndim) 11 dimension xq(6),qb(2),qc(2),ifl(6),iplane(3) 12 save 13 call eig66(cr,rr,ri,vr,vi) 14 xq(i)=asin(ri(i))/x2pi 15 i9=6 16 qb(1)=q(1)/x2pi 17 do 180 i=1,2 18 do 170 j=1,6 19 120 if(xq(j)) 130,190,140 20 130 if(qb(i)-0.5d0) 160,150,150 21 140 if(qb(i)-0.5d0) 150,150,160 22 150 continue 23 tst=abs(abs(qb(i))-abs(xq(j))) 24 160 continue 25 170 continue 26 iplane(i)=k 27 180 continue 28 190 continue 29 n1=iplane(3) 30 if(i9.eq.6) then 31 z=vr(1,n1)*vi(2,n1)-vr(2,n1)*vi(1,n1)+vr(3,n1)*vi(4,n1)-vr(4,n1) 32 endif 33 sai(6,i)=vi(i,n1)/z 34 call dacond6(a2,zero) 35 end 36