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