\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 r500903 = x;
double r500904 = y;
double r500905 = r500903 * r500904;
double r500906 = z;
double r500907 = r500906 * r500906;
double r500908 = r500905 + r500907;
double r500909 = r500908 + r500907;
double r500910 = r500909 + r500907;
return r500910;
}
double f(double x, double y, double z) {
double r500911 = 3.0;
double r500912 = z;
double r500913 = r500911 * r500912;
double r500914 = r500913 * r500912;
double r500915 = x;
double r500916 = y;
double r500917 = r500915 * r500916;
double r500918 = r500914 + r500917;
return r500918;
}




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