x + \frac{y \cdot \left(z - t\right)}{z - a}\mathsf{fma}\left(\frac{z - t}{z - a}, y, x\right)double f(double x, double y, double z, double t, double a) {
double r25180807 = x;
double r25180808 = y;
double r25180809 = z;
double r25180810 = t;
double r25180811 = r25180809 - r25180810;
double r25180812 = r25180808 * r25180811;
double r25180813 = a;
double r25180814 = r25180809 - r25180813;
double r25180815 = r25180812 / r25180814;
double r25180816 = r25180807 + r25180815;
return r25180816;
}
double f(double x, double y, double z, double t, double a) {
double r25180817 = z;
double r25180818 = t;
double r25180819 = r25180817 - r25180818;
double r25180820 = a;
double r25180821 = r25180817 - r25180820;
double r25180822 = r25180819 / r25180821;
double r25180823 = y;
double r25180824 = x;
double r25180825 = fma(r25180822, r25180823, r25180824);
return r25180825;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
| Original | 10.3 |
|---|---|
| Target | 1.2 |
| Herbie | 1.3 |
Initial program 10.3
Simplified3.1
rmApplied clear-num3.3
rmApplied fma-udef3.3
Simplified3.1
rmApplied associate-/r/1.3
Applied fma-def1.3
Final simplification1.3
herbie shell --seed 2019163 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTicks from plot-0.2.3.4, A"
:herbie-target
(+ x (/ y (/ (- z a) (- z t))))
(+ x (/ (* y (- z t)) (- z a))))