\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 r1301463 = 1.0;
double r1301464 = t;
double r1301465 = 2e-16;
double r1301466 = r1301464 * r1301465;
double r1301467 = r1301463 + r1301466;
double r1301468 = r1301467 * r1301467;
double r1301469 = -1.0;
double r1301470 = 2.0;
double r1301471 = r1301470 * r1301466;
double r1301472 = r1301469 - r1301471;
double r1301473 = r1301468 + r1301472;
return r1301473;
}
double f(double t) {
double r1301474 = t;
double r1301475 = 3.9999999999999997e-32;
double r1301476 = r1301474 * r1301475;
double r1301477 = r1301474 * r1301476;
return r1301477;
}




Bits error versus t
Results
| Original | 61.8 |
|---|---|
| Target | 50.6 |
| Herbie | 0.3 |
Initial program 61.8
Simplified61.8
Taylor expanded around -inf 0.3
Simplified0.3
Taylor expanded around 0 0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2019152 +o rules:numerics
(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)))))