\left(\frac{1}{x + 1} - \frac{2}{x}\right) + \frac{1}{x - 1}\frac{\frac{1}{x + 1}}{x} \cdot \frac{2}{x - 1}double f(double x) {
double r100500 = 1.0;
double r100501 = x;
double r100502 = r100501 + r100500;
double r100503 = r100500 / r100502;
double r100504 = 2.0;
double r100505 = r100504 / r100501;
double r100506 = r100503 - r100505;
double r100507 = r100501 - r100500;
double r100508 = r100500 / r100507;
double r100509 = r100506 + r100508;
return r100509;
}
double f(double x) {
double r100510 = 1.0;
double r100511 = x;
double r100512 = 1.0;
double r100513 = r100511 + r100512;
double r100514 = r100510 / r100513;
double r100515 = r100514 / r100511;
double r100516 = 2.0;
double r100517 = r100511 - r100512;
double r100518 = r100516 / r100517;
double r100519 = r100515 * r100518;
return r100519;
}




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