\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)\left({t}^{2} \cdot \sqrt{3.999999999999999676487027278085939408227 \cdot 10^{-32}}\right) \cdot \sqrt{3.999999999999999676487027278085939408227 \cdot 10^{-32}}double f(double t) {
double r95108 = 1.0;
double r95109 = t;
double r95110 = 2e-16;
double r95111 = r95109 * r95110;
double r95112 = r95108 + r95111;
double r95113 = r95112 * r95112;
double r95114 = -1.0;
double r95115 = 2.0;
double r95116 = r95115 * r95111;
double r95117 = r95114 - r95116;
double r95118 = r95113 + r95117;
return r95118;
}
double f(double t) {
double r95119 = t;
double r95120 = 2.0;
double r95121 = pow(r95119, r95120);
double r95122 = 3.9999999999999997e-32;
double r95123 = sqrt(r95122);
double r95124 = r95121 * r95123;
double r95125 = r95124 * r95123;
return r95125;
}




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.3
rmApplied add-sqr-sqrt0.3
Applied associate-*l*0.3
rmApplied *-un-lft-identity0.3
Applied associate-*l*0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2019303 +o rules:numerics
(FPCore (t)
:name "fma_test1"
:precision binary64
:pre (<= 0.900000000000000022 t 1.1000000000000001)
: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)))))