"无需分配不必要数组的Numpy运算符结果总和"

Question

我有两个numpy布尔数组（a和b）。我需要找出它们的元素有多少是相等的。目前，我做len(a) - (a ^ b).sum()，但据我所知，xor操作创建了一个全新的numpy数组。如何在不创建不必要的临时数组的情况下有效地实现这种期望的行为？在

我试过使用numexpr，但是我不能让它正常工作。它不支持True为1，False为0的概念，因此我必须使用ne.evaluate("sum(where(a==b, 1, 0))")，这大约需要两倍的时间。在

编辑：我忘记了，其中一个数组实际上是另一个大小不同的数组的视图，这两个数组都应该被认为是不可变的。两个数组都是二维的，大小都在25x40左右。在

是的，这是我程序的瓶颈，值得优化。在

Answer 1

首先，您可以跳过A*B步骤：

>>> a
array([ True, False,  True, False,  True], dtype=bool)
>>> b
array([False,  True,  True, False,  True], dtype=bool)
>>> np.sum(~(a^b))
3

如果你不介意摧毁阵法a或b，我不确定你会比这个更快：

^{pr2}$