update_z2_Sj function

Collapsed Gibbs for zz using matrix inversion lemma