Parfor для Python

Я ищу окончательный ответ на parfor MATLAB для Python (Scipy, Numpy).

Есть ли решение, подобное parfor? Если не, в чем сложность его создания?

ОБНОВЛЕНИЕ: Вот типичный код числовых вычислений, который мне нужно ускорить

import numpy as np

N = 2000
output = np.zeros([N,N])
for i in range(N):
    for j in range(N):
        output[i,j] = HeavyComputationThatIsThreadSafe(i,j)

Пример сложной вычислительной функции:

import scipy.optimize

def HeavyComputationThatIsThreadSafe(i,j):
    n = i * j

    return scipy.optimize.anneal(lambda x: np.sum((x-np.arange(n)**2)), np.random.random((n,1)))[0][0,0]