I have a numpy array list of matrices, i.e a "3d matrix" if that exists.
x = np.array([[1,2,3],[4,5,6],[7,8,9]]) print(x[:,:,None]*x[:,None,:]) [[[