\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 r139662 = 1.0;
double r139663 = x;
double r139664 = r139663 + r139662;
double r139665 = r139662 / r139664;
double r139666 = 2.0;
double r139667 = r139666 / r139663;
double r139668 = r139665 - r139667;
double r139669 = r139663 - r139662;
double r139670 = r139662 / r139669;
double r139671 = r139668 + r139670;
return r139671;
}
double f(double x) {
double r139672 = 1.0;
double r139673 = x;
double r139674 = 1.0;
double r139675 = r139673 + r139674;
double r139676 = r139675 * r139673;
double r139677 = r139672 / r139676;
double r139678 = 2.0;
double r139679 = r139673 - r139674;
double r139680 = r139678 / r139679;
double r139681 = r139677 * r139680;
return r139681;
}




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