\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 r62846 = 1.0;
double r62847 = t;
double r62848 = 2e-16;
double r62849 = r62847 * r62848;
double r62850 = r62846 + r62849;
double r62851 = r62850 * r62850;
double r62852 = -1.0;
double r62853 = 2.0;
double r62854 = r62853 * r62849;
double r62855 = r62852 - r62854;
double r62856 = r62851 + r62855;
return r62856;
}
double f(double t) {
double r62857 = t;
double r62858 = 3.9999999999999997e-32;
double r62859 = r62857 * r62858;
double r62860 = r62857 * r62859;
return r62860;
}




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