\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(t \cdot \sqrt{3.999999999999999676487027278085939408227 \cdot 10^{-32}}\right) \cdot t\right)double f(double t) {
double r93923 = 1.0;
double r93924 = t;
double r93925 = 2e-16;
double r93926 = r93924 * r93925;
double r93927 = r93923 + r93926;
double r93928 = r93927 * r93927;
double r93929 = -1.0;
double r93930 = 2.0;
double r93931 = r93930 * r93926;
double r93932 = r93929 - r93931;
double r93933 = r93928 + r93932;
return r93933;
}
double f(double t) {
double r93934 = 3.9999999999999997e-32;
double r93935 = sqrt(r93934);
double r93936 = t;
double r93937 = r93936 * r93935;
double r93938 = r93937 * r93936;
double r93939 = r93935 * r93938;
return r93939;
}




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 unpow20.4
Applied associate-*r*0.3
Simplified0.3
Final simplification0.3
herbie shell --seed 2019209 +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)))))