\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}\;y \le -9.1777242374334205 \cdot 10^{59}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{1}{\frac{z \cdot \frac{c}{y}}{x}}\right) - 4 \cdot \frac{a \cdot t}{c}\\
\mathbf{elif}\;y \le 1.2889311644196529 \cdot 10^{-7}:\\
\;\;\;\;\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}\;y \le 8.81398444710026263 \cdot 10^{115}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{\frac{x}{z \cdot c}}{\frac{1}{y}}\right) - 4 \cdot \frac{a \cdot t}{c}\\
\mathbf{elif}\;y \le 5.87751849036165333 \cdot 10^{168}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \left(a \cdot \frac{t}{c}\right)\\
\mathbf{elif}\;y \le 6.1051022150981975 \cdot 10^{253}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + \left(9 \cdot \frac{x}{z}\right) \cdot \frac{y}{c}\right) - 4 \cdot \frac{a \cdot t}{c}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \left(\left(\sqrt[3]{\frac{x}{\frac{z \cdot c}{y}}} \cdot \sqrt[3]{\frac{x}{\frac{z \cdot c}{y}}}\right) \cdot \sqrt[3]{\frac{x}{\frac{z \cdot c}{y}}}\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 ((y <= -9.17772423743342e+59)) {
VAR = (((b / (z * c)) + (9.0 * (1.0 / ((z * (c / y)) / x)))) - (4.0 * ((a * t) / c)));
} else {
double VAR_1;
if ((y <= 1.2889311644196529e-07)) {
VAR_1 = (((b / (z * c)) + (9.0 * ((x * y) / (z * c)))) - (4.0 * (a / (c / t))));
} else {
double VAR_2;
if ((y <= 8.813984447100263e+115)) {
VAR_2 = (((b / (z * c)) + (9.0 * ((x / (z * c)) / (1.0 / y)))) - (4.0 * ((a * t) / c)));
} else {
double VAR_3;
if ((y <= 5.877518490361653e+168)) {
VAR_3 = (((b / (z * c)) + (9.0 * ((x * y) / (z * c)))) - (4.0 * (a * (t / c))));
} else {
double VAR_4;
if ((y <= 6.105102215098197e+253)) {
VAR_4 = (((b / (z * c)) + ((9.0 * (x / z)) * (y / c))) - (4.0 * ((a * t) / c)));
} else {
VAR_4 = (((b / (z * c)) + (9.0 * ((cbrt((x / ((z * c) / y))) * cbrt((x / ((z * c) / y)))) * cbrt((x / ((z * c) / y)))))) - (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.6 |
|---|---|
| Target | 14.4 |
| Herbie | 10.4 |
if y < -9.17772423743342e+59Initial program 26.8
Taylor expanded around 0 18.8
rmApplied associate-/l*15.1
rmApplied *-un-lft-identity15.1
Applied times-frac14.9
Simplified14.9
rmApplied clear-num15.0
if -9.17772423743342e+59 < y < 1.2889311644196529e-07Initial program 17.0
Taylor expanded around 0 8.0
rmApplied associate-/l*7.8
if 1.2889311644196529e-07 < y < 8.813984447100263e+115Initial program 20.0
Taylor expanded around 0 12.0
rmApplied associate-/l*9.4
rmApplied div-inv9.4
Applied associate-/r*9.4
if 8.813984447100263e+115 < y < 5.877518490361653e+168Initial program 27.6
Taylor expanded around 0 17.2
rmApplied *-un-lft-identity17.2
Applied times-frac17.0
Simplified17.0
if 5.877518490361653e+168 < y < 6.105102215098197e+253Initial program 28.2
Taylor expanded around 0 20.0
rmApplied times-frac14.2
Applied associate-*r*14.2
if 6.105102215098197e+253 < y Initial program 33.1
Taylor expanded around 0 23.9
rmApplied associate-/l*21.5
rmApplied add-cube-cbrt21.7
Final simplification10.4
herbie shell --seed 2020102
(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)))