\frac{x \cdot \left(y - z\right)}{y}\begin{array}{l}
\mathbf{if}\;\frac{x \cdot \left(y - z\right)}{y} \leq -4.317300523670314 \cdot 10^{+292}:\\
\;\;\;\;x \cdot \left(1 - \frac{z}{y}\right)\\
\mathbf{elif}\;\frac{x \cdot \left(y - z\right)}{y} \leq -8066392121713.778 \lor \neg \left(\frac{x \cdot \left(y - z\right)}{y} \leq 70307706833107.55\right) \land \frac{x \cdot \left(y - z\right)}{y} \leq 2.2205050860855823 \cdot 10^{+306}:\\
\;\;\;\;x - \frac{x \cdot z}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{y}{y - z}}\\
\end{array}(FPCore (x y z) :precision binary64 (/ (* x (- y z)) y))
(FPCore (x y z)
:precision binary64
(if (<= (/ (* x (- y z)) y) -4.317300523670314e+292)
(* x (- 1.0 (/ z y)))
(if (or (<= (/ (* x (- y z)) y) -8066392121713.778)
(and (not (<= (/ (* x (- y z)) y) 70307706833107.55))
(<= (/ (* x (- y z)) y) 2.2205050860855823e+306)))
(- x (/ (* x z) y))
(/ x (/ y (- y z))))))double code(double x, double y, double z) {
return (x * (y - z)) / y;
}
double code(double x, double y, double z) {
double tmp;
if (((x * (y - z)) / y) <= -4.317300523670314e+292) {
tmp = x * (1.0 - (z / y));
} else if ((((x * (y - z)) / y) <= -8066392121713.778) || (!(((x * (y - z)) / y) <= 70307706833107.55) && (((x * (y - z)) / y) <= 2.2205050860855823e+306))) {
tmp = x - ((x * z) / y);
} else {
tmp = x / (y / (y - z));
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 12.4 |
|---|---|
| Target | 3.1 |
| Herbie | 0.2 |
if (/.f64 (*.f64 x (-.f64 y z)) y) < -4.31730052367031402e292Initial program 57.5
rmApplied *-un-lft-identity_binary6457.5
Applied times-frac_binary641.5
Simplified1.5
rmApplied *-un-lft-identity_binary641.5
Applied *-un-lft-identity_binary641.5
Applied times-frac_binary641.5
Simplified1.5
Simplified1.5
if -4.31730052367031402e292 < (/.f64 (*.f64 x (-.f64 y z)) y) < -8066392121713.77832 or 70307706833107.547 < (/.f64 (*.f64 x (-.f64 y z)) y) < 2.22050508608558229e306Initial program 0.2
Taylor expanded around 0 0.2
if -8066392121713.77832 < (/.f64 (*.f64 x (-.f64 y z)) y) < 70307706833107.547 or 2.22050508608558229e306 < (/.f64 (*.f64 x (-.f64 y z)) y) Initial program 14.8
rmApplied sub-neg_binary6414.8
Applied distribute-rgt-in_binary6414.8
Simplified14.8
Simplified14.8
rmApplied distribute-rgt-neg-in_binary6414.8
Applied distribute-lft-out_binary6414.8
Applied associate-/l*_binary640.1
Final simplification0.2
herbie shell --seed 2021110
(FPCore (x y z)
:name "Diagrams.Backend.Cairo.Internal:setTexture from diagrams-cairo-1.3.0.3"
:precision binary64
:herbie-target
(if (< z -2.060202331921739e+104) (- x (/ (* z x) y)) (if (< z 1.6939766013828526e+213) (/ x (/ y (- y z))) (* (- y z) (/ x y))))
(/ (* x (- y z)) y))