Internal function to sample cluster membership indicator
polyurncpp(e, muA0, sigma0, A, sigma, sigmaInv, sigma0Inv, c)
Useful links