x + y \cdot \frac{z - t}{a - t}x + y \cdot \left(\left(z - t\right) \cdot \frac{1}{a - t}\right)double f(double x, double y, double z, double t, double a) {
double r633790 = x;
double r633791 = y;
double r633792 = z;
double r633793 = t;
double r633794 = r633792 - r633793;
double r633795 = a;
double r633796 = r633795 - r633793;
double r633797 = r633794 / r633796;
double r633798 = r633791 * r633797;
double r633799 = r633790 + r633798;
return r633799;
}
double f(double x, double y, double z, double t, double a) {
double r633800 = x;
double r633801 = y;
double r633802 = z;
double r633803 = t;
double r633804 = r633802 - r633803;
double r633805 = 1.0;
double r633806 = a;
double r633807 = r633806 - r633803;
double r633808 = r633805 / r633807;
double r633809 = r633804 * r633808;
double r633810 = r633801 * r633809;
double r633811 = r633800 + r633810;
return r633811;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 1.2 |
|---|---|
| Target | 0.4 |
| Herbie | 1.2 |
Initial program 1.2
rmApplied div-inv1.2
Final simplification1.2
herbie shell --seed 2020047 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisLine from plot-0.2.3.4, B"
:precision binary64
:herbie-target
(if (< y -8.508084860551241e-17) (+ x (* y (/ (- z t) (- a t)))) (if (< y 2.894426862792089e-49) (+ x (* (* y (- z t)) (/ 1 (- a t)))) (+ x (* y (/ (- z t) (- a t))))))
(+ x (* y (/ (- z t) (- a t)))))