BUG: use special.log1p instead of np.log1p

np.log1p(np.inf) is reported to produce nans on some platforms,
see numpy issue gh-4225