Multiply high order matrices with numpy,"np.tensordot(ind, dist, axes=[1, 1])[0].T"