Average Error: 0.2 → 0.2
Time: 1.7s
Precision: 64
\[\left(\left(x - \frac{16}{116}\right) \cdot 3\right) \cdot y\]
↓
\[y \cdot \mathsf{fma}\left(3, x, -0.413793103448275856\right)\]
\left(\left(x - \frac{16}{116}\right) \cdot 3\right) \cdot y
↓
y \cdot \mathsf{fma}\left(3, x, -0.413793103448275856\right)
double code(double x, double y) {
	return (((x - (16.0 / 116.0)) * 3.0) * y);
}
↓
double code(double x, double y) {
	return (y * fma(3.0, x, -0.41379310344827586));
}

Error

Bits error versus x

Bits error versus y

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Target

Original0.2
Target0.2
Herbie0.2
\[y \cdot \left(x \cdot 3 - 0.413793103448275856\right)\]

Derivation

  1. Initial program 0.2

    \[\left(\left(x - \frac{16}{116}\right) \cdot 3\right) \cdot y\]

  2. Taylor expanded around 0 0.2

    \[\leadsto \color{blue}{3 \cdot \left(x \cdot y\right) - 0.413793103448275856 \cdot y}\]

  3. Simplified0.2

    \[\leadsto \color{blue}{y \cdot \left(3 \cdot x - 0.413793103448275856\right)}\]
  4. Using strategy rm

  5. Applied fma-neg0.2

    \[\leadsto y \cdot \color{blue}{\mathsf{fma}\left(3, x, -0.413793103448275856\right)}\]

  6. Final simplification0.2

    \[\leadsto y \cdot \mathsf{fma}\left(3, x, -0.413793103448275856\right)\]

Reproduce

herbie shell --seed 2020105 +o rules:numerics
(FPCore (x y)
  :name "Data.Colour.CIE:cieLAB from colour-2.3.3, A"
  :precision binary64

  :herbie-target
  (* y (- (* x 3) 0.41379310344827586))

  (* (* (- x (/ 16 116)) 3) y))