\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 r75031 = 1.0;
double r75032 = 2.0;
double r75033 = t;
double r75034 = r75032 * r75033;
double r75035 = r75031 + r75033;
double r75036 = r75034 / r75035;
double r75037 = r75036 * r75036;
double r75038 = r75031 + r75037;
double r75039 = r75032 + r75037;
double r75040 = r75038 / r75039;
return r75040;
}
double f(double t) {
double r75041 = 1.0;
double r75042 = 2.0;
double r75043 = t;
double r75044 = r75042 * r75043;
double r75045 = r75041 + r75043;
double r75046 = r75044 / r75045;
double r75047 = r75046 * r75046;
double r75048 = r75041 + r75047;
double r75049 = r75042 + r75047;
double r75050 = r75048 / r75049;
return r75050;
}



Bits error versus t
Results
Initial program 0.1
Final simplification0.1
herbie shell --seed 2020001
(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))))))