def _rvs(self, a):
        sz, rndm = self._size, self._random_state
        u = rndm.random_sample(size=sz)
        gm = gamma.rvs(a, size=sz, random_state=rndm)
        return gm * where(u >= 0.5, 1, -1)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/dgamma_gen._rvs/original