def _ppf(self, q, lambda_):
        vals = ceil(1--1.0/lambda_ * log1p(-q))
        vals1 = (vals-1).clip(self.a, np.inf)
        temp = self._cdf(vals1, lambda_)
        return np.where(temp >= q, vals1, vals)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_discrete_distns.py/planck_gen._ppf/swapped operands of `-`