home *** CD-ROM | disk | FTP | other *** search
/ rtsi.com / 2014.01.www.rtsi.com.tar / www.rtsi.com / OS9 / OSK / EFFO / pd0.lzh / DEMO / prime.f < prev    next >
Text File  |  1987-09-02  |  925b  |  37 lines
  1. options intsize=2, CPU=68020
  2.       program erathostenes
  3.  
  4. c *** benchmark zur berechnung der primzahlen mit dem sieb des
  5. c *** erathostenes
  6.  
  7.       implicit none
  8.       parameter size = 8190
  9.       common flags
  10.       logical*1 flags(0:size)
  11.       integer   i, prime, k, count, iter, is 
  12.       equivalence (i, d7), (k, d6), (count, d5), (prime, d4)
  13.  
  14.       type *,'10 Iterations', char(7)
  15. c      call lib$init_timer
  16.       do iter = 1, 10
  17.         count = 0
  18.         do i = 0, size
  19.           flags(i) = .true.
  20.         end do
  21.         do i = 0, size
  22.           if (flags(i)) then
  23.             prime = i + i + 3
  24.             k     = i + prime
  25.             do while (k .le. size)
  26.               flags(k) = .false.
  27.               k = k + prime
  28.             end do
  29.             count = count + 1
  30.           end if
  31.         end do
  32.       end do
  33. c      call lib$show_timer
  34.       call shell('procs', is)
  35.       type *, char(7), count,' primes.'
  36.       end
  37.