Используйте внешнее вычитание, а затем измените форму -
np.subtract.outer(A,B).reshape((A.shape)+(-1,))
Или увеличьте A
до 3D
с единственным оттенком в качестве последнего и вычтите плоский коэффициент B
-
A[...,None] - B.ravel() # or B.flat in place of B.ravel()