\frac{x \cdot y}{z}\begin{array}{l}
\mathbf{if}\;\frac{x \cdot y}{z} = -\infty:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\mathbf{elif}\;\frac{x \cdot y}{z} \le -3.8974497077 \cdot 10^{-315}:\\
\;\;\;\;\frac{x \cdot y}{z}\\
\mathbf{elif}\;\frac{x \cdot y}{z} \le 0.0:\\
\;\;\;\;\frac{x}{\frac{z}{y}}\\
\mathbf{elif}\;\frac{x \cdot y}{z} \le 1.912706477262822 \cdot 10^{267}:\\
\;\;\;\;\frac{x \cdot y}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{z}{y}}\\
\end{array}double code(double x, double y, double z) {
return ((double) (((double) (x * y)) / z));
}
double code(double x, double y, double z) {
double VAR;
if ((((double) (((double) (x * y)) / z)) <= -inf.0)) {
VAR = ((double) (x * ((double) (y / z))));
} else {
double VAR_1;
if ((((double) (((double) (x * y)) / z)) <= -3.8974497077474e-315)) {
VAR_1 = ((double) (((double) (x * y)) / z));
} else {
double VAR_2;
if ((((double) (((double) (x * y)) / z)) <= 0.0)) {
VAR_2 = ((double) (x / ((double) (z / y))));
} else {
double VAR_3;
if ((((double) (((double) (x * y)) / z)) <= 1.912706477262822e+267)) {
VAR_3 = ((double) (((double) (x * y)) / z));
} else {
VAR_3 = ((double) (x / ((double) (z / y))));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.4 |
|---|---|
| Target | 6.4 |
| Herbie | 0.7 |
if (/ (* x y) z) < -inf.0Initial program 64.0
rmApplied *-un-lft-identity64.0
Applied times-frac0.3
Simplified0.3
if -inf.0 < (/ (* x y) z) < -3.8974497077474e-315 or 0.0 < (/ (* x y) z) < 1.912706477262822e+267Initial program 0.5
if -3.8974497077474e-315 < (/ (* x y) z) < 0.0 or 1.912706477262822e+267 < (/ (* x y) z) Initial program 15.9
rmApplied associate-/l*1.2
Final simplification0.7
herbie shell --seed 2020122
(FPCore (x y z)
:name "Diagrams.Solve.Tridiagonal:solveCyclicTriDiagonal from diagrams-solve-0.1, A"
:precision binary64
:herbie-target
(if (< z -4.262230790519429e-138) (/ (* x y) z) (if (< z 1.7042130660650472e-164) (/ x (/ z y)) (* (/ x z) y)))
(/ (* x y) z))