omega_star function

Calculates parameter ω\omega^* of q(b)q^*(b) to optimize the evidence based lower bound (ELBO) in survregVB.fit.