x + \frac{y \cdot \left(z - t\right)}{a - t}\mathsf{fma}\left(\frac{z - t}{a - t}, y, x\right)double f(double x, double y, double z, double t, double a) {
double r375721 = x;
double r375722 = y;
double r375723 = z;
double r375724 = t;
double r375725 = r375723 - r375724;
double r375726 = r375722 * r375725;
double r375727 = a;
double r375728 = r375727 - r375724;
double r375729 = r375726 / r375728;
double r375730 = r375721 + r375729;
return r375730;
}
double f(double x, double y, double z, double t, double a) {
double r375731 = z;
double r375732 = t;
double r375733 = r375731 - r375732;
double r375734 = a;
double r375735 = r375734 - r375732;
double r375736 = r375733 / r375735;
double r375737 = y;
double r375738 = x;
double r375739 = fma(r375736, r375737, r375738);
return r375739;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
| Original | 11.0 |
|---|---|
| Target | 1.3 |
| Herbie | 1.4 |
Initial program 11.0
Simplified3.0
rmApplied clear-num3.3
rmApplied fma-udef3.3
Simplified3.1
rmApplied associate-/r/1.4
Applied fma-def1.4
Final simplification1.4
herbie shell --seed 2019235 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTicks from plot-0.2.3.4, B"
:precision binary64
:herbie-target
(+ x (/ y (/ (- a t) (- z t))))
(+ x (/ (* y (- z t)) (- a t))))