numpy.transpose在大小为＆gt;的矩阵中的表观错误90x90

Question

将矩形矩阵的转置添加到该矩阵应该产生围绕对角线对称的矩阵。例如：

import numpy as np
W1 = np.random.random((90, 90)) 
W1 += np.transpose(W1)
(np.transpose(W1) == W1).all()

  

真

但是这种行为似乎在更大的矩阵中崩溃了：

W2 = np.random.random((91, 91)) 
W2 += np.transpose(W2)
(np.transpose(W2) == W2).all()

  

假

差异并非微不足道：

print( np.max(W2), np.max(W2-np.transpose(W2)), np.min(W2-np.transpose(W2)))

  

2.84604132441 0.987790051685 -0.987790051685

这是numpy的缺点还是我错过了什么？