Average Error: 49.2 → 0
Time: 1.9s
Precision: 64
\[1.9 \le t \le 2.1\]
\[1.7 \cdot 10^{+308} \cdot t - 1.7 \cdot 10^{+308}\]
\[\mathsf{fma}\left(1.7 \cdot 10^{+308}, t, -1.7 \cdot 10^{+308}\right)\]
1.7 \cdot 10^{+308} \cdot t - 1.7 \cdot 10^{+308}
\mathsf{fma}\left(1.7 \cdot 10^{+308}, t, -1.7 \cdot 10^{+308}\right)
double f(double t) {
        double r774215 = 1.7e+308;
        double r774216 = t;
        double r774217 = r774215 * r774216;
        double r774218 = r774217 - r774215;
        return r774218;
}

double f(double t) {
        double r774219 = 1.7e+308;
        double r774220 = t;
        double r774221 = -r774219;
        double r774222 = fma(r774219, r774220, r774221);
        return r774222;
}

Error

Bits error versus t

Target

Original49.2
Target0
Herbie0
\[\mathsf{fma}\left(1.7 \cdot 10^{+308}, t, -1.7 \cdot 10^{+308}\right)\]

Derivation

  1. Initial program 49.2

    \[1.7 \cdot 10^{+308} \cdot t - 1.7 \cdot 10^{+308}\]
  2. Using strategy rm
  3. Applied fma-neg0

    \[\leadsto \color{blue}{\mathsf{fma}\left(1.7 \cdot 10^{+308}, t, -1.7 \cdot 10^{+308}\right)}\]
  4. Final simplification0

    \[\leadsto \mathsf{fma}\left(1.7 \cdot 10^{+308}, t, -1.7 \cdot 10^{+308}\right)\]

Reproduce

herbie shell --seed 2019156 +o rules:numerics
(FPCore (t)
  :name "fma_test2"
  :pre (<= 1.9 t 2.1)

  :herbie-target
  (fma 1.7e+308 t (- 1.7e+308))

  (- (* 1.7e+308 t) 1.7e+308))