\mathsf{log1p}\left(\frac{x}{y}\right)\mathsf{log1p}\left(\frac{x}{y}\right)double code(double x, double y) {
return ((double) log1p(((double) (x / y))));
}
double code(double x, double y) {
return ((double) log1p(((double) (x / y))));
}



Bits error versus x



Bits error versus y
Results
Initial program 5.2
Final simplification5.2
herbie shell --seed 2020153
(FPCore (x y)
:name "(log1p (/ x y))"
:precision binary64
(log1p (/ x y)))