x + \frac{y}{1.12837916709551256 \cdot e^{z} - x \cdot y}x + \frac{1}{1.12837916709551256 \cdot \frac{e^{z}}{y} - \frac{x}{1}}double f(double x, double y, double z) {
double r501982 = x;
double r501983 = y;
double r501984 = 1.1283791670955126;
double r501985 = z;
double r501986 = exp(r501985);
double r501987 = r501984 * r501986;
double r501988 = r501982 * r501983;
double r501989 = r501987 - r501988;
double r501990 = r501983 / r501989;
double r501991 = r501982 + r501990;
return r501991;
}
double f(double x, double y, double z) {
double r501992 = x;
double r501993 = 1.0;
double r501994 = 1.1283791670955126;
double r501995 = z;
double r501996 = exp(r501995);
double r501997 = y;
double r501998 = r501996 / r501997;
double r501999 = r501994 * r501998;
double r502000 = r501992 / r501993;
double r502001 = r501999 - r502000;
double r502002 = r501993 / r502001;
double r502003 = r501992 + r502002;
return r502003;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.9 |
|---|---|
| Target | 0.0 |
| Herbie | 0.0 |
Initial program 2.9
rmApplied clear-num2.9
rmApplied div-sub2.9
Simplified2.9
Simplified0.0
Final simplification0.0
herbie shell --seed 2020081 +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)))))