3 ! { dg-options "-O3 -funroll-loops" }
5 subroutine midbloc6(c,a2,a2i,q)
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)
13 call eig66(cr,rr,ri,vr,vi)
14 xq(i)=asin(ri(i))/x2pi
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
23 tst=abs(abs(qb(i))-abs(xq(j)))
31 z=vr(1,n1)*vi(2,n1)-vr(2,n1)*vi(1,n1)+vr(3,n1)*vi(4,n1)-vr(4,n1)