在python中稀疏矩阵的并行计算

我有一个大的[numpy]向量X和一个比较函数f(x,y)。我需要找到X的所有成对元素，对于它们f(X[I],X[j])<T的某个阈值T。这很好：

good_inds = {}
for i in range(0,len(X)):
   for j in range(x+1,len(X)):
       score = f(X[i],X[j])
       if score<T:
           good_inds[x,y] = score

实际上，这会构建一个字典，该字典表示稀疏矩阵。问题是它相当慢，我希望并行化此过程。请指教。