Average Error: 49.1 → 0
Time: 4.1s
Precision: 64
\[1.9 \le t \le 2.1\]
\[1.7 \cdot 10^{+308} \cdot t - 1.7 \cdot 10^{+308}\]
↓
\[\mathsf{fma}\left(\left( 1.7 \cdot 10^{+308} \right), t, \left(-1.7 \cdot 10^{+308}\right)\right)\]
1.7 \cdot 10^{+308} \cdot t - 1.7 \cdot 10^{+308}
↓
\mathsf{fma}\left(\left( 1.7 \cdot 10^{+308} \right), t, \left(-1.7 \cdot 10^{+308}\right)\right)
double f(double t) {
        double r8639614 = 1.7e+308;
        double r8639615 = t;
        double r8639616 = r8639614 * r8639615;
        double r8639617 = r8639616 - r8639614;
        return r8639617;
}
↓
double f(double t) {
        double r8639618 = 1.7e+308;
        double r8639619 = t;
        double r8639620 = -r8639618;
        double r8639621 = fma(r8639618, r8639619, r8639620);
        return r8639621;
}

Error

Bits error versus t

Target

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

Derivation

  1. Initial program 49.1

    \[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(\left( 1.7 \cdot 10^{+308} \right), t, \left(-1.7 \cdot 10^{+308}\right)\right)}\]

  4. Final simplification0

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

Reproduce

herbie shell --seed 2019128 +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))