\frac{e^{x \cdot \log \left(\frac{x}{x + y}\right)}}{x}\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}double code(double x, double y) {
return ((double) (((double) exp(((double) (x * ((double) log(((double) (x / ((double) (x + y)))))))))) / x));
}
double code(double x, double y) {
return ((double) (((double) (((double) (((double) pow(((double) (((double) cbrt(x)) / ((double) cbrt(((double) (x + y)))))), x)) * ((double) pow(((double) (((double) cbrt(x)) / ((double) cbrt(((double) (x + y)))))), x)))) * ((double) pow(((double) (((double) cbrt(x)) / ((double) cbrt(((double) (x + y)))))), x)))) / x));
}




Bits error versus x




Bits error versus y
Results
| Original | 10.7 |
|---|---|
| Target | 7.8 |
| Herbie | 5.6 |
Initial program 10.7
Simplified10.7
rmApplied add-cube-cbrt27.0
Applied add-cube-cbrt10.7
Applied times-frac10.7
Applied unpow-prod-down6.7
rmApplied times-frac6.7
Applied unpow-prod-down5.6
Final simplification5.6
herbie shell --seed 2020155
(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.0 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.0 y)) x))))
(/ (exp (* x (log (/ x (+ x y))))) x))