\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 r43733 = 1.0;
double r43734 = t;
double r43735 = 2e-16;
double r43736 = r43734 * r43735;
double r43737 = r43733 + r43736;
double r43738 = r43737 * r43737;
double r43739 = -1.0;
double r43740 = 2.0;
double r43741 = r43740 * r43736;
double r43742 = r43739 - r43741;
double r43743 = r43738 + r43742;
return r43743;
}
double f(double t) {
double r43744 = t;
double r43745 = 3.9999999999999997e-32;
double r43746 = r43744 * r43745;
double r43747 = r43744 * r43746;
return r43747;
}




Bits error versus t
Results
| Original | 61.8 |
|---|---|
| Target | 50.6 |
| Herbie | 0.3 |
Initial program 61.8
Taylor expanded around 0 0.3
rmApplied sqr-pow0.3
Applied associate-*r*0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2020060
(FPCore (t)
:name "fma_test1"
:precision binary64
: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)))))