x + \frac{y}{1.128379167095512558560699289955664426088 \cdot e^{z} - x \cdot y}\frac{1}{\mathsf{fma}\left(\frac{1.128379167095512558560699289955664426088}{y}, e^{z}, -x\right)} + xdouble f(double x, double y, double z) {
double r350683 = x;
double r350684 = y;
double r350685 = 1.1283791670955126;
double r350686 = z;
double r350687 = exp(r350686);
double r350688 = r350685 * r350687;
double r350689 = r350683 * r350684;
double r350690 = r350688 - r350689;
double r350691 = r350684 / r350690;
double r350692 = r350683 + r350691;
return r350692;
}
double f(double x, double y, double z) {
double r350693 = 1.0;
double r350694 = 1.1283791670955126;
double r350695 = y;
double r350696 = r350694 / r350695;
double r350697 = z;
double r350698 = exp(r350697);
double r350699 = x;
double r350700 = -r350699;
double r350701 = fma(r350696, r350698, r350700);
double r350702 = r350693 / r350701;
double r350703 = r350702 + r350699;
return r350703;
}




Bits error versus x




Bits error versus y




Bits error versus z
| Original | 3.0 |
|---|---|
| Target | 0.1 |
| Herbie | 0.1 |
Initial program 3.0
rmApplied clear-num3.0
rmApplied add-cube-cbrt3.0
Applied associate-/l*3.0
Simplified0.1
Final simplification0.1
herbie shell --seed 2019209 +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.12837916709551256 y) (exp z)) x)))
(+ x (/ y (- (* 1.12837916709551256 (exp z)) (* x y)))))