A = A/np.maximum(A.sum(1).reshape(d,1),1) # normalisation
Pi = Pi/Pi.sum()