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




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.1 |
|---|---|
| Target | 1.3 |
| Herbie | 0.3 |
if x < 5884704.858232215Initial program 0.1
Simplified0.1
rmApplied add-sqr-sqrt0.1
Applied add-sqr-sqrt0.1
Applied difference-of-squares0.1
Applied associate-*l*0.1
Simplified0.1
if 5884704.858232215 < x Initial program 10.4
Simplified10.4
Taylor expanded around inf 10.5
Simplified0.4
Final simplification0.3
herbie shell --seed 2020191
(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)))