@pbrod: I can't reproduce the `np.log1p` bug.  Using both 1.7.1 and 1.8.0, `np.log1p(np.inf)` return `inf`.  Perhaps it is platform dependent? (I'm using Ubuntu 12.04 (64 bit).) Have you reported the problem on the numpy issue tracker (https://github.com/numpy/numpy/issues)?