x + \frac{\left(y - z\right) \cdot t}{a - z}\left(\frac{1}{\frac{a - z}{y}} - \frac{z}{a - z}\right) \cdot t + xdouble f(double x, double y, double z, double t, double a) {
double r675621 = x;
double r675622 = y;
double r675623 = z;
double r675624 = r675622 - r675623;
double r675625 = t;
double r675626 = r675624 * r675625;
double r675627 = a;
double r675628 = r675627 - r675623;
double r675629 = r675626 / r675628;
double r675630 = r675621 + r675629;
return r675630;
}
double f(double x, double y, double z, double t, double a) {
double r675631 = 1.0;
double r675632 = a;
double r675633 = z;
double r675634 = r675632 - r675633;
double r675635 = y;
double r675636 = r675634 / r675635;
double r675637 = r675631 / r675636;
double r675638 = r675633 / r675634;
double r675639 = r675637 - r675638;
double r675640 = t;
double r675641 = r675639 * r675640;
double r675642 = x;
double r675643 = r675641 + r675642;
return r675643;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 10.5 |
|---|---|
| Target | 0.5 |
| Herbie | 1.3 |
Initial program 10.5
Simplified1.2
rmApplied div-sub1.2
rmApplied fma-udef1.2
rmApplied clear-num1.3
Final simplification1.3
herbie shell --seed 2019351 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTick from plot-0.2.3.4, A"
:precision binary64
:herbie-target
(if (< t -1.0682974490174067e-39) (+ x (* (/ (- y z) (- a z)) t)) (if (< t 3.9110949887586375e-141) (+ x (/ (* (- y z) t) (- a z))) (+ x (* (/ (- y z) (- a z)) t))))
(+ x (/ (* (- y z) t) (- a z))))