\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 r28344 = 1.0;
double r28345 = 2.0;
double r28346 = t;
double r28347 = r28345 / r28346;
double r28348 = r28344 / r28346;
double r28349 = r28344 + r28348;
double r28350 = r28347 / r28349;
double r28351 = r28345 - r28350;
double r28352 = r28351 * r28351;
double r28353 = r28344 + r28352;
double r28354 = r28345 + r28352;
double r28355 = r28353 / r28354;
return r28355;
}
double f(double t) {
double r28356 = 1.0;
double r28357 = 2.0;
double r28358 = t;
double r28359 = r28358 * r28356;
double r28360 = r28356 + r28359;
double r28361 = r28357 / r28360;
double r28362 = r28357 - r28361;
double r28363 = r28362 * r28362;
double r28364 = r28356 + r28363;
double r28365 = r28357 + r28363;
double r28366 = r28364 / r28365;
return r28366;
}



Bits error versus t
Results
Initial program 0.0
Simplified0.0
Final simplification0.0
herbie shell --seed 2019322
(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))))))))