\frac{x \cdot y}{z}\begin{array}{l}
\mathbf{if}\;x \cdot y \le -1.122825852725341740510374421793401440482 \cdot 10^{174}:\\
\;\;\;\;\frac{x}{\frac{z}{y}}\\
\mathbf{elif}\;x \cdot y \le -1.187138743731467067277490994856280614263 \cdot 10^{-188}:\\
\;\;\;\;\frac{x \cdot y}{z}\\
\mathbf{elif}\;x \cdot y \le 2.131572119135182027868179571511010877863 \cdot 10^{-317}:\\
\;\;\;\;x \cdot \frac{y}{z}\\
\mathbf{elif}\;x \cdot y \le 1.379689488789308046398090942052262660111 \cdot 10^{229}:\\
\;\;\;\;\frac{x \cdot y}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{z} \cdot y\\
\end{array}double f(double x, double y, double z) {
double r674269 = x;
double r674270 = y;
double r674271 = r674269 * r674270;
double r674272 = z;
double r674273 = r674271 / r674272;
return r674273;
}
double f(double x, double y, double z) {
double r674274 = x;
double r674275 = y;
double r674276 = r674274 * r674275;
double r674277 = -1.1228258527253417e+174;
bool r674278 = r674276 <= r674277;
double r674279 = z;
double r674280 = r674279 / r674275;
double r674281 = r674274 / r674280;
double r674282 = -1.187138743731467e-188;
bool r674283 = r674276 <= r674282;
double r674284 = r674276 / r674279;
double r674285 = 2.1315721191352e-317;
bool r674286 = r674276 <= r674285;
double r674287 = r674275 / r674279;
double r674288 = r674274 * r674287;
double r674289 = 1.379689488789308e+229;
bool r674290 = r674276 <= r674289;
double r674291 = r674274 / r674279;
double r674292 = r674291 * r674275;
double r674293 = r674290 ? r674284 : r674292;
double r674294 = r674286 ? r674288 : r674293;
double r674295 = r674283 ? r674284 : r674294;
double r674296 = r674278 ? r674281 : r674295;
return r674296;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 6.4 |
|---|---|
| Target | 6.5 |
| Herbie | 0.5 |
if (* x y) < -1.1228258527253417e+174Initial program 22.8
rmApplied associate-/l*1.7
if -1.1228258527253417e+174 < (* x y) < -1.187138743731467e-188 or 2.1315721191352e-317 < (* x y) < 1.379689488789308e+229Initial program 0.3
if -1.187138743731467e-188 < (* x y) < 2.1315721191352e-317Initial program 13.5
rmApplied *-un-lft-identity13.5
Applied times-frac0.6
Simplified0.6
if 1.379689488789308e+229 < (* x y) Initial program 31.3
rmApplied associate-/l*0.3
rmApplied associate-/r/0.8
Final simplification0.5
herbie shell --seed 2019354
(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))