\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}\;c \le -1.8461967142404801 \cdot 10^{167}:\\
\;\;\;\;\left(\frac{1}{z} \cdot \frac{b}{c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \left(\left(\left(\sqrt[3]{a} \cdot \sqrt[3]{a}\right) \cdot \left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)\right) \cdot \frac{\sqrt[3]{a}}{\frac{c}{\sqrt[3]{t}}}\right)\\
\mathbf{elif}\;c \le -3.5147769750141271 \cdot 10^{30}:\\
\;\;\;\;\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}\\
\mathbf{elif}\;c \le -6.3887029057961972 \cdot 10^{-159}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \frac{\frac{a}{c}}{\frac{1}{t}}\\
\mathbf{elif}\;c \le 2.11190068261123611 \cdot 10^{-110}:\\
\;\;\;\;\frac{\left(\left(x \cdot 9\right) \cdot y - \left(z \cdot 4\right) \cdot \left(t \cdot a\right)\right) + b}{z \cdot c}\\
\mathbf{elif}\;c \le 8.6607696898741896 \cdot 10^{-82}:\\
\;\;\;\;\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{elif}\;c \le 1.97963596605295414 \cdot 10^{180}:\\
\;\;\;\;\left(\frac{1}{z} \cdot \frac{b}{c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \left(\left(\left(\sqrt[3]{a} \cdot \sqrt[3]{a}\right) \cdot \left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)\right) \cdot \frac{\sqrt[3]{a}}{\frac{c}{\sqrt[3]{t}}}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\frac{z \cdot c}{y}}\right) - 4 \cdot \frac{a}{\frac{c}{t}}\\
\end{array}double code(double x, double y, double z, double t, double a, double b, double c) {
return ((double) (((double) (((double) (((double) (((double) (x * 9.0)) * y)) - ((double) (((double) (((double) (z * 4.0)) * t)) * a)))) + b)) / ((double) (z * c))));
}
double code(double x, double y, double z, double t, double a, double b, double c) {
double VAR;
if ((c <= -1.84619671424048e+167)) {
VAR = ((double) (((double) (((double) (((double) (1.0 / z)) * ((double) (b / c)))) + ((double) (9.0 * ((double) (((double) (x * y)) / ((double) (z * c)))))))) - ((double) (4.0 * ((double) (((double) (((double) (((double) cbrt(a)) * ((double) cbrt(a)))) * ((double) (((double) cbrt(t)) * ((double) cbrt(t)))))) * ((double) (((double) cbrt(a)) / ((double) (c / ((double) cbrt(t))))))))))));
} else {
double VAR_1;
if ((c <= -3.514776975014127e+30)) {
VAR_1 = ((double) (((double) (((double) (b / ((double) (z * c)))) + ((double) (9.0 * ((double) (((double) (x / z)) * ((double) (y / c)))))))) - ((double) (4.0 * ((double) (((double) (a * t)) / c))))));
} else {
double VAR_2;
if ((c <= -6.388702905796197e-159)) {
VAR_2 = ((double) (((double) (((double) (b / ((double) (z * c)))) + ((double) (9.0 * ((double) (((double) (x * y)) / ((double) (z * c)))))))) - ((double) (4.0 * ((double) (((double) (a / c)) / ((double) (1.0 / t))))))));
} else {
double VAR_3;
if ((c <= 2.111900682611236e-110)) {
VAR_3 = ((double) (((double) (((double) (((double) (((double) (x * 9.0)) * y)) - ((double) (((double) (z * 4.0)) * ((double) (t * a)))))) + b)) / ((double) (z * c))));
} else {
double VAR_4;
if ((c <= 8.66076968987419e-82)) {
VAR_4 = ((double) (((double) (((double) (b / ((double) (z * c)))) + ((double) (((double) (9.0 * ((double) (x / z)))) * ((double) (y / c)))))) - ((double) (4.0 * ((double) (((double) (a * t)) / c))))));
} else {
double VAR_5;
if ((c <= 1.9796359660529541e+180)) {
VAR_5 = ((double) (((double) (((double) (((double) (1.0 / z)) * ((double) (b / c)))) + ((double) (9.0 * ((double) (((double) (x * y)) / ((double) (z * c)))))))) - ((double) (4.0 * ((double) (((double) (((double) (((double) cbrt(a)) * ((double) cbrt(a)))) * ((double) (((double) cbrt(t)) * ((double) cbrt(t)))))) * ((double) (((double) cbrt(a)) / ((double) (c / ((double) cbrt(t))))))))))));
} else {
VAR_5 = ((double) (((double) (((double) (b / ((double) (z * c)))) + ((double) (9.0 * ((double) (x / ((double) (((double) (z * c)) / y)))))))) - ((double) (4.0 * ((double) (a / ((double) (c / t))))))));
}
VAR_4 = VAR_5;
}
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.0 |
|---|---|
| Target | 14.0 |
| Herbie | 9.5 |
if c < -1.84619671424048e+167 or 8.66076968987419e-82 < c < 1.9796359660529541e+180Initial program 20.7
Taylor expanded around 0 12.7
rmApplied associate-/l*9.8
rmApplied add-cube-cbrt10.1
Applied *-un-lft-identity10.1
Applied times-frac10.1
Applied add-cube-cbrt10.2
Applied times-frac9.6
Simplified9.6
rmApplied *-un-lft-identity9.6
Applied times-frac8.6
if -1.84619671424048e+167 < c < -3.514776975014127e+30Initial program 20.4
Taylor expanded around 0 10.9
rmApplied times-frac9.9
if -3.514776975014127e+30 < c < -6.388702905796197e-159Initial program 14.8
Taylor expanded around 0 5.1
rmApplied associate-/l*7.8
rmApplied div-inv7.8
Applied associate-/r*8.3
if -6.388702905796197e-159 < c < 2.111900682611236e-110Initial program 13.9
rmApplied associate-*l*9.1
if 2.111900682611236e-110 < c < 8.66076968987419e-82Initial program 17.7
Taylor expanded around 0 7.0
rmApplied times-frac10.1
Applied associate-*r*10.1
if 1.9796359660529541e+180 < c Initial program 27.8
Taylor expanded around 0 18.8
rmApplied associate-/l*15.7
rmApplied associate-/l*12.8
Final simplification9.5
herbie shell --seed 2020122
(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)))