![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to whetstone.f CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [Power 6/32 Unix Tahoe 4.2BSD] / cci / usr / src / usr.bin / f77 / f77pass1 / tst|
Return to whetstone.f CVS log | Up to [Power 6/32 Unix Tahoe 4.2BSD] / cci / usr / src / usr.bin / f77 / f77pass1 / tst |
1.1 root 1: c whetstone benchmark
2: c single precision variant of program
3: common t,t1,t2,e1(4),j,k,l
4: t=0.499975e00
5: t1=0.50025e00
6: t2=2.0e00
7: c execute program twice, once with loopcount i=10, and again with i=20.
8: i=10
9: c***** start of timed interval *****
10: 10 btime=secnds(0.0)
11: c
12: c******** test 1 ********
13: c
14: isave=i
15: n1=0
16: n2=12*i
17: n3=14*i
18: n4=345*i
19: n5=0
20: n6=210*i
21: n7=32*i
22: n8=899*i
23: n9=616*i
24: n10=0
25: n11=93*i
26: n12=0
27: x1=1.0e0
28: x2=-1.0e0
29: x3=-1.0e0
30: x4=-1.0e0
31: if(n1)19,19,11
32: 11 do 18 i=1,n1,1
33: x1=(x1+x2+x3-x4)*t
34: x2=(x1+x2-x3+x4)*t
35: x4=(-x1+x2+x3+x4)*t
36: x3=(x1-x2+x3+x4)*t
37: 18 continue
38: 19 continue
39: call pout(n1,n1,n1,x1,x2,x3,x4)
40: c
41: c******** test 2 ********
42: c
43: e1(2)=-1.0e0
44: e1(3)=-1.0e0
45: e1(4)=-1.0e0
46: if(n2)29,29,21
47: 21 do 28 i=1,n2,1
48: e1(1)=(e1(1)+e1(2)+e1(3)-e1(4))*t
49: e1(2)=(e1(1)+e1(2)-e1(3)+e1(4))*t
50: e1(3)=(e1(1)-e1(2)+e1(3)+e1(4))*t
51: e1(4)=(-e1(1)+e1(2)+e1(3)+e1(4))*t
52: 28 continue
53: 29 continue
54: call pout(n2,n3,n2,e1(1),e1(2),e1(3),e1(4))
55: c
56: c******** test 3 ********
57: c
58: 31 do 38 i=1,n3,1
59: 38 call pa(e1)
60: 39 continue
61: call pout(n3,n2,n2,e1(1),e1(2),e1(3),e1(4))
62: c
63: c******** test 4 ********
64: c
65: if(n4)49,49,41
66: 41 do 48 i=1,n4,1
67: if (j-1)43,42,43
68: 42 j=2
69: goto44
70: 43 j=3
71: 44 if(j-2)46,46,45
72: 45 j=0
73: goto47
74: 46 j=1
75: 47 if(j-1)411,412,412
76: 411 j=1
77: goto48
78: 412 j=0
79: 48 continue
80: 49 continue
81: call pout(n4,j,j,x1,x2,x3,x4)
82: c
83: c******** test 5 ********
84: c
85: j=1
86: k=2
87: l=3
88: if(n6)69,69,61
89: 61 do 68 i=1,n6,1
90: j=j*(k-j)*(l-k)
91: k=l*k-(l-j)*k
92: l=(l-k)*(k+j)
93: e1(l-1)=j+k+l
94: e1(k-1)=j+k+l
95: 68 continue
96: 69 continue
97: call pout(n6,j,k,e1(1),e1(2),e1(3),e1(4))
98: c
99: c******** test 6 ********
100: c
101: x=0.5e0
102: y=0.5e0
103: if(n7)79,79,71
104: 71 do 78 i=1,n7,1
105: x=t*atan(t2*sin(x)+cos(x)/cos(x+y)+cos(x-y)-1.0e0)
106: y=t*atan(t2*sin(y)+cos(y)/cos(x+y)+cos(x-y)-1.0e0)
107: 78 continue
108: 79 continue
109: call pout(n7,j,k,x,x,y,y)
110: c
111: c******** test 7 ********
112: c
113: x=1.0e0
114: y=1.0e0
115: z=1.0e0
116: if(n8)89,89,81
117: 81 do 88 i=1,n8,1
118: 88 call p3(x,y,z)
119: 89 continue
120: call pout(n8,j,k,x,y,z,z)
121: c
122: c******** test 8 ********
123: c
124: j=1
125: k=2
126: l=3
127: e1(1)=1.0e0
128: e1(2)=2.0e0
129: e1(3)=3.0e0
130: if(n9)99,99,91
131: 91 do 98 i=1,n9,1
132: 98 call p0
133: 99 continue
134: call pout(n9,j,k,e1(1),e1(2),e1(3),e1(4))
135: c
136: c******** test 9 ********
137: c
138: j=2
139: k=3
140: if(n10)109,109,101
141: 101 do 108 i=1,n10,1
142: j=j+k
143: k=j+k
144: j=j-k
145: k=k-j-j
146: 108 continue
147: 109 continue
148: call pout(n10,j,k,x1,x2,x3,x4)
149: c
150: c******** test 10 ********
151: c
152: x=0.75e0
153: if(n11)119,119,111
154: 111 do 118 i=1,n11,1
155: 118 x=sqrt(exp(alog(x)/t1))
156: 119 continue
157: c
158: c******** end of timed interval ********
159: c save time for loopcount i=10 and repeat for i=20
160: delta=secnds(btime)-btime
161: if(isave.eq.20) go to 120
162: time1=delta
163: i=20
164: goto 10
165: 120 time2=delta
166: c calculate whetstone kips
167: wkip=(1.0/(time2-time1))*1000
168: write(6, 121) wkip
169: 121 format (f10.2,' whetstone kips per second.')
170: stop
171: end
172:
173: subroutine pa(e)
174: common t,t1,t2
175: dimension e(4)
176: j=0
177: 1 e(1)=(e(1)+e(2)+e(3)-e(4))*t
178: e(2)=(e(1)+e(2)-e(3)+e(4))*t
179: e(3)=(e(1)-e(2)+e(3)+e(4))*t
180: e(4)=(-e(1)+e(2)+e(3)+e(4))/t2
181: j=j+1
182: if(j-6)1,2,2
183: 2 continue
184: return
185: end
186: subroutine p0
187: common t,t1,t2,e1(4),j,k,l
188: e1(j)=e1(k)
189: e1(k)=e1(l)
190: e1(l)=e1(j)
191: return
192: end
193: subroutine p3(x,y,z)
194: common t,t1,t2
195: x1=x
196: y1=y
197: x1=t*(x1+y1)
198: y1=t*(x1+y1)
199: z=(x1+y1)/t2
200: return
201: end
202: subroutine pout(n,j,k,x1,x2,x3,x4)
203: write (6,1) n,j,k,x1,x2,x3,x4
204: 1 format (1h ,3i7,4e12.4)
205: return
206: end
207:
208: subroutine prtres(wk)
209: return
210: end
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.