\frac{x + y}{1 - \frac{y}{z}}\begin{array}{l}
\mathbf{if}\;\frac{x + y}{1 - \frac{y}{z}} \le -4.9425289194373917 \cdot 10^{-297}:\\
\;\;\;\;\frac{x + y}{1 - \frac{y}{z}}\\
\mathbf{elif}\;\frac{x + y}{1 - \frac{y}{z}} \le -0.0:\\
\;\;\;\;\frac{1}{\sqrt{1} + \frac{\sqrt{y}}{\sqrt{z}}} \cdot \frac{x + y}{\sqrt{1} - \frac{\sqrt{y}}{\sqrt{z}}}\\
\mathbf{else}:\\
\;\;\;\;\left(x + y\right) \cdot \frac{1}{1 - \frac{y}{z}}\\
\end{array}double code(double x, double y, double z) {
return ((double) (((double) (x + y)) / ((double) (1.0 - ((double) (y / z))))));
}
double code(double x, double y, double z) {
double VAR;
if ((((double) (((double) (x + y)) / ((double) (1.0 - ((double) (y / z)))))) <= -4.942528919437392e-297)) {
VAR = ((double) (((double) (x + y)) / ((double) (1.0 - ((double) (y / z))))));
} else {
double VAR_1;
if ((((double) (((double) (x + y)) / ((double) (1.0 - ((double) (y / z)))))) <= -0.0)) {
VAR_1 = ((double) (((double) (1.0 / ((double) (((double) sqrt(1.0)) + ((double) (((double) sqrt(y)) / ((double) sqrt(z)))))))) * ((double) (((double) (x + y)) / ((double) (((double) sqrt(1.0)) - ((double) (((double) sqrt(y)) / ((double) sqrt(z))))))))));
} else {
VAR_1 = ((double) (((double) (x + y)) * ((double) (1.0 / ((double) (1.0 - ((double) (y / z))))))));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 7.5 |
|---|---|
| Target | 4.1 |
| Herbie | 6.2 |
if (/ (+ x y) (- 1.0 (/ y z))) < -4.942528919437392e-297Initial program 0.1
if -4.942528919437392e-297 < (/ (+ x y) (- 1.0 (/ y z))) < -0.0Initial program 59.2
rmApplied add-sqr-sqrt61.2
Applied add-sqr-sqrt62.8
Applied times-frac62.8
Applied add-sqr-sqrt62.8
Applied difference-of-squares62.8
Applied *-un-lft-identity62.8
Applied times-frac48.9
if -0.0 < (/ (+ x y) (- 1.0 (/ y z))) Initial program 0.1
rmApplied div-inv0.1
Final simplification6.2
herbie shell --seed 2020128
(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))))