\frac{\frac{1}{x}}{y \cdot \left(1 + z \cdot z\right)}\left(\frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{y} \cdot \frac{\frac{\sqrt[3]{1}}{x}}{\sqrt{1 + z \cdot z}}\right) \cdot \frac{1}{\sqrt{1 + z \cdot z}}double f(double x, double y, double z) {
double r387509 = 1.0;
double r387510 = x;
double r387511 = r387509 / r387510;
double r387512 = y;
double r387513 = z;
double r387514 = r387513 * r387513;
double r387515 = r387509 + r387514;
double r387516 = r387512 * r387515;
double r387517 = r387511 / r387516;
return r387517;
}
double f(double x, double y, double z) {
double r387518 = 1.0;
double r387519 = cbrt(r387518);
double r387520 = r387519 * r387519;
double r387521 = y;
double r387522 = r387520 / r387521;
double r387523 = x;
double r387524 = r387519 / r387523;
double r387525 = z;
double r387526 = r387525 * r387525;
double r387527 = r387518 + r387526;
double r387528 = sqrt(r387527);
double r387529 = r387524 / r387528;
double r387530 = r387522 * r387529;
double r387531 = 1.0;
double r387532 = r387531 / r387528;
double r387533 = r387530 * r387532;
return r387533;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.1 |
|---|---|
| Target | 5.5 |
| Herbie | 6.0 |
Initial program 6.1
rmApplied *-un-lft-identity6.1
Applied add-cube-cbrt6.1
Applied times-frac6.1
Applied times-frac6.4
Simplified6.4
rmApplied add-sqr-sqrt6.4
Applied associate-/r*6.4
rmApplied div-inv6.4
Applied associate-*r*6.0
Final simplification6.0
herbie shell --seed 2020047
(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)))))