#!/usr/bin/env python

from multiprocessing import Pool

import numpy as np

# points per call

N = 5e6

# number of calls

M = 100

# threads

C = 8

def find_pi(m):

return np.add.reduce(x1)*4.0/N

if __name__ == '__main__':

p = Pool(C)

results = p.map(find_pi, range(M))

print reduce(lambda x,y: (x+y)/2, results)