x \cdot \left(1 + y \cdot y\right)
1 \cdot x + y \cdot \left(y \cdot x\right)
double f(double x, double y) {
double r325947 = x;
double r325948 = 1.0;
double r325949 = y;
double r325950 = r325949 * r325949;
double r325951 = r325948 + r325950;
double r325952 = r325947 * r325951;
return r325952;
}
double f(double x, double y) {
double r325953 = 1.0;
double r325954 = x;
double r325955 = r325953 * r325954;
double r325956 = y;
double r325957 = r325956 * r325954;
double r325958 = r325956 * r325957;
double r325959 = r325955 + r325958;
return r325959;
}




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
Simplified5.3
rmApplied associate-*l*0.1
Final simplification0.1
herbie shell --seed 2019303 +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))))