The log1p functions compute ln(1+y) accurately, even for tiny y.