It does not seem to be directly related to this PR --- `expm1` is only used in `_cdf`, and the argument does not go to positive inf. (BTW `boxcox` has been reimplemeted in terms of `expm1`, is it affected? ping @WarrenWeckesser)
But I agree, it might be that `expm1` & `log1p` have to be replaced everywhere in the scipy codebase. I'll open a ticket.