double code(double x, double y, double z) {
return ((double) (((double) (((double) (((double) (((double) (x - 0.5)) * ((double) log(x)))) - x)) + 0.91893853320467)) + ((double) (((double) (((double) (((double) (((double) (((double) (y + 0.0007936500793651)) * z)) - 0.0027777777777778)) * z)) + 0.083333333333333)) / x))));
}
double code(double x, double y, double z) {
double VAR;
if (((z <= -1.8860184916081755e+24) || !(z <= 3.529602132269039e+16))) {
VAR = ((double) (((double) (((double) (((double) (((double) (x - 0.5)) * ((double) log(x)))) - x)) + 0.91893853320467)) + ((double) (((double) (((double) (z * ((double) (z / x)))) * ((double) (y + 0.0007936500793651)))) - ((double) (z * ((double) (0.0027777777777778 / x))))))));
} else {
VAR = ((double) (((double) (0.91893853320467 + ((double) (((double) (((double) (((double) (x - 0.5)) * ((double) (((double) log(((double) pow(x, 0.3333333333333333)))) * 2.0)))) + ((double) (((double) (x - 0.5)) * ((double) log(((double) cbrt(x)))))))) - x)))) + ((double) (((double) (((double) (z * ((double) (((double) (z * ((double) (y + 0.0007936500793651)))) - 0.0027777777777778)))) + 0.083333333333333)) / x))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.0 |
|---|---|
| Target | 1.3 |
| Herbie | 0.4 |
if z < -1.88601849160817548e24 or 35296021322690392 < z Initial program 22.7
Taylor expanded around inf 23.1
Simplified0.7
if -1.88601849160817548e24 < z < 35296021322690392Initial program 0.3
rmApplied add-cube-cbrt0.3
Applied log-prod0.3
Applied distribute-lft-in0.3
Simplified0.3
rmApplied pow1/30.3
Final simplification0.4
herbie shell --seed 2020180
(FPCore (x y z)
:name "Numeric.SpecFunctions:$slogFactorial from math-functions-0.1.5.2, B"
:precision binary64
:herbie-target
(+ (+ (+ (* (- x 0.5) (log x)) (- 0.91893853320467 x)) (/ 0.083333333333333 x)) (* (/ z x) (- (* z (+ y 0.0007936500793651)) 0.0027777777777778)))
(+ (+ (- (* (- x 0.5) (log x)) x) 0.91893853320467) (/ (+ (* (- (* (+ y 0.0007936500793651) z) 0.0027777777777778) z) 0.083333333333333) x)))