\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 r162562 = 1.0;
double r162563 = x;
double r162564 = r162563 + r162562;
double r162565 = r162562 / r162564;
double r162566 = 2.0;
double r162567 = r162566 / r162563;
double r162568 = r162565 - r162567;
double r162569 = r162563 - r162562;
double r162570 = r162562 / r162569;
double r162571 = r162568 + r162570;
return r162571;
}
double f(double x) {
double r162572 = 1.0;
double r162573 = x;
double r162574 = 1.0;
double r162575 = r162573 + r162574;
double r162576 = r162575 * r162573;
double r162577 = r162572 / r162576;
double r162578 = 2.0;
double r162579 = r162573 - r162574;
double r162580 = r162578 / r162579;
double r162581 = r162577 * r162580;
return r162581;
}




Bits error versus x
Results
| Original | 10.6 |
|---|---|
| Target | 0.3 |
| Herbie | 0.1 |
Initial program 10.6
rmApplied frac-sub26.2
Applied frac-add25.6
Taylor expanded around 0 0.3
rmApplied *-un-lft-identity0.3
Applied times-frac0.1
Final simplification0.1
herbie shell --seed 2020045
(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))))