\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\frac{1}{\left(x + 1\right) \cdot x} \cdot \frac{2}{x - 1}double f(double x) {
double r95842 = 1.0;
double r95843 = x;
double r95844 = r95843 + r95842;
double r95845 = r95842 / r95844;
double r95846 = 2.0;
double r95847 = r95846 / r95843;
double r95848 = r95845 - r95847;
double r95849 = r95843 - r95842;
double r95850 = r95842 / r95849;
double r95851 = r95848 + r95850;
return r95851;
}
double f(double x) {
double r95852 = 1.0;
double r95853 = x;
double r95854 = 1.0;
double r95855 = r95853 + r95854;
double r95856 = r95855 * r95853;
double r95857 = r95852 / r95856;
double r95858 = 2.0;
double r95859 = r95853 - r95854;
double r95860 = r95858 / r95859;
double r95861 = r95857 * r95860;
return r95861;
}




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