\frac{\frac{1}{x}}{y \cdot \left(1 + z \cdot z\right)}\begin{array}{l}
\mathbf{if}\;y \cdot \left(1 + z \cdot z\right) \le -1.19195819535175784 \cdot 10^{32} \lor \neg \left(y \cdot \left(1 + z \cdot z\right) \le 7.24453765483583405 \cdot 10^{88}\right):\\
\;\;\;\;\frac{\frac{1}{\sqrt{1 + z \cdot z}}}{y \cdot \left(\sqrt{1 + z \cdot z} \cdot x\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{x}}{\left(y \cdot \sqrt{1 + z \cdot z}\right) \cdot \sqrt{1 + z \cdot z}}\\
\end{array}double code(double x, double y, double z) {
return ((double) (((double) (1.0 / x)) / ((double) (y * ((double) (1.0 + ((double) (z * z))))))));
}
double code(double x, double y, double z) {
double VAR;
if (((((double) (y * ((double) (1.0 + ((double) (z * z)))))) <= -1.1919581953517578e+32) || !(((double) (y * ((double) (1.0 + ((double) (z * z)))))) <= 7.244537654835834e+88))) {
VAR = ((double) (((double) (1.0 / ((double) sqrt(((double) (1.0 + ((double) (z * z)))))))) / ((double) (y * ((double) (((double) sqrt(((double) (1.0 + ((double) (z * z)))))) * x))))));
} else {
VAR = ((double) (((double) (1.0 / x)) / ((double) (((double) (y * ((double) sqrt(((double) (1.0 + ((double) (z * z)))))))) * ((double) sqrt(((double) (1.0 + ((double) (z * z))))))))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.3 |
|---|---|
| Target | 5.5 |
| Herbie | 5.0 |
if (* y (+ 1.0 (* z z))) < -1.1919581953517578e+32 or 7.244537654835834e+88 < (* y (+ 1.0 (* z z))) Initial program 9.4
rmApplied *-un-lft-identity9.4
Applied *-un-lft-identity9.4
Applied times-frac9.4
Applied times-frac7.5
Simplified7.5
rmApplied associate-*l/7.4
Simplified7.4
rmApplied add-sqr-sqrt7.4
Applied div-inv7.4
Applied times-frac7.4
Applied associate-/l*7.4
Simplified7.4
if -1.1919581953517578e+32 < (* y (+ 1.0 (* z z))) < 7.244537654835834e+88Initial program 0.3
rmApplied add-sqr-sqrt0.3
Applied associate-*r*0.3
Final simplification5.0
herbie shell --seed 2020126
(FPCore (x y z)
:name "Statistics.Distribution.CauchyLorentz:$cdensity from math-functions-0.1.5.2"
:precision binary64
:herbie-target
(if (< (* y (+ 1.0 (* z z))) (- INFINITY)) (/ (/ 1.0 y) (* (+ 1.0 (* z z)) x)) (if (< (* y (+ 1.0 (* z z))) 8.680743250567252e+305) (/ (/ 1.0 x) (* (+ 1.0 (* z z)) y)) (/ (/ 1.0 y) (* (+ 1.0 (* z z)) x))))
(/ (/ 1.0 x) (* y (+ 1.0 (* z z)))))