\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 r66897 = 1.0;
double r66898 = x;
double r66899 = r66898 + r66897;
double r66900 = r66897 / r66899;
double r66901 = 2.0;
double r66902 = r66901 / r66898;
double r66903 = r66900 - r66902;
double r66904 = r66898 - r66897;
double r66905 = r66897 / r66904;
double r66906 = r66903 + r66905;
return r66906;
}
double f(double x) {
double r66907 = 2.0;
double r66908 = x;
double r66909 = 1.0;
double r66910 = r66908 + r66909;
double r66911 = r66910 * r66908;
double r66912 = r66907 / r66911;
double r66913 = r66908 - r66909;
double r66914 = r66912 / r66913;
return r66914;
}




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