如何比较持有 numpy.ndarray 的数据类的相等性(bool(a==b) 引发 ValueError)?

jim*_*qaz 5 python numpy python-dataclasses

如果我创建一个包含 Numpy ndarray 的 Python 数据类,我将无法再使用自动生成的数据类__eq__

import numpy as np

@dataclass
class Instr:
    foo: np.ndarray
    bar: np.ndarray

arr = np.array([1])
arr2 = np.array([1, 2])
print(Instr(arr, arr) == Instr(arr2, arr2))
Run Code Online (Sandbox Code Playgroud)

ValueError:包含多个元素的数组的真值不明确。使用 a.any() 或 a.all()

这是因为ndarray.__eq__ 有时返回ndarray真值,通过比较a[0]b[0],等等等等到2。这是相当复杂的,不直观的时间越长,而实际上只提高在阵列不同的形状,或者有错误不同的价值观什么的。

如何安全地比较@dataclasses 持有 Numpy 数组?


@dataclass的实现__eq__是使用生成的eval()。它的源在堆栈跟踪中丢失,无法使用 来查看inspect,但它实际上使用的是调用 bool(foo)的元组比较

import dis
dis.dis(Instr.__eq__)
Run Code Online (Sandbox Code Playgroud)

摘抄:

  3          12 LOAD_FAST                0 (self)
             14 LOAD_ATTR                1 (foo)
             16 LOAD_FAST                0 (self)
             18 LOAD_ATTR                2 (bar)
             20 BUILD_TUPLE              2
             22 LOAD_FAST                1 (other)
             24 LOAD_ATTR                1 (foo)
             26 LOAD_FAST                1 (other)
             28 LOAD_ATTR                2 (bar)
             30 BUILD_TUPLE              2
             32 COMPARE_OP               2 (==)
             34 RETURN_VALUE
Run Code Online (Sandbox Code Playgroud)

FHT*_*ell 5

解决方案是放入您自己的__eq__方法并进行设置,eq=False以便数据类不会生成自己的方法(尽管检查最后一步不是必需的文档,但我认为无论如何还是明确的很好)。

import numpy as np

def array_eq(arr1, arr2):
    return (isinstance(arr1, np.ndarray) and
            isinstance(arr2, np.ndarray) and
            arr1.shape == arr2.shape and
            (arr1 == arr2).all())

@dataclass(eq=False)
class Instr:

    foo: np.ndarray
    bar: np.ndarray

    def __eq__(self, other):
        if not isinstance(other, Instr):
            return NotImplemented
        return array_eq(self.foo, other.foo) and array_eq(self.bar, other.bar)
Run Code Online (Sandbox Code Playgroud)

编辑

通用数据类的通用快速解决方案,其中一些值是 numpy 数组,而另一些不是

import numpy as np
from dataclasses import dataclass, astuple

def array_safe_eq(a, b) -> bool:
    """Check if a and b are equal, even if they are numpy arrays"""
    if a is b:
        return True
    if isinstance(a, np.ndarray) and isinstance(b, np.ndarray):
        return a.shape == b.shape and (a == b).all()
    try:
        return a == b
    except TypeError:
        return NotImplemented

def dc_eq(dc1, dc2) -> bool:
   """checks if two dataclasses which hold numpy arrays are equal"""
   if dc1 is dc2:
        return True
   if dc1.__class__ is not dc2.__class__:
       return NotImplmeneted  # better than False
   t1 = astuple(dc1)
   t2 = astuple(dc2)
   return all(array_safe_eq(a1, a2) for a1, a2 in zip(t1, t2))

# usage
@dataclass(eq=False)
class T:

   a: int
   b: np.ndarray
   c: np.ndarray

   def __eq__(self, other):
        return dc_eq(self, other)
Run Code Online (Sandbox Code Playgroud)