update_SigmaINV_faster_Sj function

Gibbs sampling for Σ1\Sigma^{-1} per component