大规模numpy矩阵中三角形不等式的检验

当然,使用(未测试的)可以很容易地引导最里面的循环:

for i in range(N):
    for j in range(i):
        assert all(D[i,j] <= D[i,:] + D[:,j])

对于双矢量化,您可以循环通过 k (也未测试):

for k in range(N):
    row = D[k,:].reshape(1, N)
    col = D[:,k].reshape(N, 1)
    assert all(D <= row + col)

( row + col 生成一个与 D )