\left(\left(x \cdot y + z \cdot z\right) + z \cdot z\right) + z \cdot z
\left(3 \cdot \left|z\right|\right) \cdot \sqrt{{z}^{2}} + x \cdot ydouble f(double x, double y, double z) {
double r520856 = x;
double r520857 = y;
double r520858 = r520856 * r520857;
double r520859 = z;
double r520860 = r520859 * r520859;
double r520861 = r520858 + r520860;
double r520862 = r520861 + r520860;
double r520863 = r520862 + r520860;
return r520863;
}
double f(double x, double y, double z) {
double r520864 = 3.0;
double r520865 = z;
double r520866 = fabs(r520865);
double r520867 = r520864 * r520866;
double r520868 = 2.0;
double r520869 = pow(r520865, r520868);
double r520870 = sqrt(r520869);
double r520871 = r520867 * r520870;
double r520872 = x;
double r520873 = y;
double r520874 = r520872 * r520873;
double r520875 = r520871 + r520874;
return r520875;
}




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
Taylor expanded around 0 0.1
Simplified0.1
rmApplied fma-udef0.1
rmApplied add-sqr-sqrt0.1
Applied associate-*r*0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2020064 +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)))