how can i transform blocks into a blockdiagonal matrix (numpy),"r = np.kron(np.diag([1, 2, 3]), np.ones((3, 3), dtype='int'))"