\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 r57497 = 1.0;
double r57498 = t;
double r57499 = 2e-16;
double r57500 = r57498 * r57499;
double r57501 = r57497 + r57500;
double r57502 = r57501 * r57501;
double r57503 = -1.0;
double r57504 = 2.0;
double r57505 = r57504 * r57500;
double r57506 = r57503 - r57505;
double r57507 = r57502 + r57506;
return r57507;
}
double f(double t) {
double r57508 = t;
double r57509 = 3.9999999999999997e-32;
double r57510 = r57508 * r57509;
double r57511 = r57508 * r57510;
return r57511;
}




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
rmApplied associate-*l*0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2020018
(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)))))