\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)t \cdot \left(t \cdot 3.9999999999999997 \cdot 10^{-32}\right)double f(double t) {
double r3014250 = 1.0;
double r3014251 = t;
double r3014252 = 2e-16;
double r3014253 = r3014251 * r3014252;
double r3014254 = r3014250 + r3014253;
double r3014255 = r3014254 * r3014254;
double r3014256 = -1.0;
double r3014257 = 2.0;
double r3014258 = r3014257 * r3014253;
double r3014259 = r3014256 - r3014258;
double r3014260 = r3014255 + r3014259;
return r3014260;
}
double f(double t) {
double r3014261 = t;
double r3014262 = 3.9999999999999997e-32;
double r3014263 = r3014261 * r3014262;
double r3014264 = r3014261 * r3014263;
return r3014264;
}




Bits error versus t
Results
| Original | 61.8 |
|---|---|
| Target | 50.6 |
| Herbie | 0.3 |
Initial program 61.8
Simplified0.4
rmApplied associate-*l*0.3
Taylor expanded around 0 0.3
Final simplification0.3
herbie shell --seed 2019162 +o rules:numerics
(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)))))