Numpy Вычитание широковещательного массива

Question

data = np.array([[1,2,3], [2,4,6],[2,4,6]])
cs = np.array([[5,5,5],[6,6,6]])
data - cs[:, np.newaxis]

data измерение 3x3
cs измерение 2x3
cs[:, np.newaxis] измерение 2x1x3 Кто-нибудь может объяснить, как работает вычитание?