\frac{e^{x \cdot \log \left(\frac{x}{x + y}\right)}}{x}\begin{array}{l}
\mathbf{if}\;x \le -2.575795415940148063364153542342468795089 \cdot 10^{92} \lor \neg \left(x \le 2.134762495789531033807406856794841587543\right):\\
\;\;\;\;\frac{e^{-y}}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{e^{x \cdot \left(\left(\log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{x + y}}\right) + \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{x + y}}\right)\right) + \log \left(\frac{\sqrt[3]{x}}{\sqrt[3]{x + y}}\right)\right)}}{x}\\
\end{array}double f(double x, double y) {
double r274618 = x;
double r274619 = y;
double r274620 = r274618 + r274619;
double r274621 = r274618 / r274620;
double r274622 = log(r274621);
double r274623 = r274618 * r274622;
double r274624 = exp(r274623);
double r274625 = r274624 / r274618;
return r274625;
}
double f(double x, double y) {
double r274626 = x;
double r274627 = -2.575795415940148e+92;
bool r274628 = r274626 <= r274627;
double r274629 = 2.134762495789531;
bool r274630 = r274626 <= r274629;
double r274631 = !r274630;
bool r274632 = r274628 || r274631;
double r274633 = y;
double r274634 = -r274633;
double r274635 = exp(r274634);
double r274636 = r274635 / r274626;
double r274637 = cbrt(r274626);
double r274638 = r274626 + r274633;
double r274639 = cbrt(r274638);
double r274640 = r274637 / r274639;
double r274641 = log(r274640);
double r274642 = r274641 + r274641;
double r274643 = r274642 + r274641;
double r274644 = r274626 * r274643;
double r274645 = exp(r274644);
double r274646 = r274645 / r274626;
double r274647 = r274632 ? r274636 : r274646;
return r274647;
}




Bits error versus x




Bits error versus y
Results
| Original | 11.0 |
|---|---|
| Target | 7.5 |
| Herbie | 0.2 |
if x < -2.575795415940148e+92 or 2.134762495789531 < x Initial program 11.1
Taylor expanded around inf 0.0
Simplified0.0
if -2.575795415940148e+92 < x < 2.134762495789531Initial program 10.8
rmApplied add-cube-cbrt15.9
Applied add-cube-cbrt10.9
Applied times-frac10.9
Applied log-prod2.4
rmApplied times-frac2.4
Applied log-prod0.3
Final simplification0.2
herbie shell --seed 2019208 +o rules:numerics
(FPCore (x y)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, F"
:precision binary64
:herbie-target
(if (< y -3.73118442066479561e94) (/ (exp (/ -1 y)) x) (if (< y 2.81795924272828789e37) (/ (pow (/ x (+ y x)) x) x) (if (< y 2.347387415166998e178) (log (exp (/ (pow (/ x (+ y x)) x) x))) (/ (exp (/ -1 y)) x))))
(/ (exp (* x (log (/ x (+ x y))))) x))