\frac{e^{x \cdot \log \left(\frac{x}{x + y}\right)}}{x}\begin{array}{l}
\mathbf{if}\;x \le -3.371891914955366 \cdot 10^{41} \lor \neg \left(x \le 6.5526303521982673 \cdot 10^{-6}\right):\\
\;\;\;\;\frac{e^{-y}}{x}\\
\mathbf{else}:\\
\;\;\;\;\frac{\left({\left(\frac{\sqrt[3]{x}}{\sqrt[3]{x + y}}\right)}^{x} \cdot {\left(\frac{\sqrt[3]{x}}{\sqrt[3]{x + y}}\right)}^{x}\right) \cdot {\left(\frac{\sqrt[3]{x}}{\sqrt[3]{x + y}}\right)}^{x}}{x}\\
\end{array}double f(double x, double y) {
double r370633 = x;
double r370634 = y;
double r370635 = r370633 + r370634;
double r370636 = r370633 / r370635;
double r370637 = log(r370636);
double r370638 = r370633 * r370637;
double r370639 = exp(r370638);
double r370640 = r370639 / r370633;
return r370640;
}
double f(double x, double y) {
double r370641 = x;
double r370642 = -3.371891914955366e+41;
bool r370643 = r370641 <= r370642;
double r370644 = 6.552630352198267e-06;
bool r370645 = r370641 <= r370644;
double r370646 = !r370645;
bool r370647 = r370643 || r370646;
double r370648 = y;
double r370649 = -r370648;
double r370650 = exp(r370649);
double r370651 = r370650 / r370641;
double r370652 = cbrt(r370641);
double r370653 = r370641 + r370648;
double r370654 = cbrt(r370653);
double r370655 = r370652 / r370654;
double r370656 = pow(r370655, r370641);
double r370657 = r370656 * r370656;
double r370658 = r370657 * r370656;
double r370659 = r370658 / r370641;
double r370660 = r370647 ? r370651 : r370659;
return r370660;
}




Bits error versus x




Bits error versus y
Results
| Original | 10.9 |
|---|---|
| Target | 8.2 |
| Herbie | 0.2 |
if x < -3.371891914955366e+41 or 6.552630352198267e-06 < x Initial program 11.0
Simplified11.0
Taylor expanded around inf 0.3
if -3.371891914955366e+41 < x < 6.552630352198267e-06Initial program 10.8
Simplified10.7
rmApplied add-cube-cbrt12.7
Applied add-cube-cbrt10.8
Applied times-frac10.8
Applied unpow-prod-down2.4
rmApplied times-frac2.4
Applied unpow-prod-down0.2
Final simplification0.2
herbie shell --seed 2020047 +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.7311844206647956e+94) (/ (exp (/ -1 y)) x) (if (< y 2.817959242728288e+37) (/ (pow (/ x (+ y x)) x) x) (if (< y 2.347387415166998e+178) (log (exp (/ (pow (/ x (+ y x)) x) x))) (/ (exp (/ -1 y)) x))))
(/ (exp (* x (log (/ x (+ x y))))) x))