\frac{1 + \frac{2 \cdot t}{1 + t} \cdot \frac{2 \cdot t}{1 + t}}{2 + \frac{2 \cdot t}{1 + t} \cdot \frac{2 \cdot t}{1 + t}}\frac{1 + \frac{2 \cdot t}{1 + t} \cdot \frac{2 \cdot t}{1 + t}}{2 + \frac{2 \cdot t}{1 + t} \cdot \frac{2 \cdot t}{1 + t}}double f(double t) {
double r57338 = 1.0;
double r57339 = 2.0;
double r57340 = t;
double r57341 = r57339 * r57340;
double r57342 = r57338 + r57340;
double r57343 = r57341 / r57342;
double r57344 = r57343 * r57343;
double r57345 = r57338 + r57344;
double r57346 = r57339 + r57344;
double r57347 = r57345 / r57346;
return r57347;
}
double f(double t) {
double r57348 = 1.0;
double r57349 = 2.0;
double r57350 = t;
double r57351 = r57349 * r57350;
double r57352 = r57348 + r57350;
double r57353 = r57351 / r57352;
double r57354 = r57353 * r57353;
double r57355 = r57348 + r57354;
double r57356 = r57349 + r57354;
double r57357 = r57355 / r57356;
return r57357;
}



Bits error versus t
Results
Initial program 0.1
Final simplification0.1
herbie shell --seed 2020064 +o rules:numerics
(FPCore (t)
:name "Kahan p13 Example 1"
:precision binary64
(/ (+ 1 (* (/ (* 2 t) (+ 1 t)) (/ (* 2 t) (+ 1 t)))) (+ 2 (* (/ (* 2 t) (+ 1 t)) (/ (* 2 t) (+ 1 t))))))