x \cdot \left(1 - x \cdot y\right)
x \cdot \left(1 - x \cdot y\right)
double f(double x, double y) {
double r56552 = x;
double r56553 = 1.0;
double r56554 = y;
double r56555 = r56552 * r56554;
double r56556 = r56553 - r56555;
double r56557 = r56552 * r56556;
return r56557;
}
double f(double x, double y) {
double r56558 = x;
double r56559 = 1.0;
double r56560 = y;
double r56561 = r56558 * r56560;
double r56562 = r56559 - r56561;
double r56563 = r56558 * r56562;
return r56563;
}



Bits error versus x



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