\frac{e^{x \cdot \log \left(\frac{x}{x + y}\right)}}{x}\frac{{\left(\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{\sqrt[3]{x + y} \cdot \sqrt[3]{x + y}}\right)}^{x} \cdot {\left(\frac{\sqrt[3]{x}}{\sqrt[3]{x + y}}\right)}^{x}}{x}double code(double x, double y) {
return (((double) exp(((double) (x * ((double) log((x / ((double) (x + y))))))))) / x);
}
double code(double x, double y) {
return (((double) (((double) pow((((double) (((double) cbrt(x)) * ((double) cbrt(x)))) / ((double) (((double) cbrt(((double) (x + y)))) * ((double) cbrt(((double) (x + y))))))), x)) * ((double) pow((((double) cbrt(x)) / ((double) cbrt(((double) (x + y))))), x)))) / x);
}




Bits error versus x




Bits error versus y
Results
| Original | 11.1 |
|---|---|
| Target | 8.2 |
| Herbie | 6.7 |
Initial program 11.1
Simplified11.1
rmApplied add-cube-cbrt28.2
Applied add-cube-cbrt11.1
Applied times-frac11.1
Applied unpow-prod-down6.7
Final simplification6.7
herbie shell --seed 2020182
(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))