\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 r101867 = 1.0;
double r101868 = x;
double r101869 = r101868 + r101867;
double r101870 = r101867 / r101869;
double r101871 = 2.0;
double r101872 = r101871 / r101868;
double r101873 = r101870 - r101872;
double r101874 = r101868 - r101867;
double r101875 = r101867 / r101874;
double r101876 = r101873 + r101875;
return r101876;
}
double f(double x) {
double r101877 = 1.0;
double r101878 = x;
double r101879 = 1.0;
double r101880 = r101878 + r101879;
double r101881 = r101877 / r101880;
double r101882 = r101881 / r101878;
double r101883 = 2.0;
double r101884 = r101878 - r101879;
double r101885 = r101883 / r101884;
double r101886 = r101882 * r101885;
return r101886;
}




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))))