\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)\sqrt{3.999999999999999676487027278085939408227 \cdot 10^{-32}} \cdot \left(\left(\left|t\right| \cdot \sqrt{3.999999999999999676487027278085939408227 \cdot 10^{-32}}\right) \cdot \sqrt{{t}^{2}}\right)double f(double t) {
double r44126 = 1.0;
double r44127 = t;
double r44128 = 2e-16;
double r44129 = r44127 * r44128;
double r44130 = r44126 + r44129;
double r44131 = r44130 * r44130;
double r44132 = -1.0;
double r44133 = 2.0;
double r44134 = r44133 * r44129;
double r44135 = r44132 - r44134;
double r44136 = r44131 + r44135;
return r44136;
}
double f(double t) {
double r44137 = 3.9999999999999997e-32;
double r44138 = sqrt(r44137);
double r44139 = t;
double r44140 = fabs(r44139);
double r44141 = r44140 * r44138;
double r44142 = 2.0;
double r44143 = pow(r44139, r44142);
double r44144 = sqrt(r44143);
double r44145 = r44141 * r44144;
double r44146 = r44138 * r44145;
return r44146;
}




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.4
rmApplied add-sqr-sqrt0.4
Applied associate-*r*0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2019235 +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)))))