\frac{x \cdot y}{z}\begin{array}{l}
\mathbf{if}\;x \cdot y = -\infty:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\mathbf{elif}\;x \cdot y \le -2.7755259961563269 \cdot 10^{-232}:\\
\;\;\;\;\frac{x \cdot y}{z}\\
\mathbf{elif}\;x \cdot y \le 6.23637277415903404 \cdot 10^{-305}:\\
\;\;\;\;\frac{x}{\frac{z}{y}}\\
\mathbf{elif}\;x \cdot y \le 2.791679912717515 \cdot 10^{186}:\\
\;\;\;\;\frac{x \cdot y}{z}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\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) (x * y)) <= -inf.0)) {
VAR = ((double) (x * ((double) (y / z))));
} else {
double VAR_1;
if ((((double) (x * y)) <= -2.775525996156327e-232)) {
VAR_1 = ((double) (((double) (x * y)) / z));
} else {
double VAR_2;
if ((((double) (x * y)) <= 6.236372774159034e-305)) {
VAR_2 = ((double) (x / ((double) (z / y))));
} else {
double VAR_3;
if ((((double) (x * y)) <= 2.7916799127175147e+186)) {
VAR_3 = ((double) (((double) (x * y)) / z));
} else {
VAR_3 = ((double) (x * ((double) (y / z))));
}
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.0 |
| Herbie | 0.3 |
if (* x y) < -inf.0 or 2.7916799127175147e+186 < (* x y) Initial program 36.0
rmApplied *-un-lft-identity36.0
Applied times-frac0.8
Simplified0.8
if -inf.0 < (* x y) < -2.775525996156327e-232 or 6.236372774159034e-305 < (* x y) < 2.7916799127175147e+186Initial program 0.2
if -2.775525996156327e-232 < (* x y) < 6.236372774159034e-305Initial program 14.7
rmApplied associate-/l*0.2
Final simplification0.3
herbie shell --seed 2020123 +o rules:numerics
(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))