\frac{\frac{1}{x}}{y \cdot \left(1 + z \cdot z\right)}\frac{1}{y \cdot \left(1 \cdot x\right) + \left(y \cdot z\right) \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 (* 1.0 x)) (* (* y z) (* x z)))))
double code(double x, double y, double z) {
return ((1.0 / x) / ((double) (y * ((double) (1.0 + ((double) (z * z)))))));
}
double code(double x, double y, double z) {
return (1.0 / ((double) (((double) (y * ((double) (1.0 * x)))) + ((double) (((double) (y * z)) * ((double) (x * z)))))));
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.2 |
|---|---|
| Target | 5.5 |
| Herbie | 1.7 |
Initial program Error: 6.2 bits
SimplifiedError: 6.4 bits
rmApplied distribute-lft-inError: 6.4 bits
Applied distribute-lft-inError: 6.4 bits
SimplifiedError: 6.4 bits
SimplifiedError: 3.6 bits
rmApplied associate-*r*Error: 1.7 bits
Final simplificationError: 1.7 bits
herbie shell --seed 2020204
(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)))))