numpy / scipy数据集合之间的对距离

Question

我有一个L个阵列的列表samples，形状N乘3表示空间中的三维点。我想计算列表中每对i，j数组之间的成对距离矩阵。我现在正在做的只是循环夫妻i和j并使用cdist

for i in xrange(L-1):
    for j in xrange(i+1,L):
        dists=cdist(sample[i],sample[j])

然后，我使用获得的距离矩阵来执行进一步的操作。

有没有办法避免双循环？也许通过einsum通过合适的N乘3乘L阵列？