p25.f000064400003130000310000000023650617341261000120370ustar00krivecgrp000000000000000 PROGRAM P25 C C EACH STEP TIMED SEPARATELY. C ONLY ONE LOOP CONFIGURATION TESTED (BECAUSE RRPSK 1.0 PROVED C STABILITY ON LOOP REARRANGEMENTS). C I-LOOP INSIDE (FOR COMPILERS WHICH DO NOT EXCHANGE LOOPS). C V1, 92/11/26. C IMPLICIT REAL*8 (A-H,O-Z) PARAMETER (ND = 400, NIN = 5, NOUT = 6) DIMENSION A(ND,ND), B(ND,ND), C(ND,ND) C WRITE (NOUT,200) READ (NIN,100) N WRITE (NOUT,201) ND, N CALL TEMPD(T0, T1, TD, NOUT, .TRUE.) DO 12 J = 1,N DO 10 I = 1,N A(I,J) = (J * I) B(I,J) = (J + I) 10 CONTINUE 12 CONTINUE CALL TEMPD(T1, T2, TD, NOUT, .TRUE.) DO 24 J = 1,N DO 18 I = 1,N C(I,J) = 0.D0 18 CONTINUE 24 CONTINUE CALL TEMPD(T2, T3, TD, NOUT, .TRUE.) DO 26 J = 1,N DO 22 K = 1,N DO 20 I = 1,N C(I,J) = C(I,J) + A(I,K) * B(K,J) 20 CONTINUE 22 CONTINUE 26 CONTINUE CALL TEMPD(T3, T4, TD, NOUT, .TRUE.) WRITE (NOUT,202) N, ND, ((C(I,J), J = 1,4), I = 1,4) C 100 FORMAT (I4) 200 FORMAT (1H , 8HP25 V1 ) 201 FORMAT (1H , 8HND, N , 2I8) 202 FORMAT (1H , 2I8, /, (1H , 4E16.8)) END = 0.D0 18 CONTINUE 24 CONTINUE CALL TEMPD(T2, T3, TD, NOUT, .TRUE.) DO 26 J = 1,N DO 22 K = 1,N DO 20 I = 1,N C(I,J) = C(I,J) + A(I,K) * B(K,J) 20 CONTINUE 22 CONTINUE 26 CONTIp25v.f000064400003130000310000000022220617341264000122200ustar00krivecgrp000000000000000 PROGRAM P25 C C EACH STEP TIMED SEPARATELY. C ONLY ONE LOOP CONFIGURATION TESTED (BECAUSE RRPSK 1.0 PROVED C STABILITY ON LOOP REARRANGEMENTS). C I-LOOP INSIDE (FOR COMPILERS WHICH DO NOT EXCHANGE LOOPS). C V1, 92/11/26. C IMPLICIT REAL*8 (A-H,O-Z) PARAMETER (ND = 4000, NIN = 5, NOUT = 6) DIMENSION A(ND,ND), B(ND,ND), C(ND,ND) C WRITE (NOUT,200) READ (NIN,100) N WRITE (NOUT,201) ND, N CALL TEMPD(T0, T1, TD, NOUT, .TRUE.) DO 12 J = 1,N DO 10 I = 1,N A(I,J) = (J * I) B(I,J) = (J + I) 10 CONTINUE 12 CONTINUE CALL TEMPD(T1, T2, TD, NOUT, .TRUE.) DO 24 J = 1,N DO 18 I = 1,N C(I,J) = 0.D0 18 CONTINUE 24 CONTINUE CALL TEMPD(T2, T3, TD, NOUT, .TRUE.) CALL DGEMM('N', 'N', N, N, N, 1.0D0, & A, ND, B, ND, 0.D0, C, ND) CALL TEMPD(T3, T4, TD, NOUT, .TRUE.) WRITE (NOUT,202) N, ND, ((C(I,J), J = 1,4), I = 1,4) C 100 FORMAT (I4) 200 FORMAT (1H , 8HP25 V1 ) 201 FORMAT (1H , 8HND, N , 2I8) 202 FORMAT (1H , 2I8, /, (1H , 4E16.8)) END EMPD(T1, T2, TD, NOUT, .TRUE.) DO 24 J = 1,N DO 18 I = 1,N C(I,J) = 0.D0 18 CONTINUE 24 CONTINUE CALL TEMPD(T2, T3, TD, NOUT, .TRUE.) CALL DGEMM('N', 'N', N, N, N, 1.0D0, & A, ND, B, ND, 0.D0, C, ND) CALL TEMPD(T3, T4, TD, NOUT, .TRUE.) WRITE (NOUT,202) N, ND, ((C(I,J), J = 1,4), I tempd_hp.f000064400003130000310000000027040617341265700132410ustar00krivecgrp000000000000000 SUBROUTINE TEMPD(TOLD, TNEW, TDIF, NOUT, LPR) INTEGER NOUT REAL*8 TOLD, TNEW, TDIF LOGICAL LPR C C SGI PCL TIMING ROUTINE, USING ETIME. C ETIME REPORTS ELAPSED EXECUTION TIME (USER, SYSTEM) SINCE THE C LAST CALL TO ITSELF; THEREFORE THIS PROGRAM WORKS OK ONLY IF THERE C ARE NO OTHER CALLS TO ETIME EXCEPT FROM THIS PROGRAM. C (ETIME MAY NOT MEASURE THE CPU TIME ITSELF.) C C INPUT: TOLD, NOUT, LPR (ALL UNCHANGED). C OUTPUT: TNEW, TDIF. C LOGICAL INIT, TEST REAL*4 ETIME, TARRAY DIMENSION TARRAY(2) REAL*8 TSUM, SDIF DATA INIT / .TRUE. /, TEST / .FALSE. / C IF (INIT) THEN WRITE (NOUT,2000) TSUM = ETIME(TARRAY) TDIF = TARRAY(1) TNEW = TOLD + TDIF SDIF = TARRAY(2) INIT = .FALSE. IF (TEST) WRITE (NOUT,2200) TSUM, TNEW, TDIF, TOLD, SDIF RETURN ENDIF TSUM = ETIME(TARRAY) TDIF = TARRAY(1) TNEW = TOLD + TDIF SDIF = TARRAY(2) IF (TEST) WRITE (NOUT,2200) TSUM, TNEW, TDIF, TOLD, SDIF IF (LPR) WRITE (NOUT,2100) TNEW, TDIF, TOLD, SDIF RETURN C 2000 FORMAT (' -TEMPD- 3.0, 95/12/20. SGI PC TIMING ROUTINE.') 2100 FORMAT (' TIME', F12.2, A ' DIF', F12.2, B ' REF', F12.2, C ' DIF SYS', F12.2) 2200 FORMAT (' SUM USER DIF OLD SYS', 5F8.2) END TSUM = ETIME(TARRAY) TDIF = TARRAY(1) TNEW =tempd_sgi.f000066600003130000310000000032510617341265700134160ustar00krivecgrp000000000000000 SUBROUTINE TEMPD(TOLD, TNEW, TDIF, NOUT, LPR) INTEGER NOUT REAL*8 TOLD, TNEW, TDIF LOGICAL LPR C C SGI PCL TIMING ROUTINE, USING DTIME. C DTIME REPORTS ELAPSED EXECUTION TIME (USER, SYSTEM) SINCE THE C LAST CALL TO ITSELF; THEREFORE THIS PROGRAM WORKS OK ONLY IF THERE C ARE NO OTHER CALLS TO DTIME EXCEPT FROM THIS PROGRAM. C (DTIME MAY NOT MEASURE THE CPU TIME ITSELF.) C C INPUT: TOLD, NOUT, LPR (ALL UNCHANGED). C OUTPUT: TNEW (CUMULATIVE TIME), TDIF (DIFFERENCE FROM TOLD). C TCUM SAVES THE CUMULATIVE TIME. C LOGICAL INIT, TEST REAL*4 DTIME, TARRAY DIMENSION TARRAY(2) REAL*8 TSUM, TCUM, SCUM, SDIF DATA INIT / .TRUE. /, TEST / .FALSE. / DATA TCUM / 0.D0 /, SCUM / 0.D0 / C IF (INIT) THEN WRITE (NOUT,2000) INIT = .FALSE. C TSUM = DTIME(TARRAY) TCUM = TCUM + TARRAY(1) SCUM = SCUM + TARRAY(2) C TNEW = TCUM TDIF = TNEW - TOLD IF (TEST) WRITE (NOUT,2200) TSUM, TNEW, TDIF, TOLD, SCUM, SDIF RETURN ENDIF C TSUM = DTIME(TARRAY) TCUM = TCUM + TARRAY(1) SCUM = SCUM + TARRAY(2) C TNEW = TCUM TDIF = TNEW - TOLD IF (TEST) WRITE (NOUT,2200) TSUM, TNEW, TDIF, TOLD, SCUM, SDIF IF (LPR) WRITE (NOUT,2100) TNEW, TDIF, TOLD, SCUM RETURN C 2000 FORMAT (' -TEMPD- 4.0, 96/01/19. SGI PC TIMING ROUTINE.') 2100 FORMAT (' USER', F12.2, A ' DIF', F12.2, B ' REF', F12.2, C ' SYS', F12.2) 2200 FORMAT (' SUM USER DIF OLD SYS SDIF', 6F8.2) END TARRAY(1) SCUM = SCUM + TARRAY(2) C TNEW = TCUM TDIF = TNEW - TOLD IF (TEST) WRITE (NOUT,2200) TSUM, TNEW, TDIF, TOLD, SCUM, SDIF IF (LPR) WRITE (NOUT,2100) TNEW, TDIF, TOLD, SCUM RETURN C 2000 FORMAT (' -TEMPD- 4.0, 96/01/19. SGI PC TIMING ROUTINE.') 2100 FORMAT (' USER', F12.2, A