\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 r3007876 = 1.0;
double r3007877 = t;
double r3007878 = 2e-16;
double r3007879 = r3007877 * r3007878;
double r3007880 = r3007876 + r3007879;
double r3007881 = r3007880 * r3007880;
double r3007882 = -1.0;
double r3007883 = 2.0;
double r3007884 = r3007883 * r3007879;
double r3007885 = r3007882 - r3007884;
double r3007886 = r3007881 + r3007885;
return r3007886;
}
double f(double t) {
double r3007887 = 3.9999999999999997e-32;
double r3007888 = t;
double r3007889 = r3007887 * r3007888;
double r3007890 = r3007889 * r3007888;
return r3007890;
}




Bits error versus t
Results
| Original | 61.8 |
|---|---|
| Target | 50.6 |
| Herbie | 0.3 |
Initial program 61.8
Simplified50.6
Taylor expanded around 0 0.3
Simplified0.3
Taylor expanded around 0 0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2019165
(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)))))