double code(double x) {
return ((double) (((double) (((double) (1.0 / ((double) (x + 1.0)))) - ((double) (2.0 / x)))) + ((double) (1.0 / ((double) (x - 1.0))))));
}
double code(double x) {
double VAR;
if (((x <= -838439.3309523697) || !(x <= 511.7848760495989))) {
VAR = ((double) (((double) (2.0 / ((double) pow(x, 7.0)))) + ((double) (((double) (2.0 / ((double) pow(x, 5.0)))) + ((double) (2.0 / ((double) pow(x, 3.0))))))));
} else {
VAR = ((double) (((double) (((double) (((double) (x - 1.0)) * ((double) (((double) (x * 1.0)) - ((double) (2.0 * ((double) (x + 1.0)))))))) + ((double) (1.0 * ((double) (x * ((double) (x + 1.0)))))))) / ((double) (x * ((double) (((double) (x * x)) - ((double) (1.0 * 1.0))))))));
}
return VAR;
}




Bits error versus x
Results
| Original | 9.5 |
|---|---|
| Target | 0.3 |
| Herbie | 0.3 |
if x < -838439.330952369724 or 511.78487604959889 < x Initial program 18.8
Taylor expanded around inf 0.5
Simplified0.5
if -838439.330952369724 < x < 511.78487604959889Initial program 0.2
rmApplied frac-sub0.2
Applied frac-add0.0
Simplified0.0
Simplified0.0
Final simplification0.3
herbie shell --seed 2020181
(FPCore (x)
:name "3frac (problem 3.3.3)"
:precision binary64
:herbie-target
(/ 2.0 (* x (- (* x x) 1.0)))
(+ (- (/ 1.0 (+ x 1.0)) (/ 2.0 x)) (/ 1.0 (- x 1.0))))