如何有效地计算中位数成对平方欧氏距离

Question

我在寻找输入数组的中位数成对平方欧氏距离

我知道我可以使用pdist来完成：

from scipy.spatial.distance import pdist
dists = pdist(inp,'sqeuclidean')
np.median(dists)

但是我需要在非常大的输入阵列上重复这个过程很多次。我怀疑这可能会非常快，因为成对距离的计算是O（n^2）——除非有聪明的方法。但是我知道pdist的代码不是特别有效（源代码是嵌套循环的），所以我想知道是否有办法加快它，特别是考虑到我只需要中间值