\frac{x + y}{1 - \frac{y}{z}}\begin{array}{l}
\mathbf{if}\;\frac{x + y}{1 - \frac{y}{z}} \leq -1.4483950108988509 \cdot 10^{-297} \lor \neg \left(\frac{x + y}{1 - \frac{y}{z}} \leq 0\right):\\
\;\;\;\;\frac{x + y}{1 - \frac{y}{z}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{1 + \frac{\sqrt{y}}{\sqrt{z}}}}{\frac{1 - \frac{\sqrt{y}}{\sqrt{z}}}{x + y}}\\
\end{array}(FPCore (x y z) :precision binary64 (/ (+ x y) (- 1.0 (/ y z))))
(FPCore (x y z)
:precision binary64
(if (or (<= (/ (+ x y) (- 1.0 (/ y z))) -1.4483950108988509e-297)
(not (<= (/ (+ x y) (- 1.0 (/ y z))) 0.0)))
(/ (+ x y) (- 1.0 (/ y z)))
(/
(/ 1.0 (+ 1.0 (/ (sqrt y) (sqrt z))))
(/ (- 1.0 (/ (sqrt y) (sqrt z))) (+ x y)))))double code(double x, double y, double z) {
return (x + y) / (1.0 - (y / z));
}
double code(double x, double y, double z) {
double tmp;
if ((((x + y) / (1.0 - (y / z))) <= -1.4483950108988509e-297) || !(((x + y) / (1.0 - (y / z))) <= 0.0)) {
tmp = (x + y) / (1.0 - (y / z));
} else {
tmp = (1.0 / (1.0 + (sqrt(y) / sqrt(z)))) / ((1.0 - (sqrt(y) / sqrt(z))) / (x + y));
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 7.8 |
|---|---|
| Target | 4.2 |
| Herbie | 6.3 |
if (/.f64 (+.f64 x y) (-.f64 1 (/.f64 y z))) < -1.44839501089885086e-297 or 0.0 < (/.f64 (+.f64 x y) (-.f64 1 (/.f64 y z))) Initial program 0.1
if -1.44839501089885086e-297 < (/.f64 (+.f64 x y) (-.f64 1 (/.f64 y z))) < 0.0Initial program 59.4
rmApplied clear-num_binary64_984959.4
rmApplied *-un-lft-identity_binary64_985059.4
Applied add-sqr-sqrt_binary64_987161.4
Applied add-sqr-sqrt_binary64_987162.7
Applied times-frac_binary64_985662.7
Applied *-un-lft-identity_binary64_985062.7
Applied difference-of-squares_binary64_981962.7
Applied times-frac_binary64_985647.7
Applied associate-/r*_binary64_979647.6
Simplified47.6
Final simplification6.3
herbie shell --seed 2020281
(FPCore (x y z)
:name "Graphics.Rendering.Chart.Backend.Diagrams:calcFontMetrics from Chart-diagrams-1.5.1, A"
:precision binary64
:herbie-target
(if (< y -3.7429310762689856e+171) (* (/ (+ y x) (- y)) z) (if (< y 3.5534662456086734e+168) (/ (+ x y) (- 1.0 (/ y z))) (* (/ (+ y x) (- y)) z)))
(/ (+ x y) (- 1.0 (/ y z))))