x + \frac{y \cdot \left(z - x\right)}{t}x - \left(\frac{x \cdot y}{t} - \frac{y \cdot z}{t}\right)double f(double x, double y, double z, double t) {
double r234870 = x;
double r234871 = y;
double r234872 = z;
double r234873 = r234872 - r234870;
double r234874 = r234871 * r234873;
double r234875 = t;
double r234876 = r234874 / r234875;
double r234877 = r234870 + r234876;
return r234877;
}
double f(double x, double y, double z, double t) {
double r234878 = x;
double r234879 = y;
double r234880 = r234878 * r234879;
double r234881 = t;
double r234882 = r234880 / r234881;
double r234883 = z;
double r234884 = r234879 * r234883;
double r234885 = r234884 / r234881;
double r234886 = r234882 - r234885;
double r234887 = r234878 - r234886;
return r234887;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.6 |
|---|---|
| Target | 2.1 |
| Herbie | 6.6 |
Initial program 6.6
rmApplied associate-/l*6.0
rmApplied associate-/r/2.1
rmApplied add-cube-cbrt2.6
Applied add-cube-cbrt2.8
Applied times-frac2.8
Applied associate-*l*1.0
Final simplification6.6
herbie shell --seed 1978988140
(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)))