Sigma_star_cluster function

Calculates parameter Σ\Sigma^* of q(β)q^*(\beta) to optimize the evidence based lower bound (ELBO) in survregVB.frailty.fit.