x - \frac{\left(y \cdot 2\right) \cdot z}{\left(z \cdot 2\right) \cdot z - y \cdot t}x - \frac{1}{\frac{2 \cdot \frac{z}{y} - \frac{t}{z}}{2}}double f(double x, double y, double z, double t) {
double r461212 = x;
double r461213 = y;
double r461214 = 2.0;
double r461215 = r461213 * r461214;
double r461216 = z;
double r461217 = r461215 * r461216;
double r461218 = r461216 * r461214;
double r461219 = r461218 * r461216;
double r461220 = t;
double r461221 = r461213 * r461220;
double r461222 = r461219 - r461221;
double r461223 = r461217 / r461222;
double r461224 = r461212 - r461223;
return r461224;
}
double f(double x, double y, double z, double t) {
double r461225 = x;
double r461226 = 1.0;
double r461227 = 2.0;
double r461228 = z;
double r461229 = y;
double r461230 = r461228 / r461229;
double r461231 = r461227 * r461230;
double r461232 = t;
double r461233 = r461232 / r461228;
double r461234 = r461231 - r461233;
double r461235 = r461234 / r461227;
double r461236 = r461226 / r461235;
double r461237 = r461225 - r461236;
return r461237;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 11.8 |
|---|---|
| Target | 0.1 |
| Herbie | 0.1 |
Initial program 11.8
rmApplied associate-/l*6.8
rmApplied *-un-lft-identity6.8
Applied *-un-lft-identity6.8
Applied times-frac6.8
Simplified6.8
Simplified2.9
rmApplied clear-num3.0
Simplified2.9
Taylor expanded around 0 0.1
Final simplification0.1
herbie shell --seed 2020036 +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)))))