from multiprocessing import Pool
from itertools import product
output = np.zeros((N, N))
pool = Pool()
chunksize = 20
for ind, res in enumerate(pool.imap(Fun, product(range(N), range(N))),
    chunksize):
    output.flat[ind] = res
Parfor for Python