Why wouldn't the argument of `_cdf` go to `inf`?  The result should be 1.  If we use the buggy version of `np.expm1`, `_cdf(inf, c)` will return `nan`.