\left(1 + t \cdot 1.999999999999999958195573448069207123682 \cdot 10^{-16}\right) \cdot \left(1 + t \cdot 1.999999999999999958195573448069207123682 \cdot 10^{-16}\right) + \left(-1 - 2 \cdot \left(t \cdot 1.999999999999999958195573448069207123682 \cdot 10^{-16}\right)\right)t \cdot \left(t \cdot 3.999999999999999676487027278085939408227 \cdot 10^{-32}\right)double f(double t) {
double r5426231 = 1.0;
double r5426232 = t;
double r5426233 = 2e-16;
double r5426234 = r5426232 * r5426233;
double r5426235 = r5426231 + r5426234;
double r5426236 = r5426235 * r5426235;
double r5426237 = -1.0;
double r5426238 = 2.0;
double r5426239 = r5426238 * r5426234;
double r5426240 = r5426237 - r5426239;
double r5426241 = r5426236 + r5426240;
return r5426241;
}
double f(double t) {
double r5426242 = t;
double r5426243 = 3.9999999999999997e-32;
double r5426244 = r5426242 * r5426243;
double r5426245 = r5426242 * r5426244;
return r5426245;
}




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.4
Simplified0.4
rmApplied associate-*l*0.3
Final simplification0.3
herbie shell --seed 2019174 +o rules:numerics
(FPCore (t)
:name "fma_test1"
:pre (<= 0.9 t 1.1)
:herbie-target
(fma (+ 1.0 (* t 2e-16)) (+ 1.0 (* t 2e-16)) (- -1.0 (* 2.0 (* t 2e-16))))
(+ (* (+ 1.0 (* t 2e-16)) (+ 1.0 (* t 2e-16))) (- -1.0 (* 2.0 (* t 2e-16)))))