\frac{\frac{1}{x}}{y \cdot \left(1 + z \cdot z\right)}\left(\left(\sqrt[3]{\frac{1}{y} \cdot \frac{\sqrt[3]{\frac{1}{x}} \cdot \sqrt[3]{\frac{1}{x}}}{\sqrt{1 + z \cdot z}}} \cdot \sqrt[3]{\frac{1}{y} \cdot \frac{\sqrt[3]{\frac{1}{x}} \cdot \sqrt[3]{\frac{1}{x}}}{\sqrt{1 + z \cdot z}}}\right) \cdot \sqrt[3]{\frac{1}{y} \cdot \frac{\sqrt[3]{\frac{1}{x}} \cdot \sqrt[3]{\frac{1}{x}}}{\sqrt{1 + z \cdot z}}}\right) \cdot \frac{\sqrt[3]{\frac{1}{x}}}{\sqrt{1 + z \cdot z}}double f(double x, double y, double z) {
double r349660 = 1.0;
double r349661 = x;
double r349662 = r349660 / r349661;
double r349663 = y;
double r349664 = z;
double r349665 = r349664 * r349664;
double r349666 = r349660 + r349665;
double r349667 = r349663 * r349666;
double r349668 = r349662 / r349667;
return r349668;
}
double f(double x, double y, double z) {
double r349669 = 1.0;
double r349670 = y;
double r349671 = r349669 / r349670;
double r349672 = 1.0;
double r349673 = x;
double r349674 = r349672 / r349673;
double r349675 = cbrt(r349674);
double r349676 = r349675 * r349675;
double r349677 = z;
double r349678 = r349677 * r349677;
double r349679 = r349669 + r349678;
double r349680 = sqrt(r349679);
double r349681 = r349676 / r349680;
double r349682 = r349671 * r349681;
double r349683 = cbrt(r349682);
double r349684 = r349683 * r349683;
double r349685 = r349684 * r349683;
double r349686 = r349675 / r349680;
double r349687 = r349685 * r349686;
return r349687;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.5 |
|---|---|
| Target | 5.7 |
| Herbie | 5.9 |
Initial program 6.5
rmApplied div-inv6.5
Applied times-frac6.5
rmApplied add-sqr-sqrt6.5
Applied add-cube-cbrt7.1
Applied times-frac7.1
Applied associate-*r*5.7
rmApplied add-cube-cbrt5.9
Final simplification5.9
herbie shell --seed 2020025
(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)))))