\frac{x + \frac{y \cdot z - x}{t \cdot z - x}}{x + 1.0}\begin{array}{l}
\mathbf{if}\;z \le -1.1409575287436362 \cdot 10^{+77}:\\
\;\;\;\;\frac{x + \frac{y}{t}}{x + 1.0}\\
\mathbf{elif}\;z \le 1.2840647236389548 \cdot 10^{+57}:\\
\;\;\;\;\frac{x + \frac{1}{\frac{t \cdot z - x}{y \cdot z - x}}}{x + 1.0}\\
\mathbf{else}:\\
\;\;\;\;\frac{x + \frac{y}{t}}{x + 1.0}\\
\end{array}double f(double x, double y, double z, double t) {
double r34567985 = x;
double r34567986 = y;
double r34567987 = z;
double r34567988 = r34567986 * r34567987;
double r34567989 = r34567988 - r34567985;
double r34567990 = t;
double r34567991 = r34567990 * r34567987;
double r34567992 = r34567991 - r34567985;
double r34567993 = r34567989 / r34567992;
double r34567994 = r34567985 + r34567993;
double r34567995 = 1.0;
double r34567996 = r34567985 + r34567995;
double r34567997 = r34567994 / r34567996;
return r34567997;
}
double f(double x, double y, double z, double t) {
double r34567998 = z;
double r34567999 = -1.1409575287436362e+77;
bool r34568000 = r34567998 <= r34567999;
double r34568001 = x;
double r34568002 = y;
double r34568003 = t;
double r34568004 = r34568002 / r34568003;
double r34568005 = r34568001 + r34568004;
double r34568006 = 1.0;
double r34568007 = r34568001 + r34568006;
double r34568008 = r34568005 / r34568007;
double r34568009 = 1.2840647236389548e+57;
bool r34568010 = r34567998 <= r34568009;
double r34568011 = 1.0;
double r34568012 = r34568003 * r34567998;
double r34568013 = r34568012 - r34568001;
double r34568014 = r34568002 * r34567998;
double r34568015 = r34568014 - r34568001;
double r34568016 = r34568013 / r34568015;
double r34568017 = r34568011 / r34568016;
double r34568018 = r34568001 + r34568017;
double r34568019 = r34568018 / r34568007;
double r34568020 = r34568010 ? r34568019 : r34568008;
double r34568021 = r34568000 ? r34568008 : r34568020;
return r34568021;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 7.2 |
|---|---|
| Target | 0.4 |
| Herbie | 3.4 |
if z < -1.1409575287436362e+77 or 1.2840647236389548e+57 < z Initial program 18.0
Taylor expanded around inf 7.7
if -1.1409575287436362e+77 < z < 1.2840647236389548e+57Initial program 0.7
rmApplied clear-num0.8
rmApplied *-un-lft-identity0.8
Applied associate-/r*0.8
Simplified0.8
Final simplification3.4
herbie shell --seed 2019163
(FPCore (x y z t)
:name "Diagrams.Trail:splitAtParam from diagrams-lib-1.3.0.3, A"
:herbie-target
(/ (+ x (- (/ y (- t (/ x z))) (/ x (- (* t z) x)))) (+ x 1.0))
(/ (+ x (/ (- (* y z) x) (- (* t z) x))) (+ x 1.0)))