in submix.py [0:0]
def _lambda_solver_iteration(self, p, p_prime, p0):
lamb = 1.0
eps = float('inf')
while eps > self.target_consumption:
lamb = self.lambda_dec_factor*lamb
mix_star = self.mix((p + p_prime)/2, p0, lamb=lamb)
p_mix = self.mix(p, p0, lamb=lamb)
p_prime_mix = self.mix(p_prime, p0, lamb=lamb)
eps = max(self.renyi_dp(mix_star,p_mix, alpha=self.alpha),
self.renyi_dp(mix_star,p_prime_mix, alpha=self.alpha))
return mix_star, lamb