\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\frac{2}{x - 1} \cdot \frac{1}{\left(x + 1\right) \cdot x}double f(double x) {
double r7647671 = 1.0;
double r7647672 = x;
double r7647673 = r7647672 + r7647671;
double r7647674 = r7647671 / r7647673;
double r7647675 = 2.0;
double r7647676 = r7647675 / r7647672;
double r7647677 = r7647674 - r7647676;
double r7647678 = r7647672 - r7647671;
double r7647679 = r7647671 / r7647678;
double r7647680 = r7647677 + r7647679;
return r7647680;
}
double f(double x) {
double r7647681 = 2.0;
double r7647682 = x;
double r7647683 = 1.0;
double r7647684 = r7647682 - r7647683;
double r7647685 = r7647681 / r7647684;
double r7647686 = 1.0;
double r7647687 = r7647682 + r7647683;
double r7647688 = r7647687 * r7647682;
double r7647689 = r7647686 / r7647688;
double r7647690 = r7647685 * r7647689;
return r7647690;
}




Bits error versus x
Results
| Original | 9.7 |
|---|---|
| Target | 0.2 |
| Herbie | 0.1 |
Initial program 9.7
rmApplied frac-sub26.2
Applied frac-add25.7
Taylor expanded around 0 0.2
rmApplied *-un-lft-identity0.2
Applied times-frac0.1
Final simplification0.1
herbie shell --seed 2019171
(FPCore (x)
:name "3frac (problem 3.3.3)"
:herbie-target
(/ 2.0 (* x (- (* x x) 1.0)))
(+ (- (/ 1.0 (+ x 1.0)) (/ 2.0 x)) (/ 1.0 (- x 1.0))))