\frac{\left(1 - x\right) \cdot \left(3 - x\right)}{y \cdot 3}\left(1 - x\right) \cdot \left(\frac{3 - x}{3} \cdot \frac{1}{y}\right)double f(double x, double y) {
double r221353 = 1.0;
double r221354 = x;
double r221355 = r221353 - r221354;
double r221356 = 3.0;
double r221357 = r221356 - r221354;
double r221358 = r221355 * r221357;
double r221359 = y;
double r221360 = r221359 * r221356;
double r221361 = r221358 / r221360;
return r221361;
}
double f(double x, double y) {
double r221362 = 1.0;
double r221363 = x;
double r221364 = r221362 - r221363;
double r221365 = 3.0;
double r221366 = r221365 - r221363;
double r221367 = r221366 / r221365;
double r221368 = 1.0;
double r221369 = y;
double r221370 = r221368 / r221369;
double r221371 = r221367 * r221370;
double r221372 = r221364 * r221371;
return r221372;
}




Bits error versus x




Bits error versus y
Results
| Original | 5.7 |
|---|---|
| Target | 0.1 |
| Herbie | 0.2 |
Initial program 5.7
rmApplied times-frac0.1
rmApplied div-inv0.2
Applied associate-*l*0.2
Simplified0.1
rmApplied div-inv0.2
Final simplification0.2
herbie shell --seed 2019362 +o rules:numerics
(FPCore (x y)
:name "Diagrams.TwoD.Arc:bezierFromSweepQ1 from diagrams-lib-1.3.0.3"
:precision binary64
:herbie-target
(* (/ (- 1 x) y) (/ (- 3 x) 3))
(/ (* (- 1 x) (- 3 x)) (* y 3)))