\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 r133892 = 1.0;
double r133893 = x;
double r133894 = r133893 + r133892;
double r133895 = r133892 / r133894;
double r133896 = 2.0;
double r133897 = r133896 / r133893;
double r133898 = r133895 - r133897;
double r133899 = r133893 - r133892;
double r133900 = r133892 / r133899;
double r133901 = r133898 + r133900;
return r133901;
}
double f(double x) {
double r133902 = 2.0;
double r133903 = x;
double r133904 = 1.0;
double r133905 = r133903 + r133904;
double r133906 = r133905 * r133903;
double r133907 = r133902 / r133906;
double r133908 = r133903 - r133904;
double r133909 = r133907 / r133908;
return r133909;
}




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