\frac{x + y}{1 - \frac{y}{z}}\begin{array}{l}
\mathbf{if}\;\frac{x + y}{1 - \frac{y}{z}} \leq -9.261366165399573 \cdot 10^{-284}:\\
\;\;\;\;\frac{x + y}{1 - \frac{y}{z}}\\
\mathbf{elif}\;\frac{x + y}{1 - \frac{y}{z}} \leq 0:\\
\;\;\;\;\frac{\sqrt{x + y}}{1 + \frac{\sqrt{y}}{\sqrt{z}}} \cdot \frac{\sqrt{x + y}}{1 - \frac{\sqrt{y}}{\sqrt{z}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x + y}{1 - \frac{y}{z}}\\
\end{array}(FPCore (x y z) :precision binary64 (/ (+ x y) (- 1.0 (/ y z))))
(FPCore (x y z)
:precision binary64
(if (<= (/ (+ x y) (- 1.0 (/ y z))) -9.261366165399573e-284)
(/ (+ x y) (- 1.0 (/ y z)))
(if (<= (/ (+ x y) (- 1.0 (/ y z))) 0.0)
(*
(/ (sqrt (+ x y)) (+ 1.0 (/ (sqrt y) (sqrt z))))
(/ (sqrt (+ x y)) (- 1.0 (/ (sqrt y) (sqrt z)))))
(/ (+ x y) (- 1.0 (/ y z))))))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))) <= -9.261366165399573e-284) {
tmp = (x + y) / (1.0 - (y / z));
} else if (((x + y) / (1.0 - (y / z))) <= 0.0) {
tmp = (sqrt(x + y) / (1.0 + (sqrt(y) / sqrt(z)))) * (sqrt(x + y) / (1.0 - (sqrt(y) / sqrt(z))));
} else {
tmp = (x + y) / (1.0 - (y / z));
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 7.7 |
|---|---|
| Target | 3.8 |
| Herbie | 6.5 |
if (/.f64 (+.f64 x y) (-.f64 1 (/.f64 y z))) < -9.2613661653995734e-284 or -0.0 < (/.f64 (+.f64 x y) (-.f64 1 (/.f64 y z))) Initial program 0.1
if -9.2613661653995734e-284 < (/.f64 (+.f64 x y) (-.f64 1 (/.f64 y z))) < -0.0Initial program 58.6
rmApplied add-sqr-sqrt_binary6461.0
Applied add-sqr-sqrt_binary6462.5
Applied times-frac_binary6462.5
Applied add-sqr-sqrt_binary6462.5
Applied difference-of-squares_binary6462.5
Applied add-sqr-sqrt_binary6462.7
Applied times-frac_binary6449.2
Simplified49.2
Simplified49.2
Final simplification6.5
herbie shell --seed 2020273
(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))))