\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\frac{1}{\left(x + 1\right) \cdot x} \cdot \frac{2}{x - 1}double f(double x) {
double r98538 = 1.0;
double r98539 = x;
double r98540 = r98539 + r98538;
double r98541 = r98538 / r98540;
double r98542 = 2.0;
double r98543 = r98542 / r98539;
double r98544 = r98541 - r98543;
double r98545 = r98539 - r98538;
double r98546 = r98538 / r98545;
double r98547 = r98544 + r98546;
return r98547;
}
double f(double x) {
double r98548 = 1.0;
double r98549 = x;
double r98550 = 1.0;
double r98551 = r98549 + r98550;
double r98552 = r98551 * r98549;
double r98553 = r98548 / r98552;
double r98554 = 2.0;
double r98555 = r98549 - r98550;
double r98556 = r98554 / r98555;
double r98557 = r98553 * r98556;
return r98557;
}




Bits error versus x
Results
| Original | 9.7 |
|---|---|
| Target | 0.3 |
| Herbie | 0.1 |
Initial program 9.7
rmApplied frac-sub25.8
Applied frac-add25.3
Simplified25.7
Taylor expanded around 0 0.3
rmApplied *-un-lft-identity0.3
Applied times-frac0.1
Final simplification0.1
herbie shell --seed 2019326 +o rules:numerics
(FPCore (x)
:name "3frac (problem 3.3.3)"
:precision binary64
:herbie-target
(/ 2 (* x (- (* x x) 1)))
(+ (- (/ 1 (+ x 1)) (/ 2 x)) (/ 1 (- x 1))))