def _ppf(self, q, mu):
        vals = ceil(special.pdtrik(q, mu))
        vals1 = np.maximum(vals - 1, 0)
        temp = special.pdtr(vals1, mu)
        return np.where(temp >= q, vals1, vals)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_discrete_distns.py/poisson_gen._ppf/original