As suggested by rhettinger in http://bugs.python.org/msg294438:

I agree that the comment should be changed. While we at it, perhaps sync-up with expovariate() code and eliminate the ``u <= 1e-7`` test:

Instead of:

        elif alpha == 1.0:
            # expovariate(1)
            u = random()
            while u <= 1e-7:
                u = random()
            return -_log(u) * beta

Use this instead:

        elif alpha == 1.0:
            # expovariate(1.0 / beta)
            return -_log(1.0 - random()) * beta

New changeset 63d152232e1742660f481c04a811f824b91f6790 by Raymond Hettinger (leodema) in branch 'master':
bpo-30561: Sync-up expovariate() and gammavariate code (GH-1934)
https://github.com/python/cpython/commit/63d152232e1742660f481c04a811f824b91f6790