\frac{x \cdot y}{z}\begin{array}{l}
\mathbf{if}\;x \cdot y \le -2.6967839648114804 \cdot 10^{268}:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\mathbf{elif}\;x \cdot y \le -7.97061052598067788 \cdot 10^{-294}:\\
\;\;\;\;\frac{x \cdot y}{z}\\
\mathbf{elif}\;x \cdot y \le 6.0589118027638559 \cdot 10^{-195}:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\mathbf{elif}\;x \cdot y \le 1.29204454885475858 \cdot 10^{212}:\\
\;\;\;\;\frac{x \cdot y}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{z}{y}}\\
\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) <= -2.6967839648114804e+268)) {
temp = (x * (y / z));
} else {
double temp_1;
if (((x * y) <= -7.970610525980678e-294)) {
temp_1 = ((x * y) / z);
} else {
double temp_2;
if (((x * y) <= 6.058911802763856e-195)) {
temp_2 = (x * (y / z));
} else {
double temp_3;
if (((x * y) <= 1.2920445488547586e+212)) {
temp_3 = ((x * y) / z);
} else {
temp_3 = (x / (z / y));
}
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 | 5.8 |
|---|---|
| Target | 5.9 |
| Herbie | 0.3 |
if (* x y) < -2.6967839648114804e+268 or -7.970610525980678e-294 < (* x y) < 6.058911802763856e-195Initial program 15.3
rmApplied *-un-lft-identity15.3
Applied times-frac0.5
Simplified0.5
if -2.6967839648114804e+268 < (* x y) < -7.970610525980678e-294 or 6.058911802763856e-195 < (* x y) < 1.2920445488547586e+212Initial program 0.2
if 1.2920445488547586e+212 < (* x y) Initial program 28.8
rmApplied associate-/l*1.2
Final simplification0.3
herbie shell --seed 2020057 +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))