Copyright Digital Equipment Corp. All rights reserved.
The log1p functions compute ln(1+y) accurately, even for tiny y.