\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(t \cdot \left(t \cdot \sqrt{3.9999999999999997 \cdot 10^{-32}}\right)\right) \cdot \sqrt{3.9999999999999997 \cdot 10^{-32}}double f(double t) {
double r70725 = 1.0;
double r70726 = t;
double r70727 = 2e-16;
double r70728 = r70726 * r70727;
double r70729 = r70725 + r70728;
double r70730 = r70729 * r70729;
double r70731 = -1.0;
double r70732 = 2.0;
double r70733 = r70732 * r70728;
double r70734 = r70731 - r70733;
double r70735 = r70730 + r70734;
return r70735;
}
double f(double t) {
double r70736 = t;
double r70737 = 3.9999999999999997e-32;
double r70738 = sqrt(r70737);
double r70739 = r70736 * r70738;
double r70740 = r70736 * r70739;
double r70741 = r70740 * r70738;
return r70741;
}




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 add-sqr-sqrt0.4
Applied associate-*l*0.4
rmApplied sqr-pow0.4
Applied associate-*r*0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2020065
(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)))))