\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 r162066 = 1.0;
double r162067 = x;
double r162068 = r162067 + r162066;
double r162069 = r162066 / r162068;
double r162070 = 2.0;
double r162071 = r162070 / r162067;
double r162072 = r162069 - r162071;
double r162073 = r162067 - r162066;
double r162074 = r162066 / r162073;
double r162075 = r162072 + r162074;
return r162075;
}
double f(double x) {
double r162076 = 1.0;
double r162077 = x;
double r162078 = 1.0;
double r162079 = r162077 + r162078;
double r162080 = r162079 * r162077;
double r162081 = r162076 / r162080;
double r162082 = 2.0;
double r162083 = r162077 - r162078;
double r162084 = r162082 / r162083;
double r162085 = r162081 * r162084;
return r162085;
}




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