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
)