x + \frac{y}{1.12837916709551256 \cdot e^{z} - x \cdot y}x + \frac{1}{1 \cdot \left(\frac{1.12837916709551256 \cdot \left(\sqrt[3]{e^{z}} \cdot \sqrt[3]{e^{z}}\right)}{\frac{y}{\sqrt[3]{e^{z}}}} - x\right)}double code(double x, double y, double z) {
return (x + (y / ((1.1283791670955126 * exp(z)) - (x * y))));
}
double code(double x, double y, double z) {
return (x + (1.0 / (1.0 * (((1.1283791670955126 * (cbrt(exp(z)) * cbrt(exp(z)))) / (y / cbrt(exp(z)))) - x))));
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.8 |
|---|---|
| Target | 0.1 |
| Herbie | 0.1 |
Initial program 2.8
rmApplied clear-num2.8
rmApplied *-un-lft-identity2.8
Applied *-un-lft-identity2.8
Applied times-frac2.8
Simplified2.8
Simplified0.1
rmApplied add-cube-cbrt0.1
Applied associate-/l*0.1
Applied associate-*r/0.1
Final simplification0.1
herbie shell --seed 2020066 +o rules:numerics
(FPCore (x y z)
:name "Numeric.SpecFunctions:invErfc from math-functions-0.1.5.2, A"
:precision binary64
:herbie-target
(+ x (/ 1 (- (* (/ 1.1283791670955126 y) (exp z)) x)))
(+ x (/ y (- (* 1.1283791670955126 (exp z)) (* x y)))))