home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / adaptor.zip / adapt.zip / adaptor / examples / simple / prime.f < prev    next >
Text File  |  1993-06-17  |  658b  |  29 lines
  1. C     data parallel source program
  2.  
  3.       program prime
  4.       integer n, s, k
  5.       logical*1 a(:)
  6. !hpf$ distribute a(block)
  7.       print *, 'Input n for counting primes in range 2 to n : '
  8.       read *, n
  9.       allocate (a(1:n))
  10.       call cm_timer_clear (1)
  11.       call cm_timer_start (1)
  12.       a = .true.
  13.       a(1) = .false.
  14.       k = 2
  15.       do while (k*k .le. n)
  16.          a(k*k:n:k) = .false.
  17.          k = k + 1
  18.          do while (.not. a(k))
  19.            k = k + 1
  20.          end do
  21.       end do
  22.       s = count (a)
  23.       call cm_timer_stop (1)
  24.       print *, 'There are ',s,' primes until ', n
  25.       call cm_timer_print (1)
  26.       deallocate (a)
  27.       end
  28.  
  29.