https://issues.dlang.org/show_bug.cgi?id=23677 --- Comment #2 from John Hall <john.michael.hall at gmail.com> --- I double-checked and log1p doesn't provide the same accuracy as the CEPHES implementation for doubles. --