Lots of stuff.
1. Implement the S_star cpu sampling function
2. Notice that it's really slow, try to optimize S calculation by overloading to allow partial recalculations.
3. Realize S calculation is already blazing fast
4. Decide to commit what I have so far, and dive into the full conditional for S star