Update overflow bound in log1p(f)

And add tests in this region.

Fix functions in region between previous
overflow bound and maximum float.

Accuracy is maintained by relying on high
accuracy log (in branch), however performance
drops slightly. Ideally the core approximation
should also work above threshold.