import numpy as np
a = np.array([[1, 0, 1], [0, 0, 0]])
b = np.array([[1, 0, 1], [1, 1, 0], [1, 1, 0]])
print(a.shape[0])
print(a.shape[1])
print(np.array([[np.sum(a[j] ^ b[i])for i in range(0, b.shape[0])] for j in range(0, a.shape[0])]))
print(np.array([[np.sum(a[1]^[1,0,1])], [np.sum(a[0]^[1, 1, 0])], [np.sum(a[0]^[1, 1, 0])]]))