\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(t \cdot \left(t \cdot \sqrt{3.999999999999999676487027278085939408227 \cdot 10^{-32}}\right)\right)double f(double t) {
double r4065704 = 1.0;
double r4065705 = t;
double r4065706 = 2e-16;
double r4065707 = r4065705 * r4065706;
double r4065708 = r4065704 + r4065707;
double r4065709 = r4065708 * r4065708;
double r4065710 = -1.0;
double r4065711 = 2.0;
double r4065712 = r4065711 * r4065707;
double r4065713 = r4065710 - r4065712;
double r4065714 = r4065709 + r4065713;
return r4065714;
}
double f(double t) {
double r4065715 = 3.9999999999999997e-32;
double r4065716 = sqrt(r4065715);
double r4065717 = t;
double r4065718 = r4065717 * r4065716;
double r4065719 = r4065717 * r4065718;
double r4065720 = r4065716 * r4065719;
return r4065720;
}




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