\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\frac{\frac{2}{\left(x + 1\right) \cdot x}}{x - 1}double f(double x) {
double r4532490 = 1.0;
double r4532491 = x;
double r4532492 = r4532491 + r4532490;
double r4532493 = r4532490 / r4532492;
double r4532494 = 2.0;
double r4532495 = r4532494 / r4532491;
double r4532496 = r4532493 - r4532495;
double r4532497 = r4532491 - r4532490;
double r4532498 = r4532490 / r4532497;
double r4532499 = r4532496 + r4532498;
return r4532499;
}
double f(double x) {
double r4532500 = 2.0;
double r4532501 = x;
double r4532502 = 1.0;
double r4532503 = r4532501 + r4532502;
double r4532504 = r4532503 * r4532501;
double r4532505 = r4532500 / r4532504;
double r4532506 = r4532501 - r4532502;
double r4532507 = r4532505 / r4532506;
return r4532507;
}




Bits error versus x
Results
| Original | 9.9 |
|---|---|
| Target | 0.3 |
| Herbie | 0.1 |
Initial program 9.9
rmApplied frac-sub25.7
Applied frac-add25.3
Taylor expanded around 0 0.3
rmApplied associate-/r*0.1
Final simplification0.1
herbie shell --seed 2019174 +o rules:numerics
(FPCore (x)
:name "3frac (problem 3.3.3)"
:herbie-target
(/ 2.0 (* x (- (* x x) 1.0)))
(+ (- (/ 1.0 (+ x 1.0)) (/ 2.0 x)) (/ 1.0 (- x 1.0))))