\left(\left(x \cdot y + z \cdot z\right) + z \cdot z\right) + z \cdot z
\left(3 \cdot z\right) \cdot z + x \cdot y
double f(double x, double y, double z) {
double r496106 = x;
double r496107 = y;
double r496108 = r496106 * r496107;
double r496109 = z;
double r496110 = r496109 * r496109;
double r496111 = r496108 + r496110;
double r496112 = r496111 + r496110;
double r496113 = r496112 + r496110;
return r496113;
}
double f(double x, double y, double z) {
double r496114 = 3.0;
double r496115 = z;
double r496116 = r496114 * r496115;
double r496117 = r496116 * r496115;
double r496118 = x;
double r496119 = y;
double r496120 = r496118 * r496119;
double r496121 = r496117 + r496120;
return r496121;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 0.1 |
|---|---|
| Target | 0.1 |
| Herbie | 0.1 |
Initial program 0.1
Simplified0.1
rmApplied fma-udef0.1
rmApplied associate-*r*0.1
Final simplification0.1
herbie shell --seed 2020062 +o rules:numerics
(FPCore (x y z)
:name "Linear.Quaternion:$c/ from linear-1.19.1.3, A"
:precision binary64
:herbie-target
(+ (* (* 3 z) z) (* y x))
(+ (+ (+ (* x y) (* z z)) (* z z)) (* z z)))