\frac{x + y}{1 - \frac{y}{z}}\begin{array}{l}
\mathbf{if}\;\frac{x + y}{1 - \frac{y}{z}} \leq -3.826021058198235 \cdot 10^{-268} \lor \neg \left(\frac{x + y}{1 - \frac{y}{z}} \leq 0\right):\\
\;\;\;\;\frac{x + y}{1 - \frac{y}{z}}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x + y\right) \cdot \frac{\sqrt[3]{1} \cdot \sqrt[3]{1}}{1 + \frac{\sqrt{y}}{\sqrt{z}}}\right) \cdot \frac{\sqrt[3]{1}}{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.826021058198235e-268)
(not (<= (/ (+ x y) (- 1.0 (/ y z))) 0.0)))
(/ (+ x y) (- 1.0 (/ y z)))
(*
(* (+ x y) (/ (* (cbrt 1.0) (cbrt 1.0)) (+ 1.0 (/ (sqrt y) (sqrt z)))))
(/ (cbrt 1.0) (- 1.0 (/ (sqrt y) (sqrt 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))) <= -3.826021058198235e-268) || !(((x + y) / (1.0 - (y / z))) <= 0.0)) {
tmp = (x + y) / (1.0 - (y / z));
} else {
tmp = ((x + y) * ((cbrt(1.0) * cbrt(1.0)) / (1.0 + (sqrt(y) / sqrt(z))))) * (cbrt(1.0) / (1.0 - (sqrt(y) / sqrt(z))));
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 7.5 |
|---|---|
| Target | 3.9 |
| Herbie | 6.5 |
if (/.f64 (+.f64 x y) (-.f64 1 (/.f64 y z))) < -3.8260210581982348e-268 or -0.0 < (/.f64 (+.f64 x y) (-.f64 1 (/.f64 y z))) Initial program 0.1
if -3.8260210581982348e-268 < (/.f64 (+.f64 x y) (-.f64 1 (/.f64 y z))) < -0.0Initial program 56.6
rmApplied div-inv_binary6456.6
rmApplied add-sqr-sqrt_binary6459.0
Applied add-sqr-sqrt_binary6461.7
Applied times-frac_binary6461.7
Applied *-un-lft-identity_binary6461.7
Applied difference-of-squares_binary6461.7
Applied add-cube-cbrt_binary6461.7
Applied times-frac_binary6460.9
Applied associate-*r*_binary6448.8
Final simplification6.5
herbie shell --seed 2020260
(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))))