\frac{x \cdot y}{z}\begin{array}{l}
\mathbf{if}\;x \cdot y \leq -1.2583117062414792 \cdot 10^{+189}:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\mathbf{elif}\;x \cdot y \leq -3.465479670824521 \cdot 10^{-205} \lor \neg \left(x \cdot y \leq 2.9882567291840723 \cdot 10^{-243}\right) \land x \cdot y \leq 7.971146091977197 \cdot 10^{+161}:\\
\;\;\;\;\frac{x \cdot y}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{z}{y}}\\
\end{array}(FPCore (x y z) :precision binary64 (/ (* x y) z))
(FPCore (x y z)
:precision binary64
(if (<= (* x y) -1.2583117062414792e+189)
(* x (/ y z))
(if (or (<= (* x y) -3.465479670824521e-205)
(and (not (<= (* x y) 2.9882567291840723e-243))
(<= (* x y) 7.971146091977197e+161)))
(/ (* x y) z)
(/ x (/ z y)))))double code(double x, double y, double z) {
return (x * y) / z;
}
double code(double x, double y, double z) {
double tmp;
if ((x * y) <= -1.2583117062414792e+189) {
tmp = x * (y / z);
} else if (((x * y) <= -3.465479670824521e-205) || (!((x * y) <= 2.9882567291840723e-243) && ((x * y) <= 7.971146091977197e+161))) {
tmp = (x * y) / z;
} else {
tmp = x / (z / y);
}
return tmp;
}









Bits error versus x









Bits error versus y









Bits error versus z
Results
| Original | 6.4 |
|---|---|
| Target | 6.4 |
| Herbie | 0.6 |
| Alternative 1 | |
|---|---|
| Error | 6.1 |
| Cost | 1283 |
| Alternative 2 | |
|---|---|
| Error | 6.2 |
| Cost | 320 |
| Alternative 3 | |
|---|---|
| Error | 6.1 |
| Cost | 320 |
| Alternative 4 | |
|---|---|
| Error | 49.5 |
| Cost | 64 |
| Alternative 5 | |
|---|---|
| Error | 61.8 |
| Cost | 64 |

if (*.f64 x y) < -1.25831170624147922e189Initial program 24.5
rmApplied *-un-lft-identity_binary64_2053824.5
Applied times-frac_binary64_205441.9
Simplified1.9
Simplified1.9
if -1.25831170624147922e189 < (*.f64 x y) < -3.46547967082452108e-205 or 2.9882567291840723e-243 < (*.f64 x y) < 7.97114609197719727e161Initial program 0.2
Simplified0.2
if -3.46547967082452108e-205 < (*.f64 x y) < 2.9882567291840723e-243 or 7.97114609197719727e161 < (*.f64 x y) Initial program 14.0
rmApplied associate-/l*_binary64_204830.9
Simplified0.9
Final simplification0.6
herbie shell --seed 2021044
(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))