\frac{\frac{1}{x}}{y \cdot \left(1 + z \cdot z\right)}\frac{\frac{1}{\sqrt[3]{x} \cdot \sqrt[3]{x}}}{\left|\sqrt[3]{1 + z \cdot z}\right| \cdot y} \cdot \frac{\frac{\frac{1}{\sqrt[3]{x}}}{\sqrt{1 + z \cdot z}}}{\sqrt{\sqrt[3]{1 + z \cdot z}}}double f(double x, double y, double z) {
double r327305 = 1.0;
double r327306 = x;
double r327307 = r327305 / r327306;
double r327308 = y;
double r327309 = z;
double r327310 = r327309 * r327309;
double r327311 = r327305 + r327310;
double r327312 = r327308 * r327311;
double r327313 = r327307 / r327312;
return r327313;
}
double f(double x, double y, double z) {
double r327314 = 1.0;
double r327315 = x;
double r327316 = cbrt(r327315);
double r327317 = r327316 * r327316;
double r327318 = r327314 / r327317;
double r327319 = 1.0;
double r327320 = z;
double r327321 = r327320 * r327320;
double r327322 = r327319 + r327321;
double r327323 = cbrt(r327322);
double r327324 = fabs(r327323);
double r327325 = y;
double r327326 = r327324 * r327325;
double r327327 = r327318 / r327326;
double r327328 = r327319 / r327316;
double r327329 = sqrt(r327322);
double r327330 = r327328 / r327329;
double r327331 = sqrt(r327323);
double r327332 = r327330 / r327331;
double r327333 = r327327 * r327332;
return r327333;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.9 |
|---|---|
| Target | 6.0 |
| Herbie | 6.4 |
Initial program 6.9
rmApplied add-sqr-sqrt6.9
Applied associate-*r*6.9
rmApplied associate-/r*6.3
rmApplied add-cube-cbrt6.4
Applied sqrt-prod6.4
Applied add-cube-cbrt6.9
Applied *-un-lft-identity6.9
Applied times-frac6.9
Applied times-frac6.7
Applied times-frac6.6
Simplified6.4
Final simplification6.4
herbie shell --seed 2020018
(FPCore (x y z)
:name "Statistics.Distribution.CauchyLorentz:$cdensity from math-functions-0.1.5.2"
:precision binary64
:herbie-target
(if (< (* y (+ 1 (* z z))) #f) (/ (/ 1 y) (* (+ 1 (* z z)) x)) (if (< (* y (+ 1 (* z z))) 8.680743250567252e+305) (/ (/ 1 x) (* (+ 1 (* z z)) y)) (/ (/ 1 y) (* (+ 1 (* z z)) x))))
(/ (/ 1 x) (* y (+ 1 (* z z)))))