def _lambda_solver_iteration()

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