x \cdot \left(1 + y \cdot y\right)
x \cdot 1 + \left(y \cdot x\right) \cdot y
double f(double x, double y) {
double r460283 = x;
double r460284 = 1.0;
double r460285 = y;
double r460286 = r460285 * r460285;
double r460287 = r460284 + r460286;
double r460288 = r460283 * r460287;
return r460288;
}
double f(double x, double y) {
double r460289 = x;
double r460290 = 1.0;
double r460291 = r460289 * r460290;
double r460292 = y;
double r460293 = r460292 * r460289;
double r460294 = r460293 * r460292;
double r460295 = r460291 + r460294;
return r460295;
}




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
Taylor expanded around 0 5.3
Simplified0.1
Final simplification0.1
herbie shell --seed 2020045 +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))))