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




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 7.5 |
|---|---|
| Target | 4.0 |
| Herbie | 6.3 |
if (/ (+ x y) (- 1.0 (/ y z))) < -1.2144667589559522e-305 or 0.0 < (/ (+ x y) (- 1.0 (/ y z))) Initial program 0.1
if -1.2144667589559522e-305 < (/ (+ x y) (- 1.0 (/ y z))) < 0.0Initial program 59.8
rmApplied add-sqr-sqrt61.8
Applied add-sqr-sqrt63.0
Applied times-frac63.0
Applied add-sqr-sqrt63.0
Applied difference-of-squares63.0
Applied add-sqr-sqrt63.0
Applied times-frac50.1
Final simplification6.3
herbie shell --seed 2020148
(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))))