\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)\sqrt{1} \cdot \left(3.9999999999999997 \cdot 10^{-32} \cdot {t}^{2}\right)double f(double t) {
double r60339 = 1.0;
double r60340 = t;
double r60341 = 2e-16;
double r60342 = r60340 * r60341;
double r60343 = r60339 + r60342;
double r60344 = r60343 * r60343;
double r60345 = -1.0;
double r60346 = 2.0;
double r60347 = r60346 * r60342;
double r60348 = r60345 - r60347;
double r60349 = r60344 + r60348;
return r60349;
}
double f(double t) {
double r60350 = 1.0;
double r60351 = sqrt(r60350);
double r60352 = 3.9999999999999997e-32;
double r60353 = t;
double r60354 = 2.0;
double r60355 = pow(r60353, r60354);
double r60356 = r60352 * r60355;
double r60357 = r60351 * r60356;
return r60357;
}




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 add-sqr-sqrt0.3
Applied associate-*l*0.4
rmApplied *-un-lft-identity0.4
Applied sqrt-prod0.4
Applied associate-*l*0.4
Simplified0.3
Final simplification0.3
herbie shell --seed 2020059
(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)))))