\frac{x \cdot y}{z}\begin{array}{l}
\mathbf{if}\;x \cdot y \le -1.21374163384962988 \cdot 10^{271}:\\
\;\;\;\;\frac{x}{\frac{z}{y}}\\
\mathbf{elif}\;x \cdot y \le -3.4666700231195389 \cdot 10^{-139}:\\
\;\;\;\;1 \cdot \frac{x \cdot y}{z}\\
\mathbf{elif}\;x \cdot y \le 0.0:\\
\;\;\;\;\frac{x}{\frac{z}{y}}\\
\mathbf{elif}\;x \cdot y \le 1.1437129577976631 \cdot 10^{286}:\\
\;\;\;\;1 \cdot \frac{x \cdot y}{z}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\end{array}double code(double x, double y, double z) {
return ((x * y) / z);
}
double code(double x, double y, double z) {
double temp;
if (((x * y) <= -1.2137416338496299e+271)) {
temp = (x / (z / y));
} else {
double temp_1;
if (((x * y) <= -3.466670023119539e-139)) {
temp_1 = (1.0 * ((x * y) / z));
} else {
double temp_2;
if (((x * y) <= 0.0)) {
temp_2 = (x / (z / y));
} else {
double temp_3;
if (((x * y) <= 1.1437129577976631e+286)) {
temp_3 = (1.0 * ((x * y) / z));
} else {
temp_3 = (x * (y / z));
}
temp_2 = temp_3;
}
temp_1 = temp_2;
}
temp = temp_1;
}
return temp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.2 |
|---|---|
| Target | 6.2 |
| Herbie | 0.5 |
if (* x y) < -1.2137416338496299e+271 or -3.466670023119539e-139 < (* x y) < 0.0Initial program 15.3
rmApplied associate-/l*1.0
if -1.2137416338496299e+271 < (* x y) < -3.466670023119539e-139 or 0.0 < (* x y) < 1.1437129577976631e+286Initial program 0.3
rmApplied clear-num0.7
rmApplied *-un-lft-identity0.7
Applied add-cube-cbrt0.7
Applied times-frac0.7
Simplified0.7
Simplified0.3
if 1.1437129577976631e+286 < (* x y) Initial program 52.9
rmApplied *-un-lft-identity52.9
Applied times-frac0.2
Simplified0.2
Final simplification0.5
herbie shell --seed 2020066 +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))