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 ((x <= 477006.7565026866)) {
VAR = ((double) (((double) (((double) (x - 0.5)) * ((double) log(((double) sqrt(x)))))) + ((double) (((double) (((double) (((double) (z * ((double) (((double) (z * ((double) (y + 0.0007936500793651)))) - 0.0027777777777778)))) + 0.083333333333333)) / x)) + ((double) (((double) (((double) (x - 0.5)) * ((double) log(((double) sqrt(x)))))) + ((double) (0.91893853320467 - x))))))));
} else {
VAR = ((double) (((double) (((double) (x - 0.5)) * ((double) log(((double) (((double) cbrt(x)) * ((double) cbrt(x)))))))) + ((double) (((double) (((double) (((double) (y + 0.0007936500793651)) * ((double) (z * ((double) (z / x)))))) - ((double) (0.0027777777777778 * ((double) (z / x)))))) + ((double) (((double) (0.91893853320467 - x)) + ((double) (((double) (x - 0.5)) * ((double) log(((double) cbrt(x))))))))))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.0 |
|---|---|
| Target | 1.2 |
| Herbie | 0.4 |
if x < 477006.756502686592Initial program 0.1
Simplified0.1
rmApplied add-sqr-sqrt0.1
Applied log-prod0.1
Applied distribute-lft-in0.1
Applied associate-+l+0.1
Simplified0.1
if 477006.756502686592 < x Initial program 10.0
Simplified10.0
rmApplied add-cube-cbrt10.0
Applied log-prod10.0
Applied distribute-lft-in10.0
Applied associate-+l+10.0
Simplified10.0
Taylor expanded around inf 10.1
Simplified0.5
Final simplification0.4
herbie shell --seed 2020190
(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)))