gather code blocks into functions and annotate code

[rmcminds]

relatively self-explanatory, unlike the current code. things have gotten complex and are difficult even for me to follow. In particular, the block of code that partitions variance among subfactors is a mess.