\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{\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)}double f(double t) {
double r28869 = 1.0;
double r28870 = 2.0;
double r28871 = t;
double r28872 = r28870 / r28871;
double r28873 = r28869 / r28871;
double r28874 = r28869 + r28873;
double r28875 = r28872 / r28874;
double r28876 = r28870 - r28875;
double r28877 = r28876 * r28876;
double r28878 = r28869 + r28877;
double r28879 = r28870 + r28877;
double r28880 = r28878 / r28879;
return r28880;
}
double f(double t) {
double r28881 = 1.0;
double r28882 = 2.0;
double r28883 = t;
double r28884 = r28882 / r28883;
double r28885 = r28881 / r28883;
double r28886 = r28881 + r28885;
double r28887 = r28884 / r28886;
double r28888 = r28882 - r28887;
double r28889 = r28888 * r28888;
double r28890 = r28881 + r28889;
double r28891 = r28882 + r28889;
double r28892 = r28890 / r28891;
return r28892;
}



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