x - \frac{\left(y \cdot 2\right) \cdot z}{\left(z \cdot 2\right) \cdot z - y \cdot t}x - \frac{y \cdot 2}{z \cdot 2 - \frac{\frac{t}{z}}{\frac{1}{y}}}double f(double x, double y, double z, double t) {
double r418811 = x;
double r418812 = y;
double r418813 = 2.0;
double r418814 = r418812 * r418813;
double r418815 = z;
double r418816 = r418814 * r418815;
double r418817 = r418815 * r418813;
double r418818 = r418817 * r418815;
double r418819 = t;
double r418820 = r418812 * r418819;
double r418821 = r418818 - r418820;
double r418822 = r418816 / r418821;
double r418823 = r418811 - r418822;
return r418823;
}
double f(double x, double y, double z, double t) {
double r418824 = x;
double r418825 = y;
double r418826 = 2.0;
double r418827 = r418825 * r418826;
double r418828 = z;
double r418829 = r418828 * r418826;
double r418830 = t;
double r418831 = r418830 / r418828;
double r418832 = 1.0;
double r418833 = r418832 / r418825;
double r418834 = r418831 / r418833;
double r418835 = r418829 - r418834;
double r418836 = r418827 / r418835;
double r418837 = r418824 - r418836;
return r418837;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 11.8 |
|---|---|
| Target | 0.1 |
| Herbie | 1.0 |
Initial program 11.8
rmApplied associate-/l*7.0
rmApplied div-sub7.0
Simplified2.9
Simplified2.9
rmApplied associate-/l*2.1
rmApplied div-inv2.1
Applied associate-/r*1.0
Final simplification1.0
herbie shell --seed 2020039 +o rules:numerics
(FPCore (x y z t)
:name "Numeric.AD.Rank1.Halley:findZero from ad-4.2.4"
:precision binary64
:herbie-target
(- x (/ 1 (- (/ z y) (/ (/ t 2) z))))
(- x (/ (* (* y 2) z) (- (* (* z 2) z) (* y t)))))