\frac{\left(\left(x \cdot 9\right) \cdot y - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b}{z \cdot c}\begin{array}{l}
\mathbf{if}\;x \cdot 9 \le -8.97528476763240328 \cdot 10^{118}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\frac{z \cdot c}{y}}\right) - 4 \cdot \frac{a \cdot t}{c}\\
\mathbf{elif}\;x \cdot 9 \le -4.3405214492461963 \cdot 10^{-125}:\\
\;\;\;\;\left(\frac{\frac{b}{c}}{z} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \frac{a \cdot t}{c}\\
\mathbf{elif}\;x \cdot 9 \le -3.68290799341614394 \cdot 10^{-263}:\\
\;\;\;\;\left(\frac{\sqrt[3]{b} \cdot \sqrt[3]{b}}{\frac{z \cdot c}{\sqrt[3]{b}}} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \frac{a \cdot t}{c}\\
\mathbf{elif}\;x \cdot 9 \le 6.53549206299654139 \cdot 10^{-125}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \frac{a}{\frac{c}{t}}\\
\mathbf{elif}\;x \cdot 9 \le 2.67125869365298492 \cdot 10^{83}:\\
\;\;\;\;\frac{1}{z} \cdot \frac{\left(\left(x \cdot 9\right) \cdot y - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b}{c}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \left(\frac{x}{z} \cdot \frac{y}{c}\right)\right) - 4 \cdot \frac{a \cdot t}{c}\\
\end{array}double code(double x, double y, double z, double t, double a, double b, double c) {
return (((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / (z * c));
}
double code(double x, double y, double z, double t, double a, double b, double c) {
double VAR;
if (((x * 9.0) <= -8.975284767632403e+118)) {
VAR = (((b / (z * c)) + (9.0 * (x / ((z * c) / y)))) - (4.0 * ((a * t) / c)));
} else {
double VAR_1;
if (((x * 9.0) <= -4.340521449246196e-125)) {
VAR_1 = ((((b / c) / z) + (9.0 * ((x * y) / (z * c)))) - (4.0 * ((a * t) / c)));
} else {
double VAR_2;
if (((x * 9.0) <= -3.682907993416144e-263)) {
VAR_2 = ((((cbrt(b) * cbrt(b)) / ((z * c) / cbrt(b))) + (9.0 * ((x * y) / (z * c)))) - (4.0 * ((a * t) / c)));
} else {
double VAR_3;
if (((x * 9.0) <= 6.535492062996541e-125)) {
VAR_3 = (((b / (z * c)) + (9.0 * ((x * y) / (z * c)))) - (4.0 * (a / (c / t))));
} else {
double VAR_4;
if (((x * 9.0) <= 2.671258693652985e+83)) {
VAR_4 = ((1.0 / z) * (((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / c));
} else {
VAR_4 = (((b / (z * c)) + (9.0 * ((x / z) * (y / c)))) - (4.0 * ((a * t) / c)));
}
VAR_3 = VAR_4;
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
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




Bits error versus c
Results
| Original | 20.4 |
|---|---|
| Target | 14.6 |
| Herbie | 12.3 |
if (* x 9.0) < -8.975284767632403e+118Initial program 28.4
Taylor expanded around 0 20.8
rmApplied associate-/l*12.7
if -8.975284767632403e+118 < (* x 9.0) < -4.340521449246196e-125Initial program 18.4
Taylor expanded around 0 10.5
rmApplied *-un-lft-identity10.5
Applied times-frac11.3
rmApplied associate-*l/11.2
Simplified11.2
if -4.340521449246196e-125 < (* x 9.0) < -3.682907993416144e-263Initial program 16.3
Taylor expanded around 0 7.3
rmApplied add-cube-cbrt7.6
Applied associate-/l*7.6
if -3.682907993416144e-263 < (* x 9.0) < 6.535492062996541e-125Initial program 16.6
Taylor expanded around 0 7.8
rmApplied associate-/l*7.5
if 6.535492062996541e-125 < (* x 9.0) < 2.671258693652985e+83Initial program 18.5
rmApplied *-un-lft-identity18.5
Applied times-frac18.0
if 2.671258693652985e+83 < (* x 9.0) Initial program 27.6
Taylor expanded around 0 19.2
rmApplied times-frac17.6
Final simplification12.3
herbie shell --seed 2020100
(FPCore (x y z t a b c)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, J"
:precision binary64
:herbie-target
(if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) -1.1001567408041051e-171) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) -0.0) (/ (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) z) c) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 1.1708877911747488e-53) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 2.876823679546137e+130) (- (+ (* (* 9 (/ y c)) (/ x z)) (/ b (* c z))) (* 4 (/ (* a t) c))) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 1.3838515042456319e+158) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (- (+ (* 9 (* (/ y (* c z)) x)) (/ b (* c z))) (* 4 (/ (* a t) c))))))))
(/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)))