\frac{x \cdot y}{\left(z \cdot z\right) \cdot \left(z + 1\right)}\left(\sqrt[3]{\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{z} \cdot \left(\frac{\sqrt[3]{x}}{z} \cdot \frac{y}{z + 1}\right)} \cdot \sqrt[3]{\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{z} \cdot \left(\frac{\sqrt[3]{x}}{z} \cdot \frac{y}{z + 1}\right)}\right) \cdot \sqrt[3]{\frac{\sqrt[3]{x} \cdot \sqrt[3]{x}}{z} \cdot \left(\frac{\sqrt[3]{x}}{z} \cdot \frac{y}{z + 1}\right)}double f(double x, double y, double z) {
double r380157 = x;
double r380158 = y;
double r380159 = r380157 * r380158;
double r380160 = z;
double r380161 = r380160 * r380160;
double r380162 = 1.0;
double r380163 = r380160 + r380162;
double r380164 = r380161 * r380163;
double r380165 = r380159 / r380164;
return r380165;
}
double f(double x, double y, double z) {
double r380166 = x;
double r380167 = cbrt(r380166);
double r380168 = r380167 * r380167;
double r380169 = z;
double r380170 = r380168 / r380169;
double r380171 = r380167 / r380169;
double r380172 = y;
double r380173 = 1.0;
double r380174 = r380169 + r380173;
double r380175 = r380172 / r380174;
double r380176 = r380171 * r380175;
double r380177 = r380170 * r380176;
double r380178 = cbrt(r380177);
double r380179 = r380178 * r380178;
double r380180 = r380179 * r380178;
return r380180;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 14.6 |
|---|---|
| Target | 4.1 |
| Herbie | 1.2 |
Initial program 14.6
rmApplied times-frac11.0
rmApplied add-cube-cbrt11.3
Applied times-frac6.4
Applied associate-*l*1.1
rmApplied add-cube-cbrt1.2
Final simplification1.2
herbie shell --seed 2020057
(FPCore (x y z)
:name "Statistics.Distribution.Beta:$cvariance from math-functions-0.1.5.2"
:precision binary64
:herbie-target
(if (< z 249.6182814532307) (/ (* y (/ x z)) (+ z (* z z))) (/ (* (/ (/ y z) (+ 1 z)) x) z))
(/ (* x y) (* (* z z) (+ z 1))))