C WHETSTONE BENCHMARK C C THIS FORTRAN PROGRAM IS A TRANSLATION FROM C THE ALGORITHM DESCRIBED IN "A SYNTHETIC BENCHMARK" C H.J. CURNOW AND B.A. WICHMANN , VOL 19, NUMBER 1 C THE COMPUTER JOURNAL. C COMMON T,T1,T2,E1(4),J,K,L C C INITIALIZE SERIES APPROXIMATION CONSTANTS AND GET LOOP COUNTER C T=0.499975E00 T1=0.50025E00 TYPE 998 998 FORMAT(1X,'ITERATIONS ?') T2=2.0E00 ACCEPT 999, I 999 FORMAT(I4) N1=0 N2=12*I N3=14*I N4=345*I N5=0 N6=210*I N7=32*I N8=899*I N9= 616*I N10=0 N11=93*I N12=0 C INITIALIZE LOOP ITERATIONS C C********** MODULE 1: SIMPLE IDENTIFIERS ****************** C X1=1.0E0 X2=-1.0E0 X3=-1.0E0 X4=-1.0E0 IF(N1)19,19,11 11 DO 18 I=1,N1,1 X1=(X1+X2+X3-X4)*T X2=(X1+X2-X3+X4)*T X3=(X1-X2+X3+X4)*T X4=(-X1+X2+X3+X4)*T 18 CONTINUE 19 CONTINUE C C******** MODULE 2: ARRAY ELEMENTS ************* C E1(1)=1.0E0 E1(2)=-1.0E0 E1(3)=-1.0E0 E1(4)=-1.0E0 IF(N2)29,29,21 21 DO 28 I=1,N2,1 E1(1)=(E1(1)+E1(2)+E1(3)-E1(4))*T E1(2)=(E1(1)+E1(2)-E1(3)+E1(4))*T E1(3)=(E1(1)-E1(2)+E1(3)+E1(4))*T E1(4)=(-E1(1)+E1(2)+E1(3)+E1(4))*T 28 CONTINUE 29 CONTINUE C C********** MODULE 3: ARRAY AS PARAMETER ******* C 31 DO 38 I=1,N3,1 38 CALL PA(E1) 39 CONTINUE C C********* MODULE 4: CONDITIONAL JUMPS *********** C j = 1 IF(N4)49,49,41 41 DO 48 I=1,N4,1 IF(J-1)43,42,43 42 J=2 GOTO 44 43 J=3 44 IF(J-2)46,46,45 45 J=0 GOTO 47 46 J=1 47 IF(J-1)411,412,412 411 J=1 GOTO 48 412 J=0 48 CONTINUE 49 CONTINUE C C************ MODULE 5; OMITTED ******** C************ MODULE 6: INTEGER ARITHMETIC ********** J=1 K=2 L=3 IF(N6)69,69,61 61 DO 68 I=1,N6,1 J=J*(K-J)*(L-K) K=L*K-(L-J)*K L=(L-K)*(K+J) E1(L-1)=J+K+L E1(K-1)=J*K*L 68 CONTINUE 69 CONTINUE C C********** MODULE 7: TRIG FUNCTIONS ********** C X=0.5E0 Y=0.5E0 IF(N7)79,79,71 71 DO 78 I=1,N7,1 X=T*ATAN(T2*SIN(X)*COS(X)/(COS(X+Y)+COS(X-Y)-1.0E0)) Y=T*ATAN(T2*SIN(Y)*COS(Y)/(COS(X+Y)+COS(X-Y)-1.0E0)) 78 CONTINUE 79 CONTINUE C C*********** MODULE 8: PROCEDURE CALLS ********** C X=1.0E0 Y=1.0E0 Z=1.0E0 IF(N8)89,89,81 81 DO 88 I=1,N8,1 88 CALL P3(X,Y,Z) 89 CONTINUE C C************ MODULE 9: ARRAY REFERENCES ********** C J=1 K=2 L=3 E1(1)=1.0E0 E1(2)=2.0E0 E1(3)=3.0E0 IF(N9)99,99,91 91 DO 98 I=1,N9,1 98 CALL P0 99 CONTINUE C C********* MODULE 10: INTEGER ARITHMETIC *********** C J=2 K=3 IF(N10)109,109,101 101 DO 108 I=1,N10,1 J=J+K K=J+K J=K-J K=K-J-J 108 CONTINUE 109 CONTINUE C C********** MODULE 11: STANDARD FUNCTIONS ********** C X=0.75E0 IF(N11)119,119,111 111 DO 118 I=1,N11,1 118 X=SQRT(EXP(ALOG(X)/T1)) 119 CONTINUE STOP END C C*************** END OF TIMED INTERVAL *********** C SUBROUTINE PA(E) COMMON T,T1,T2 DIMENSION E(4) J=0 1 E(1)=(E(1)+E(2)+E(3)-E(4))*T E(2)=(E(1)+E(2)-E(3)+E(4))*T E(3)=(E(1)-E(2)+E(3)+E(4))*T E(4)=(-E(1)+E(2)+E(3)+E(4))/T2 J=J+1 IF(J-6)1,2,2 2 CONTINUE RETURN END C SUBROUTINE P0 COMMON T,T1,T2,E1(4),J,K,L E1(J)=E1(K) E1(K)=E1(L) E1(L)=E1(J) RETURN END C SUBROUTINE P3(X,Y,Z) COMMON T,T1,T2 X1=X Y1=Y X1=T*(X1+Y1) Y1=T*(X1+Y1) Z=(X1+Y1)/T2 RETURN END C