update_SigmaINV_faster function

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