\frac{x + y}{1 - \frac{y}{z}}\begin{array}{l}
\mathbf{if}\;\frac{x + y}{1 - \frac{y}{z}} \le -1.8853869057044877 \cdot 10^{-283} \lor \neg \left(\frac{x + y}{1 - \frac{y}{z}} \le 0.0\right):\\
\;\;\;\;\frac{x + y}{1 - \frac{y}{z}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\sqrt{1} + \frac{\sqrt{y}}{\sqrt{z}}}{\sqrt{x + y}} \cdot \frac{\sqrt{1} - \frac{\sqrt{y}}{\sqrt{z}}}{\sqrt{x + y}}}\\
\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)))))) <= -1.8853869057044877e-283) || !(((double) (((double) (x + y)) / ((double) (1.0 - ((double) (y / z)))))) <= 0.0))) {
VAR = ((double) (((double) (x + y)) / ((double) (1.0 - ((double) (y / z))))));
} else {
VAR = ((double) (1.0 / ((double) (((double) (((double) (((double) sqrt(1.0)) + ((double) (((double) sqrt(y)) / ((double) sqrt(z)))))) / ((double) sqrt(((double) (x + y)))))) * ((double) (((double) (((double) sqrt(1.0)) - ((double) (((double) sqrt(y)) / ((double) sqrt(z)))))) / ((double) sqrt(((double) (x + y))))))))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 7.5 |
|---|---|
| Target | 4.0 |
| Herbie | 6.4 |
if (/ (+ x y) (- 1.0 (/ y z))) < -1.8853869057044877e-283 or 0.0 < (/ (+ x y) (- 1.0 (/ y z))) Initial program 0.1
if -1.8853869057044877e-283 < (/ (+ x y) (- 1.0 (/ y z))) < 0.0Initial program 58.5
rmApplied clear-num58.5
rmApplied add-sqr-sqrt61.2
Applied add-sqr-sqrt62.1
Applied add-sqr-sqrt62.2
Applied times-frac62.2
Applied add-sqr-sqrt62.2
Applied difference-of-squares62.2
Applied times-frac49.7
Final simplification6.4
herbie shell --seed 2020174
(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) (neg y)) z) (if (< y 3.5534662456086734e+168) (/ (+ x y) (- 1.0 (/ y z))) (* (/ (+ y x) (neg y)) z)))
(/ (+ x y) (- 1.0 (/ y z))))