1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
"""
Compare the speed of primes sequentially vs. using futures.
"""
import sys
import time
import math
try:
from concurrent.futures import ThreadPoolExecutor
except ImportError:
ThreadPoolExecutor = None
try:
from concurrent.futures import ProcessPoolExecutor
except ImportError:
ProcessPoolExecutor = None
from mpi4py.futures import MPIPoolExecutor
PRIMES = [
112272535095293,
112582705942171,
112272535095293,
115280095190773,
115797848077099,
117450548693743,
993960000099397,
]
def is_prime(n):
if n % 2 == 0:
return False
sqrt_n = int(math.floor(math.sqrt(n)))
for i in range(3, sqrt_n + 1, 2):
if n % i == 0:
return False
return True
def sequential():
return list(map(is_prime, PRIMES))
def with_thread_pool_executor():
if not ThreadPoolExecutor: return None
with ThreadPoolExecutor(4) as executor:
return list(executor.map(is_prime, PRIMES))
def with_process_pool_executor():
if not ProcessPoolExecutor: return None
with ProcessPoolExecutor(4) as executor:
return list(executor.map(is_prime, PRIMES))
def with_mpi_pool_executor():
with MPIPoolExecutor(4) as executor:
return list(executor.map(is_prime, PRIMES))
def main():
for name, fn in [('sequential', sequential),
('threads', with_thread_pool_executor),
('processes', with_process_pool_executor),
('mpi4py', with_mpi_pool_executor)]:
sys.stdout.write('%s: ' % name.ljust(11))
sys.stdout.flush()
start = time.time()
result = fn()
if result is None:
sys.stdout.write(' not available\n')
elif result != [True] * len(PRIMES):
sys.stdout.write(' failed\n')
else:
sys.stdout.write('%5.2f seconds\n' % (time.time() - start))
sys.stdout.flush()
if __name__ == '__main__':
main()
|
