\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{\sqrt{3.999999999999999676487027278085939408227 \cdot 10^{-32}}} \cdot \left(t \cdot \left(t \cdot {\left(\sqrt{\sqrt{3.999999999999999676487027278085939408227 \cdot 10^{-32}}}\right)}^{3}\right)\right)double f(double t) {
double r55542 = 1.0;
double r55543 = t;
double r55544 = 2e-16;
double r55545 = r55543 * r55544;
double r55546 = r55542 + r55545;
double r55547 = r55546 * r55546;
double r55548 = -1.0;
double r55549 = 2.0;
double r55550 = r55549 * r55545;
double r55551 = r55548 - r55550;
double r55552 = r55547 + r55551;
return r55552;
}
double f(double t) {
double r55553 = 3.9999999999999997e-32;
double r55554 = sqrt(r55553);
double r55555 = sqrt(r55554);
double r55556 = t;
double r55557 = 3.0;
double r55558 = pow(r55555, r55557);
double r55559 = r55556 * r55558;
double r55560 = r55556 * r55559;
double r55561 = r55555 * r55560;
return r55561;
}




Bits error versus t
Results
| Original | 61.8 |
|---|---|
| Target | 50.6 |
| Herbie | 0.4 |
Initial program 61.8
Taylor expanded around 0 0.4
rmApplied add-sqr-sqrt0.4
Applied associate-*l*0.4
rmApplied add-sqr-sqrt0.4
Applied sqrt-prod0.4
Applied associate-*l*0.4
Simplified0.4
Final simplification0.4
herbie shell --seed 2019347
(FPCore (t)
:name "fma_test1"
:precision binary64
:pre (<= 0.9 t 1.1)
: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)))))