home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Professional / OS2PRO194.ISO / os2 / prgramer / adaptor / examples / purdue / prob02.fcm < prev    next >
Encoding:
Text File  |  1993-06-26  |  1.8 KB  |  81 lines
  1.       PROGRAM PROB02
  2. C
  3. C     PROBLEM 2
  4. C
  5. C  REFERENCE:  PROBLEMS TO TEST PARALLEL AND VECTOR LANGUAGES
  6. C              CSD-TR 516, COMPUTER SCIENCE, PURDUE UNIVERSITY
  7. C              JOHN R. RICE, MAY 1, 1985
  8. C
  9. C              REVISED BY JOHN R. RICE AND J. JING, OCT. 1, 1990
  10. C
  11. C
  12. C      *************************************************
  13. C      *      Adapted for FORTRAN D benchmarking       *
  14. C      *    by  T. HAUPT  (haupt@sccs.npac.syr.edu)    *
  15. C      *                                               *
  16. C      *    Northeast Parallel Architectures Center    *
  17. C      *   at Syracuse University, Syracuse, NY, USA   *
  18. C      *************************************************
  19. C
  20. C
  21. C       VERSION SIMD/CM2-1.00
  22. C       ==================================================
  23. C
  24.       INCLUDE '/usr/include/cm/paris-configuration-fort.h'
  25.       INTEGER KASES,NK,MK
  26.       PARAMETER (KASES=4)
  27.       INTEGER N(KASES),M(KASES)
  28. cmf$  layout m (:serial)
  29. cmf$  layout n (:serial)
  30.       DATA N / 128,1024,1024,1024 /
  31.       DATA M / 128,64,512,1024 /
  32.  
  33.       DO 50 K = 1, KASES
  34.  
  35.       CALL CM_TIMER_CLEAR(0)
  36.       CALL CM_TIMER_START(0)
  37.  
  38.       DO MANY=1,25
  39.  
  40.       NK=N(K)
  41.       MK=M(K)
  42.          CALL DOIT(NK,MK,ESTAR)
  43.  
  44.  
  45.       ENDDO
  46.  
  47.       CALL CM_TIMER_STOP(0)
  48.  
  49.       PRINT 60, N(K),M(K)
  50.   60  FORMAT ('PROBLEM 2 WITH N,M =',I6,2X,I6)
  51.       PRINT *,'GIVES ESTAR =',ESTAR
  52.  
  53.       CALL CM_TIMER_PRINT(0)
  54.  
  55.    50 CONTINUE
  56.       END
  57.  
  58.  
  59.       SUBROUTINE DOIT(NK,MK,ESTAR)
  60.       INTEGER NK,MK
  61.       REAL ESTAR
  62. c     REAL, ARRAY(NK,MK) :: TEMP
  63.       REAL, ARRAY(MK,NK) :: TEMP
  64.       REAL, ARRAY(NK)    :: H
  65.  
  66.       TEMP=1.+0.5/(REAL(-IABS(SPREAD([1:MK],2,NK)
  67.      *                          -SPREAD([1:NK],1,MK)))+0.001)
  68. !HPF$ independent, local_access
  69.       do i = 1, nk 
  70.          h(i) = 1.0
  71.          do j = 1, mk
  72.             h(i) = h(i) * temp(j,i)
  73.          end do
  74.       end do
  75.       ESTAR=SUM(h)
  76.  
  77.       END
  78.  
  79.  
  80.  
  81.