Average Error: 64.0 → 0
Time: 554.0ms
Precision: 64
\[1.8999999999999999 \le t \le 2.10000000000000009\]
\[1.6999999999999999 \cdot 10^{308} \cdot t - 1.6999999999999999 \cdot 10^{308}\]
↓
\[\mathsf{fma}\left(1.6999999999999999 \cdot 10^{308}, t, -1.6999999999999999 \cdot 10^{308}\right)\]
1.6999999999999999 \cdot 10^{308} \cdot t - 1.6999999999999999 \cdot 10^{308}
↓
\mathsf{fma}\left(1.6999999999999999 \cdot 10^{308}, t, -1.6999999999999999 \cdot 10^{308}\right)
double f(double t) {
        double r64781 = 1.7e+308;
        double r64782 = t;
        double r64783 = r64781 * r64782;
        double r64784 = r64783 - r64781;
        return r64784;
}
↓
double f(double t) {
        double r64785 = 1.7e+308;
        double r64786 = t;
        double r64787 = -r64785;
        double r64788 = fma(r64785, r64786, r64787);
        return r64788;
}

Error

Bits error versus t

Target

Original64.0
Target0
Herbie0
\[\mathsf{fma}\left(1.6999999999999999 \cdot 10^{308}, t, -1.6999999999999999 \cdot 10^{308}\right)\]

Derivation

  1. Initial program 64.0

    \[1.6999999999999999 \cdot 10^{308} \cdot t - 1.6999999999999999 \cdot 10^{308}\]
  2. Using strategy rm

  3. Applied fma-neg0

    \[\leadsto \color{blue}{\mathsf{fma}\left(1.6999999999999999 \cdot 10^{308}, t, -1.6999999999999999 \cdot 10^{308}\right)}\]

  4. Final simplification0

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

Reproduce

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

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

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