x + \frac{y \cdot \left(z - t\right)}{z - a}\frac{\frac{z - t}{z - a}}{\frac{1}{y}} + xdouble f(double x, double y, double z, double t, double a) {
double r718114 = x;
double r718115 = y;
double r718116 = z;
double r718117 = t;
double r718118 = r718116 - r718117;
double r718119 = r718115 * r718118;
double r718120 = a;
double r718121 = r718116 - r718120;
double r718122 = r718119 / r718121;
double r718123 = r718114 + r718122;
return r718123;
}
double f(double x, double y, double z, double t, double a) {
double r718124 = z;
double r718125 = t;
double r718126 = r718124 - r718125;
double r718127 = a;
double r718128 = r718124 - r718127;
double r718129 = r718126 / r718128;
double r718130 = 1.0;
double r718131 = y;
double r718132 = r718130 / r718131;
double r718133 = r718129 / r718132;
double r718134 = x;
double r718135 = r718133 + r718134;
return r718135;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 10.4 |
|---|---|
| Target | 1.2 |
| Herbie | 1.4 |
Initial program 10.4
Simplified2.9
rmApplied div-inv3.0
rmApplied add-cube-cbrt3.4
Applied associate-*l*3.4
Simplified3.4
rmApplied fma-udef3.4
Simplified3.0
rmApplied div-inv3.1
Applied associate-/r*1.4
Final simplification1.4
herbie shell --seed 2019362 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTicks from plot-0.2.3.4, A"
:precision binary64
:herbie-target
(+ x (/ y (/ (- z a) (- z t))))
(+ x (/ (* y (- z t)) (- z a))))