x \cdot \left(1 - \left(1 - y\right) \cdot z\right)
\mathsf{fma}\left(y - 1, z, 1\right) \cdot xdouble code(double x, double y, double z) {
return (x * (1.0 - ((1.0 - y) * z)));
}
double code(double x, double y, double z) {
return (fma((y - 1.0), z, 1.0) * x);
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 3.4 |
|---|---|
| Target | 0.2 |
| Herbie | 3.4 |
Initial program 3.4
Simplified3.4
Final simplification3.4
herbie shell --seed 2020057 +o rules:numerics
(FPCore (x y z)
:name "Data.Colour.RGBSpace.HSV:hsv from colour-2.3.3, J"
:precision binary64
:herbie-target
(if (< (* x (- 1 (* (- 1 y) z))) -1.618195973607049e+50) (+ x (* (- 1 y) (* (- z) x))) (if (< (* x (- 1 (* (- 1 y) z))) 3.892237649663903e+134) (- (* (* x y) z) (- (* x z) x)) (+ x (* (- 1 y) (* (- z) x)))))
(* x (- 1 (* (- 1 y) z))))