Average Error: 0.0 → 0
Time: 4.6s
Precision: 64
\[x \cdot x - 1\]
\[\mathsf{fma}\left(x, x, -1\right)\]
x \cdot x - 1
\mathsf{fma}\left(x, x, -1\right)
double f(double x) {
        double r105061 = x;
        double r105062 = r105061 * r105061;
        double r105063 = 1.0;
        double r105064 = r105062 - r105063;
        return r105064;
}

double f(double x) {
        double r105065 = x;
        double r105066 = 1.0;
        double r105067 = -r105066;
        double r105068 = fma(r105065, r105065, r105067);
        return r105068;
}

Error

Bits error versus x

Derivation

  1. Initial program 0.0

    \[x \cdot x - 1\]
  2. Using strategy rm
  3. Applied fma-neg0

    \[\leadsto \color{blue}{\mathsf{fma}\left(x, x, -1\right)}\]
  4. Final simplification0

    \[\leadsto \mathsf{fma}\left(x, x, -1\right)\]

Reproduce

herbie shell --seed 2019323 +o rules:numerics
(FPCore (x)
  :name "Data.Random.Dice:roll from dice-0.1"
  :precision binary64
  (- (* x x) 1))