\frac{x \cdot y}{z}\begin{array}{l}
\mathbf{if}\;x \cdot y \le -7.725807749934587 \cdot 10^{96}:\\
\;\;\;\;\frac{x}{z} \cdot y\\
\mathbf{elif}\;x \cdot y \le -3.0807367441454416 \cdot 10^{-282}:\\
\;\;\;\;\frac{x \cdot y}{z}\\
\mathbf{elif}\;x \cdot y \le -0.0:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\mathbf{elif}\;x \cdot y \le 5.31151602867187701 \cdot 10^{255}:\\
\;\;\;\;\frac{x \cdot y}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{1}{\frac{x}{\frac{z}{y}}}}\\
\end{array}double f(double x, double y, double z) {
double r820741 = x;
double r820742 = y;
double r820743 = r820741 * r820742;
double r820744 = z;
double r820745 = r820743 / r820744;
return r820745;
}
double f(double x, double y, double z) {
double r820746 = x;
double r820747 = y;
double r820748 = r820746 * r820747;
double r820749 = -7.725807749934587e+96;
bool r820750 = r820748 <= r820749;
double r820751 = z;
double r820752 = r820746 / r820751;
double r820753 = r820752 * r820747;
double r820754 = -3.0807367441454416e-282;
bool r820755 = r820748 <= r820754;
double r820756 = r820748 / r820751;
double r820757 = -0.0;
bool r820758 = r820748 <= r820757;
double r820759 = r820747 / r820751;
double r820760 = r820746 * r820759;
double r820761 = 5.311516028671877e+255;
bool r820762 = r820748 <= r820761;
double r820763 = 1.0;
double r820764 = r820751 / r820747;
double r820765 = r820746 / r820764;
double r820766 = r820763 / r820765;
double r820767 = r820763 / r820766;
double r820768 = r820762 ? r820756 : r820767;
double r820769 = r820758 ? r820760 : r820768;
double r820770 = r820755 ? r820756 : r820769;
double r820771 = r820750 ? r820753 : r820770;
return r820771;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.3 |
|---|---|
| Target | 6.4 |
| Herbie | 0.8 |
if (* x y) < -7.725807749934587e+96Initial program 14.6
rmApplied associate-/l*3.3
rmApplied associate-/r/4.1
if -7.725807749934587e+96 < (* x y) < -3.0807367441454416e-282 or -0.0 < (* x y) < 5.311516028671877e+255Initial program 0.4
if -3.0807367441454416e-282 < (* x y) < -0.0Initial program 17.6
rmApplied *-un-lft-identity17.6
Applied times-frac0.1
Simplified0.1
if 5.311516028671877e+255 < (* x y) Initial program 40.0
rmApplied associate-/l*0.3
rmApplied clear-num0.3
rmApplied clear-num0.4
Final simplification0.8
herbie shell --seed 2020060
(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))