One more corner case here: for q\to 0, `isf` can lose precision (generic code path is `self._ppf(1-q, *args)`). Add an explicit `_isf` method, using the fact that `log1p(-(1-q)) = log(q)`?