\left(1 + t \cdot 2 \cdot 10^{-16}\right) \cdot \left(1 + t \cdot 2 \cdot 10^{-16}\right) + \left(-1 - 2 \cdot \left(t \cdot 2 \cdot 10^{-16}\right)\right)\left(3.9999999999999997 \cdot 10^{-32} \cdot t\right) \cdot tdouble f(double t) {
double r3262105 = 1.0;
double r3262106 = t;
double r3262107 = 2e-16;
double r3262108 = r3262106 * r3262107;
double r3262109 = r3262105 + r3262108;
double r3262110 = r3262109 * r3262109;
double r3262111 = -1.0;
double r3262112 = 2.0;
double r3262113 = r3262112 * r3262108;
double r3262114 = r3262111 - r3262113;
double r3262115 = r3262110 + r3262114;
return r3262115;
}
double f(double t) {
double r3262116 = 3.9999999999999997e-32;
double r3262117 = t;
double r3262118 = r3262116 * r3262117;
double r3262119 = r3262118 * r3262117;
return r3262119;
}




Bits error versus t
Results
| Original | 61.8 |
|---|---|
| Target | 50.6 |
| Herbie | 0.3 |
Initial program 61.8
Simplified0.4
rmApplied associate-*r*0.3
Taylor expanded around -inf 0.3
Final simplification0.3
herbie shell --seed 2019121
(FPCore (t)
:name "fma_test1"
:pre (<= 0.9 t 1.1)
:herbie-target
(fma (+ 1 (* t 2e-16)) (+ 1 (* t 2e-16)) (- -1 (* 2 (* t 2e-16))))
(+ (* (+ 1 (* t 2e-16)) (+ 1 (* t 2e-16))) (- -1 (* 2 (* t 2e-16)))))