Tim Peters <t...@python.org> added the comment:
BTW, if we have to "do something", how about changing return 1.0 / u ** (1.0/alpha) to the mathematically equivalent return (1.0 / u) ** (1.0/alpha) ? Not sure about Linux-y boxes, but on Windows that would raise OverflowError instead of ZeroDivisionError. Which makes more sense, because the Pareto variable we're _trying_ to produce is too large to represent. If it returns math.inf on some box instead, same thing. Or, also faster, and suffering fewer rounding errors, return u ** (-1.0 / alpha) ---------- _______________________________________ Python tracker <rep...@bugs.python.org> <https://bugs.python.org/issue41421> _______________________________________ _______________________________________________ Python-bugs-list mailing list Unsubscribe: https://mail.python.org/mailman/options/python-bugs-list/archive%40mail-archive.com
