\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)3.9999999999999997 \cdot 10^{-32} \cdot {t}^{2}double code(double t) {
return ((double) (((double) (((double) (1.0 + ((double) (t * 2e-16)))) * ((double) (1.0 + ((double) (t * 2e-16)))))) + ((double) (-1.0 - ((double) (2.0 * ((double) (t * 2e-16))))))));
}
double code(double t) {
return ((double) (3.9999999999999997e-32 * ((double) pow(t, 2.0))));
}




Bits error versus t
Results
| Original | 61.8 |
|---|---|
| Target | 50.6 |
| Herbie | 0.3 |
Initial program 61.8
Taylor expanded around 0 0.3
Final simplification0.3
herbie shell --seed 2020150
(FPCore (t)
:name "fma_test1"
:precision binary64
:pre (<= 0.9 t 1.1)
:herbie-target
(fma (+ 1.0 (* t 2e-16)) (+ 1.0 (* t 2e-16)) (- -1.0 (* 2.0 (* t 2e-16))))
(+ (* (+ 1.0 (* t 2e-16)) (+ 1.0 (* t 2e-16))) (- -1.0 (* 2.0 (* t 2e-16)))))