NumPy数组上带有+ =的意外结果

Question

我正在使用NumPy在Python中使用标准方法创建对称矩阵/数组：

x = rand(500,500)
x = (x+x.T)
all(x==x.T)
> True

现在让我们变得聪明：

x = rand(500,500)
x += x.T
all(x==x.T)
> False

等一下

x==x.T
> array([[ True,  True,  True, ..., False, False, False],
       [ True,  True,  True, ..., False, False, False],
       [ True,  True,  True, ..., False, False, False],
       ..., 
       [False, False, False, ...,  True,  True,  True],
       [False, False, False, ...,  True,  True,  True],
       [False, False, False, ...,  True,  True,  True]], dtype=bool)

左上和右下部分是对称的。如果我选择了较小的阵列怎么办？

x = rand(50,50)
x += x.T
all(x==x.T)
> True

好....

x = rand(90,90)
x += x.T
all(x==x.T)
> True

x = rand(91,91)
x += x.T
all(x==x.T)
> False

只是要确保…

x = rand(91,91)
x = (x+x.T)
all(x==x.T)
> True

这是一个错误，还是我要学习一些有关+=NumPy数组的疯狂知识？

Answer 1

该transpose操作返回该数组的
视图
，这意味着未分配新数组。反过来，这意味着您正在同时读取和修改阵列。很难说出为什么某些大小或结果区域起作用的原因，但是最有可能与numpy处理数组加法（也许是制作子矩阵的副本）和/或数组视图（也许是它创建的小尺寸）有关。一个新数组）。

该x = x + x.T操作之所以有效，是因为您在那里在那里创建了一个新数组，然后将其分配给x。