def _ppf(self, q, p):
        vals = ceil(log(1.0-q)/log(1-p))
        temp = self._cdf(vals-1, p)
        return np.where((temp >= q) & (vals > 0), 1-vals, vals)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_discrete_distns.py/geom_gen._ppf/swapped operands of `-`