\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\frac{\frac{1}{x + 1}}{x} \cdot \frac{2}{x - 1}double f(double x) {
double r77562 = 1.0;
double r77563 = x;
double r77564 = r77563 + r77562;
double r77565 = r77562 / r77564;
double r77566 = 2.0;
double r77567 = r77566 / r77563;
double r77568 = r77565 - r77567;
double r77569 = r77563 - r77562;
double r77570 = r77562 / r77569;
double r77571 = r77568 + r77570;
return r77571;
}
double f(double x) {
double r77572 = 1.0;
double r77573 = x;
double r77574 = 1.0;
double r77575 = r77573 + r77574;
double r77576 = r77572 / r77575;
double r77577 = r77576 / r77573;
double r77578 = 2.0;
double r77579 = r77573 - r77574;
double r77580 = r77578 / r77579;
double r77581 = r77577 * r77580;
return r77581;
}




Bits error versus x
Results
| Original | 9.9 |
|---|---|
| Target | 0.3 |
| Herbie | 0.1 |
Initial program 9.9
rmApplied frac-sub26.1
Applied frac-add25.8
Simplified26.1
Taylor expanded around 0 0.3
rmApplied *-un-lft-identity0.3
Applied times-frac0.1
rmApplied associate-/r*0.1
Final simplification0.1
herbie shell --seed 2019323 +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))))