\frac{1 + \left(2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\right) \cdot \left(2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\right)}{2 + \left(2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\right) \cdot \left(2 - \frac{\frac{2}{t}}{1 + \frac{1}{t}}\right)}\frac{1 + \left(2 - \frac{2}{1 + t \cdot 1}\right) \cdot \left(2 - \frac{2}{1 + t \cdot 1}\right)}{2 + \left(2 - \frac{2}{1 + t \cdot 1}\right) \cdot \left(2 - \frac{2}{1 + t \cdot 1}\right)}double f(double t) {
double r30583 = 1.0;
double r30584 = 2.0;
double r30585 = t;
double r30586 = r30584 / r30585;
double r30587 = r30583 / r30585;
double r30588 = r30583 + r30587;
double r30589 = r30586 / r30588;
double r30590 = r30584 - r30589;
double r30591 = r30590 * r30590;
double r30592 = r30583 + r30591;
double r30593 = r30584 + r30591;
double r30594 = r30592 / r30593;
return r30594;
}
double f(double t) {
double r30595 = 1.0;
double r30596 = 2.0;
double r30597 = t;
double r30598 = r30597 * r30595;
double r30599 = r30595 + r30598;
double r30600 = r30596 / r30599;
double r30601 = r30596 - r30600;
double r30602 = r30601 * r30601;
double r30603 = r30595 + r30602;
double r30604 = r30596 + r30602;
double r30605 = r30603 / r30604;
return r30605;
}



Bits error versus t
Results
Initial program 0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019326
(FPCore (t)
:name "Kahan p13 Example 2"
:precision binary64
(/ (+ 1 (* (- 2 (/ (/ 2 t) (+ 1 (/ 1 t)))) (- 2 (/ (/ 2 t) (+ 1 (/ 1 t)))))) (+ 2 (* (- 2 (/ (/ 2 t) (+ 1 (/ 1 t)))) (- 2 (/ (/ 2 t) (+ 1 (/ 1 t))))))))