\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 r2944971 = 1.0;
double r2944972 = t;
double r2944973 = 2e-16;
double r2944974 = r2944972 * r2944973;
double r2944975 = r2944971 + r2944974;
double r2944976 = r2944975 * r2944975;
double r2944977 = -1.0;
double r2944978 = 2.0;
double r2944979 = r2944978 * r2944974;
double r2944980 = r2944977 - r2944979;
double r2944981 = r2944976 + r2944980;
return r2944981;
}
double f(double t) {
double r2944982 = 3.9999999999999997e-32;
double r2944983 = t;
double r2944984 = r2944982 * r2944983;
double r2944985 = r2944984 * r2944983;
return r2944985;
}




Bits error versus t
Results
| Original | 61.8 |
|---|---|
| Target | 50.5 |
| 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 2019133
(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)))))