blob: 163a276be29abff9489efed96500c7ed977241e8 (
plain)
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
|
! { dg-do compile }
! { dg-options "-O3" }
implicit real*8 (a-h,o-z)
common clop6(3),dps(6),aml6(6,6)
dimension y1(3)
dimension dclo(3)
dimension dx(3),dy(3)
save
do 80 ii=1,itco
y1(3)=dps(1)
do 40 l=1,3
dy(l)=clop6(l)-y1(l)
40 continue
dczp=abs(dy(3))
if(dcx.le.c1m10.and.dcz.le.c1m10.and.dcxp.le.c1m10.and.dczp
+ .le.c1m10.and.dcy.le.c1m10.and.dcyp.le.c1m10) goto 90
80 continue
write(6) itco
ii=itco
90 continue
if(ii.ne.itco) then
do 65 k=1,3
do 55 j=1,3
jj=2*j
kk=2*k
dclo(k)=aml6(kk-1,jj-1)*dx(j)+dclo(k)
dclo(k)=aml6(kk-1,jj)*dy(j)+dclo(k)
55 continue
65 continue
endif
end
|
