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 r1979 = x;
double r1980 = y;
double r1981 = z;
double r1982 = t;
double r1983 = r1981 - r1982;
double r1984 = r1980 * r1983;
double r1985 = a;
double r1986 = r1981 - r1985;
double r1987 = r1984 / r1986;
double r1988 = r1979 + r1987;
return r1988;
}
double f(double x, double y, double z, double t, double a) {
double r1989 = z;
double r1990 = t;
double r1991 = r1989 - r1990;
double r1992 = a;
double r1993 = r1989 - r1992;
double r1994 = r1991 / r1993;
double r1995 = y;
double r1996 = x;
double r1997 = fma(r1994, r1995, r1996);
return r1997;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
| Original | 11.1 |
|---|---|
| Target | 1.3 |
| Herbie | 1.4 |
Initial program 11.1
Simplified2.7
rmApplied clear-num3.0
rmApplied fma-udef3.0
Simplified2.8
rmApplied associate-/r/1.4
Applied fma-def1.4
Final simplification1.4
herbie shell --seed 2020025 +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))))