\frac{x + \frac{y \cdot z}{t}}{\left(a + 1\right) + \frac{y \cdot b}{t}}\begin{array}{l}
\mathbf{if}\;t \le -7.41591302982763643 \cdot 10^{99} \lor \neg \left(t \le 2.65209413434131747 \cdot 10^{70}\right):\\
\;\;\;\;\frac{x + \frac{y}{\sqrt[3]{t} \cdot \sqrt[3]{t}} \cdot \frac{z}{\sqrt[3]{t}}}{\left(a + 1\right) + y \cdot \frac{b}{t}}\\
\mathbf{else}:\\
\;\;\;\;\left(x + \frac{y \cdot z}{t}\right) \cdot \frac{1}{\left(a + 1\right) + \frac{y \cdot b}{t}}\\
\end{array}double code(double x, double y, double z, double t, double a, double b) {
return ((double) (((double) (x + ((double) (((double) (y * z)) / t)))) / ((double) (((double) (a + 1.0)) + ((double) (((double) (y * b)) / t))))));
}
double code(double x, double y, double z, double t, double a, double b) {
double VAR;
if (((t <= -7.415913029827636e+99) || !(t <= 2.6520941343413175e+70))) {
VAR = ((double) (((double) (x + ((double) (((double) (y / ((double) (((double) cbrt(t)) * ((double) cbrt(t)))))) * ((double) (z / ((double) cbrt(t)))))))) / ((double) (((double) (a + 1.0)) + ((double) (y * ((double) (b / t))))))));
} else {
VAR = ((double) (((double) (x + ((double) (((double) (y * z)) / t)))) * ((double) (1.0 / ((double) (((double) (a + 1.0)) + ((double) (((double) (y * b)) / t))))))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b
Results
| Original | 17.2 |
|---|---|
| Target | 14.0 |
| Herbie | 13.8 |
if t < -7.415913029827636e+99 or 2.6520941343413175e+70 < t Initial program 12.5
rmApplied add-cube-cbrt12.6
Applied times-frac8.3
rmApplied *-un-lft-identity8.3
Applied times-frac3.1
Simplified3.1
if -7.415913029827636e+99 < t < 2.6520941343413175e+70Initial program 19.9
rmApplied div-inv20.0
Final simplification13.8
herbie shell --seed 2020124
(FPCore (x y z t a b)
:name "Diagrams.Solve.Tridiagonal:solveCyclicTriDiagonal from diagrams-solve-0.1, B"
:precision binary64
:herbie-target
(if (< t -1.3659085366310088e-271) (* 1.0 (* (+ x (* (/ y t) z)) (/ 1.0 (+ (+ a 1.0) (* (/ y t) b))))) (if (< t 3.036967103737246e-130) (/ z b) (* 1.0 (* (+ x (* (/ y t) z)) (/ 1.0 (+ (+ a 1.0) (* (/ y t) b)))))))
(/ (+ x (/ (* y z) t)) (+ (+ a 1.0) (/ (* y b) t))))