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




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 17.1 |
|---|---|
| Target | 0.0 |
| Herbie | 0.0 |
Initial program 17.1
Simplified0.0
rmApplied sub-neg0.0
Applied distribute-lft-in0.0
rmApplied fma-def0.0
Final simplification0.0
herbie shell --seed 2020057 +o rules:numerics
(FPCore (x y z)
:name "Linear.Quaternion:$c/ from linear-1.19.1.3, B"
:precision binary64
:herbie-target
(* (- x z) y)
(+ (- (- (* x y) (* y z)) (* y y)) (* y y)))