\frac{x + y}{1 - \frac{y}{z}}\begin{array}{l}
\mathbf{if}\;\frac{x + y}{1 - \frac{y}{z}} \leq -3.253524157806006 \cdot 10^{-277} \lor \neg \left(\frac{x + y}{1 - \frac{y}{z}} \leq 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}(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))) -3.253524157806006e-277)
(not (<= (/ (+ x y) (- 1.0 (/ y z))) 0.0)))
(/ (+ x y) (- 1.0 (/ y z)))
(*
(/ (sqrt (+ x y)) (+ (sqrt 1.0) (/ (sqrt y) (sqrt z))))
(/ (sqrt (+ x y)) (- (sqrt 1.0) (/ (sqrt y) (sqrt z)))))))double code(double x, double y, double z) {
return (((double) (x + y)) / ((double) (1.0 - (y / z))));
}
double code(double x, double y, double z) {
double tmp;
if ((((((double) (x + y)) / ((double) (1.0 - (y / z)))) <= -3.253524157806006e-277) || !((((double) (x + y)) / ((double) (1.0 - (y / z)))) <= 0.0))) {
tmp = (((double) (x + y)) / ((double) (1.0 - (y / z))));
} else {
tmp = ((double) ((((double) sqrt(((double) (x + y)))) / ((double) (((double) sqrt(1.0)) + (((double) sqrt(y)) / ((double) sqrt(z)))))) * (((double) sqrt(((double) (x + y)))) / ((double) (((double) sqrt(1.0)) - (((double) sqrt(y)) / ((double) sqrt(z))))))));
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 7.8 |
|---|---|
| Target | 4.0 |
| Herbie | 6.6 |
if (/.f64 (+.f64 x y) (-.f64 1.0 (/.f64 y z))) < -3.2535241578060062e-277 or 0.0 < (/.f64 (+.f64 x y) (-.f64 1.0 (/.f64 y z))) Initial program 0.1
if -3.2535241578060062e-277 < (/.f64 (+.f64 x y) (-.f64 1.0 (/.f64 y z))) < 0.0Initial program 57.8
rmApplied add-sqr-sqrt_binary6460.1
Applied add-sqr-sqrt_binary6462.1
Applied times-frac_binary6462.1
Applied add-sqr-sqrt_binary6462.1
Applied difference-of-squares_binary6462.1
Applied add-sqr-sqrt_binary6462.2
Applied times-frac_binary6449.2
Final simplification6.6
herbie shell --seed 2020210
(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))))