\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 r471899 = x;
double r471900 = y;
double r471901 = r471899 * r471900;
double r471902 = z;
double r471903 = r471902 * r471902;
double r471904 = r471901 + r471903;
double r471905 = r471904 + r471903;
double r471906 = r471905 + r471903;
return r471906;
}
double f(double x, double y, double z) {
double r471907 = 3.0;
double r471908 = z;
double r471909 = r471907 * r471908;
double r471910 = r471909 * r471908;
double r471911 = x;
double r471912 = y;
double r471913 = r471911 * r471912;
double r471914 = r471910 + r471913;
return r471914;
}




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 2019353 +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)))