in submix.py [0:0]
def _mix (self, p, p_prime, p0, 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)
return mix_star, max(self.renyi_dp(mix_star,p_mix, alpha=self.alpha),
self.renyi_dp(mix_star,p_prime_mix, alpha=self.alpha))