Advertisement
Guest User

Untitled

a guest
Dec 8th, 2016
132
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.68 KB | None | 0 0
  1. import numpy as np
  2. from numba import vectorize, int32, int64, bool_
  3. import time as time_
  4.  
  5.  
  6. @vectorize([bool_(int32), bool_(int64)], nopython=True, target='parallel')
  7. def is_prime(n):
  8.     i = 2
  9.     while True:
  10.         if not n % i:
  11.             return False
  12.         i += 1
  13.         if i >= n:
  14.             return True
  15.     return True
  16.  
  17.  
  18. def primes_below_p(n):
  19.     l = []
  20.     primes = is_prime(np.arange(n + 1))
  21.     for j in xrange(9, n + 1):
  22.         if primes[j - 6] and primes[j]:
  23.             l.append((j-6, j))
  24.     return l
  25.  
  26.  
  27. a = int(round(time_.time() * 1000))
  28. primes = primes_below_p(100 * 1000)
  29. b = int(round(time_.time() * 1000))
  30. print(primes)
  31. print("%s mils" % (b-a))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement