double code(double x, double y, double z, double t, double a, double b) {
return ((double) (x + ((double) (((double) (y * ((double) (((double) (((double) (((double) (((double) (((double) (((double) (((double) (z * 3.13060547623)) + 11.1667541262)) * z)) + t)) * z)) + a)) * z)) + b)))) / ((double) (((double) (((double) (((double) (((double) (((double) (((double) (z + 15.234687407)) * z)) + 31.4690115749)) * z)) + 11.9400905721)) * z)) + 0.607771387771))))));
}
double code(double x, double y, double z, double t, double a, double b) {
double VAR;
if (((z <= -1.0110210441160586e+36) || !(z <= 3.4849124286894035e+37))) {
VAR = ((double) (x + ((double) (y * ((double) (3.13060547623 + ((double) (((double) (t / ((double) (z * z)))) - ((double) (36.527041698806414 / z))))))))));
} else {
VAR = ((double) (x + ((double) (y * ((double) (((double) (1.0 / ((double) sqrt(((double) (((double) (z * ((double) (((double) (z * ((double) (((double) (z * ((double) (z + 15.234687407)))) + 31.4690115749)))) + 11.9400905721)))) + 0.607771387771)))))) * ((double) (((double) (((double) (z * ((double) (((double) (z * ((double) (t + ((double) (z * ((double) (((double) (z * 3.13060547623)) + 11.1667541262)))))))) + a)))) + b)) / ((double) sqrt(((double) (((double) (z * ((double) (((double) (z * ((double) (((double) (z * ((double) (z + 15.234687407)))) + 31.4690115749)))) + 11.9400905721)))) + 0.607771387771))))))))))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b
Results
| Original | 29.4 |
|---|---|
| Target | 1.0 |
| Herbie | 1.2 |
if z < -1.01102104411605861e36 or 3.4849124286894035e37 < z Initial program 59.8
Simplified57.2
Taylor expanded around inf 1.6
Simplified1.6
if -1.01102104411605861e36 < z < 3.4849124286894035e37Initial program 1.3
Simplified0.6
rmApplied add-sqr-sqrt1.0
Applied *-un-lft-identity1.0
Applied times-frac0.9
Final simplification1.2
herbie shell --seed 2020180
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2, D"
:precision binary64
:herbie-target
(if (< z -6.499344996252632e+53) (+ x (* (+ (- 3.13060547623 (/ 36.527041698806414 z)) (/ t (* z z))) (/ y 1.0))) (if (< z 7.066965436914287e+59) (+ x (/ y (/ (+ (* (+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721) z) 0.607771387771) (+ (* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z) b)))) (+ x (* (+ (- 3.13060547623 (/ 36.527041698806414 z)) (/ t (* z z))) (/ y 1.0)))))
(+ x (/ (* y (+ (* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z) b)) (+ (* (+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721) z) 0.607771387771))))