x \cdot \left(1 - x \cdot y\right)
x \cdot \left(1 - x \cdot y\right)
double f(double x, double y) {
double r57260 = x;
double r57261 = 1.0;
double r57262 = y;
double r57263 = r57260 * r57262;
double r57264 = r57261 - r57263;
double r57265 = r57260 * r57264;
return r57265;
}
double f(double x, double y) {
double r57266 = x;
double r57267 = 1.0;
double r57268 = y;
double r57269 = r57266 * r57268;
double r57270 = r57267 - r57269;
double r57271 = r57266 * r57270;
return r57271;
}



Bits error versus x



Bits error versus y
Results
Initial program 0.1
Final simplification0.1
herbie shell --seed 2019304 +o rules:numerics
(FPCore (x y)
:name "Numeric.SpecFunctions:log1p from math-functions-0.1.5.2, A"
:precision binary64
(* x (- 1 (* x y))))