2017-02-20 Paul Thomas <pault@gcc.gnu.org>
[official-gcc.git] / gcc / testsuite / gfortran.dg / pr49721-1.f
blob39e2ed74ef728efb646aa6f56f807c161f97b006
1 ! PR middle-end/49721
2 ! { dg-do compile }
3 ! { dg-options "-O3 -funroll-loops" }
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