\frac{\frac{1}{x}}{y \cdot \left(1 + z \cdot z\right)}\frac{\frac{1}{y}}{x + z \cdot \left(x \cdot z\right)}(FPCore (x y z) :precision binary64 (/ (/ 1.0 x) (* y (+ 1.0 (* z z)))))
(FPCore (x y z) :precision binary64 (/ (/ 1.0 y) (+ x (* z (* x z)))))
double code(double x, double y, double z) {
return (1.0 / x) / (y * (1.0 + (z * z)));
}
double code(double x, double y, double z) {
return (1.0 / y) / (x + (z * (x * z)));
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.5 |
|---|---|
| Target | 5.1 |
| Herbie | 3.8 |
Initial program 6.5
rmApplied *-un-lft-identity_binary64_109906.5
Applied add-sqr-sqrt_binary64_110126.5
Applied times-frac_binary64_109966.5
Applied times-frac_binary64_109966.4
Simplified6.4
Simplified6.6
rmApplied associate-*r*_binary64_109303.8
rmApplied un-div-inv_binary64_109883.8
Final simplification3.8
herbie shell --seed 2021015
(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)))))