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 r374273 = x;
double r374274 = y;
double r374275 = z;
double r374276 = t;
double r374277 = r374275 - r374276;
double r374278 = r374274 * r374277;
double r374279 = a;
double r374280 = r374279 - r374276;
double r374281 = r374278 / r374280;
double r374282 = r374273 + r374281;
return r374282;
}
double f(double x, double y, double z, double t, double a) {
double r374283 = z;
double r374284 = t;
double r374285 = r374283 - r374284;
double r374286 = a;
double r374287 = r374286 - r374284;
double r374288 = r374285 / r374287;
double r374289 = y;
double r374290 = x;
double r374291 = fma(r374288, r374289, r374290);
return r374291;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
| Original | 10.9 |
|---|---|
| Target | 1.2 |
| Herbie | 1.3 |
Initial program 10.9
Simplified2.9
rmApplied div-inv2.9
rmApplied add-cube-cbrt3.4
Applied associate-/r*3.4
rmApplied fma-udef3.4
Simplified1.3
rmApplied fma-def1.3
Final simplification1.3
herbie shell --seed 2019303 +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))))