x \cdot \left(1 + y \cdot y\right)
x \cdot 1 + y \cdot \left(y \cdot x\right)
double f(double x, double y) {
double r582387 = x;
double r582388 = 1.0;
double r582389 = y;
double r582390 = r582389 * r582389;
double r582391 = r582388 + r582390;
double r582392 = r582387 * r582391;
return r582392;
}
double f(double x, double y) {
double r582393 = x;
double r582394 = 1.0;
double r582395 = r582393 * r582394;
double r582396 = y;
double r582397 = r582396 * r582393;
double r582398 = r582396 * r582397;
double r582399 = r582395 + r582398;
return r582399;
}




Bits error versus x




Bits error versus y
Results
| Original | 5.3 |
|---|---|
| Target | 0.1 |
| Herbie | 0.1 |
Initial program 5.3
rmApplied distribute-lft-in5.3
Simplified5.3
rmApplied associate-*l*0.1
Final simplification0.1
herbie shell --seed 2020047 +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))))