x + y \cdot \frac{z - t}{a - t}x + \frac{y}{\frac{a - t}{z - t}}double f(double x, double y, double z, double t, double a) {
double r578089 = x;
double r578090 = y;
double r578091 = z;
double r578092 = t;
double r578093 = r578091 - r578092;
double r578094 = a;
double r578095 = r578094 - r578092;
double r578096 = r578093 / r578095;
double r578097 = r578090 * r578096;
double r578098 = r578089 + r578097;
return r578098;
}
double f(double x, double y, double z, double t, double a) {
double r578099 = x;
double r578100 = y;
double r578101 = a;
double r578102 = t;
double r578103 = r578101 - r578102;
double r578104 = z;
double r578105 = r578104 - r578102;
double r578106 = r578103 / r578105;
double r578107 = r578100 / r578106;
double r578108 = r578099 + r578107;
return r578108;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 1.4 |
|---|---|
| Target | 0.4 |
| Herbie | 1.3 |
Initial program 1.4
rmApplied clear-num1.4
rmApplied un-div-inv1.3
Final simplification1.3
herbie shell --seed 2020065 +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)))))