x \cdot \left(1 + y \cdot y\right)
x \cdot 1 + \left(x \cdot y\right) \cdot y
double f(double x, double y) {
double r451045 = x;
double r451046 = 1.0;
double r451047 = y;
double r451048 = r451047 * r451047;
double r451049 = r451046 + r451048;
double r451050 = r451045 * r451049;
return r451050;
}
double f(double x, double y) {
double r451051 = x;
double r451052 = 1.0;
double r451053 = r451051 * r451052;
double r451054 = y;
double r451055 = r451051 * r451054;
double r451056 = r451055 * r451054;
double r451057 = r451053 + r451056;
return r451057;
}




Bits error versus x




Bits error versus y
Results
| Original | 5.4 |
|---|---|
| Target | 0.1 |
| Herbie | 0.1 |
Initial program 5.4
rmApplied distribute-lft-in5.4
rmApplied associate-*r*0.1
Final simplification0.1
herbie shell --seed 2020064 +o rules:numerics
(FPCore (x y)
:name "Numeric.Integration.TanhSinh:everywhere from integration-0.2.1"
:precision binary64
:herbie-target
(+ x (* (* x y) y))
(* x (+ 1 (* y y))))