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) {
double VAR;
if (((y <= 710.9535870360606) || !((y <= 5.967650503494568e+51) || (!(y <= 3.164678549410275e+86) && (y <= 1.1647152521001902e+159))))) {
VAR = ((double) (((double) (((double) exp(((double) (x * ((double) (((double) log(((double) (((double) cbrt(x)) / ((double) cbrt(((double) (y + x)))))))) * 2.0)))))) * ((double) pow(((double) (((double) cbrt(x)) / ((double) cbrt(((double) (y + x)))))), x)))) / x));
} else {
VAR = ((double) log(((double) exp(((double) (((double) pow(((double) (x / ((double) (y + x)))), x)) / x))))));
}
return VAR;
}




Bits error versus x




Bits error versus y
Results
| Original | 11.0 |
|---|---|
| Target | 8.3 |
| Herbie | 4.9 |
if y < 710.95358703606064 or 5.96765050349456811e51 < y < 3.16467854941027495e86 or 1.16471525210019018e159 < y Initial program 8.3
Simplified8.3
rmApplied add-cube-cbrt28.4
Applied add-cube-cbrt8.3
Applied times-frac8.3
Applied unpow-prod-down4.3
Simplified4.3
rmApplied add-exp-log46.4
Applied add-exp-log34.2
Applied div-exp34.2
Applied add-exp-log46.4
Applied add-exp-log34.2
Applied div-exp34.2
Applied prod-exp34.2
Applied pow-exp33.6
Simplified3.0
if 710.95358703606064 < y < 5.96765050349456811e51 or 3.16467854941027495e86 < y < 1.16471525210019018e159Initial program 36.0
Simplified36.0
rmApplied add-log-exp22.5
Final simplification4.9
herbie shell --seed 2020180
(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))