\frac{x \cdot \left(y + z\right)}{z}\frac{\frac{x}{\frac{\sqrt[3]{z} \cdot \sqrt[3]{z}}{\sqrt[3]{y} \cdot \sqrt[3]{y}}}}{\frac{\sqrt[3]{z}}{\sqrt[3]{y}}} + xdouble f(double x, double y, double z) {
double r248734 = x;
double r248735 = y;
double r248736 = z;
double r248737 = r248735 + r248736;
double r248738 = r248734 * r248737;
double r248739 = r248738 / r248736;
return r248739;
}
double f(double x, double y, double z) {
double r248740 = x;
double r248741 = z;
double r248742 = cbrt(r248741);
double r248743 = r248742 * r248742;
double r248744 = y;
double r248745 = cbrt(r248744);
double r248746 = r248745 * r248745;
double r248747 = r248743 / r248746;
double r248748 = r248740 / r248747;
double r248749 = r248742 / r248745;
double r248750 = r248748 / r248749;
double r248751 = r248750 + r248740;
return r248751;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 12.3 |
|---|---|
| Target | 3.3 |
| Herbie | 0.9 |
Initial program 12.3
Simplified3.5
Taylor expanded around 0 4.9
Simplified4.7
rmApplied fma-udef4.7
Simplified4.9
rmApplied associate-/l*3.3
rmApplied add-cube-cbrt3.6
Applied add-cube-cbrt3.7
Applied times-frac3.7
Applied associate-/r*0.9
Final simplification0.9
herbie shell --seed 2019347 +o rules:numerics
(FPCore (x y z)
:name "Numeric.SpecFunctions:choose from math-functions-0.1.5.2"
:precision binary64
:herbie-target
(/ x (/ z (+ y z)))
(/ (* x (+ y z)) z))