为什么NumPy中的矩阵乘法不调用__add__方法？

Question

NumPy中的class matrix中的方法__mul__是如何工作的呢？ 我想实现二进制矩阵乘法，我有一个 class Binary ：

class Binary(int):
    def __init__(self, val):
        if val != 0:
            self.val = 1
        else:
            self.val = 0

    def __add__(self, other):
        print('add')
        return self.val ^ other

    def __radd__(self, other):
        print('radd')
        return self.val ^ other

我的测试：

from Binary import Binary
from numpy import matrix

i = Binary(1)
o = Binary(0)
a = matrix([i, i, o, i, i, o, o], dtype=Binary)
b = matrix([[o, o, i],
           [o, i, o],
           [o, i, i],
           [i, o, o],
           [i, o, i],
           [i, i, o],
           [i, i, i]], dtype=Binary)
print(a * b)

结果：

./test.py
[[2 1 2]]

不使用方法__add__ 。 而矩阵乘法中有求和。