def _mix()

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))