x + \frac{y \cdot \left(z - x\right)}{t}\frac{y}{t} \cdot \left(z - x\right) + xdouble f(double x, double y, double z, double t) {
double r1753 = x;
double r1754 = y;
double r1755 = z;
double r1756 = r1755 - r1753;
double r1757 = r1754 * r1756;
double r1758 = t;
double r1759 = r1757 / r1758;
double r1760 = r1753 + r1759;
return r1760;
}
double f(double x, double y, double z, double t) {
double r1761 = y;
double r1762 = t;
double r1763 = r1761 / r1762;
double r1764 = z;
double r1765 = x;
double r1766 = r1764 - r1765;
double r1767 = r1763 * r1766;
double r1768 = r1767 + r1765;
return r1768;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.4 |
|---|---|
| Target | 2.1 |
| Herbie | 2.1 |
Initial program 6.4
Simplified2.1
rmApplied fma-udef2.1
Final simplification2.1
herbie shell --seed 2020025 +o rules:numerics
(FPCore (x y z t)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, D"
:precision binary64
:herbie-target
(- x (+ (* x (/ y t)) (* (- z) (/ y t))))
(+ x (/ (* y (- z x)) t)))