However, scipy.special.log1p gives the correct answer on my system:

In [2]: scipy.special.log1p(inf)
Out[2]: inf