\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}\;\left(x \cdot 9\right) \cdot y \le -5.2031824371176434 \cdot 10^{281}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \left(\frac{y}{c} \cdot \frac{x}{z}\right)\right) - 4 \cdot \left(\left(a \cdot \left(\sqrt[3]{t} \cdot \left(\sqrt[3]{\sqrt[3]{t}} \cdot \sqrt[3]{\sqrt[3]{t}}\right)\right)\right) \cdot \left(\sqrt[3]{\sqrt[3]{t}} \cdot \frac{\sqrt[3]{t}}{c}\right)\right)\\
\mathbf{elif}\;\left(x \cdot 9\right) \cdot y \le -4.44461852404034349 \cdot 10^{-84}:\\
\;\;\;\;\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}\;\left(x \cdot 9\right) \cdot y \le 6.32049670678102 \cdot 10^{-96}:\\
\;\;\;\;\left(\frac{\frac{b}{c}}{z} + 9 \cdot \left(\frac{y}{c} \cdot \frac{x}{z}\right)\right) - 4 \cdot \left(\left(a \cdot \left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)\right) \cdot \frac{\sqrt[3]{t}}{c}\right)\\
\mathbf{elif}\;\left(x \cdot 9\right) \cdot y \le 1.2680964165413392 \cdot 10^{182}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \frac{a}{\frac{c}{t}}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x}{\frac{c}{y} \cdot z}\right) - 4 \cdot \left(\left(a \cdot \left(\sqrt[3]{t} \cdot \sqrt[3]{t}\right)\right) \cdot \frac{\sqrt[3]{t}}{c}\right)\\
\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 ((((double) (((double) (x * 9.0)) * y)) <= -5.2031824371176434e+281)) {
VAR = ((double) (((double) (((double) (b / ((double) (z * c)))) + ((double) (9.0 * ((double) (((double) (y / c)) * ((double) (x / z)))))))) - ((double) (4.0 * ((double) (((double) (a * ((double) (((double) cbrt(t)) * ((double) (((double) cbrt(((double) cbrt(t)))) * ((double) cbrt(((double) cbrt(t)))))))))) * ((double) (((double) cbrt(((double) cbrt(t)))) * ((double) (((double) cbrt(t)) / c))))))))));
} else {
double VAR_1;
if ((((double) (((double) (x * 9.0)) * y)) <= -4.4446185240403435e-84)) {
VAR_1 = ((double) (((double) (((double) (b / ((double) (z * c)))) + ((double) (9.0 * ((double) (((double) (x * y)) / ((double) (z * c)))))))) - ((double) (4.0 * ((double) (a / ((double) (c / t))))))));
} else {
double VAR_2;
if ((((double) (((double) (x * 9.0)) * y)) <= 6.32049670678102e-96)) {
VAR_2 = ((double) (((double) (((double) (((double) (b / c)) / z)) + ((double) (9.0 * ((double) (((double) (y / c)) * ((double) (x / z)))))))) - ((double) (4.0 * ((double) (((double) (a * ((double) (((double) cbrt(t)) * ((double) cbrt(t)))))) * ((double) (((double) cbrt(t)) / c))))))));
} else {
double VAR_3;
if ((((double) (((double) (x * 9.0)) * y)) <= 1.2680964165413392e+182)) {
VAR_3 = ((double) (((double) (((double) (b / ((double) (z * c)))) + ((double) (9.0 * ((double) (((double) (x * y)) / ((double) (z * c)))))))) - ((double) (4.0 * ((double) (a / ((double) (c / t))))))));
} else {
VAR_3 = ((double) (((double) (((double) (b / ((double) (z * c)))) + ((double) (9.0 * ((double) (x / ((double) (((double) (c / y)) * z)))))))) - ((double) (4.0 * ((double) (((double) (a * ((double) (((double) cbrt(t)) * ((double) cbrt(t)))))) * ((double) (((double) cbrt(t)) / c))))))));
}
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.5 |
|---|---|
| Target | 14.4 |
| Herbie | 8.0 |
if (* (* x 9.0) y) < -5.2031824371176434e+281Initial program 55.3
Taylor expanded around 0 52.4
rmApplied *-un-lft-identity52.4
Applied add-cube-cbrt52.4
Applied associate-*r*52.4
Applied times-frac52.2
Simplified52.2
rmApplied *-commutative52.2
Applied *-commutative52.2
Applied times-frac11.9
rmApplied add-cube-cbrt11.9
Applied associate-*r*11.9
Applied associate-*r*11.9
Applied associate-*l*10.8
if -5.2031824371176434e+281 < (* (* x 9.0) y) < -4.4446185240403435e-84 or 6.32049670678102e-96 < (* (* x 9.0) y) < 1.2680964165413392e+182Initial program 17.7
Taylor expanded around 0 8.1
rmApplied associate-/l*6.8
if -4.4446185240403435e-84 < (* (* x 9.0) y) < 6.32049670678102e-96Initial program 16.7
Taylor expanded around 0 7.8
rmApplied *-un-lft-identity7.8
Applied add-cube-cbrt8.3
Applied associate-*r*8.3
Applied times-frac7.3
Simplified7.3
rmApplied *-commutative7.3
Applied *-commutative7.3
Applied times-frac8.8
rmApplied *-commutative8.8
Applied associate-/r*8.4
if 1.2680964165413392e+182 < (* (* x 9.0) y) Initial program 36.8
Taylor expanded around 0 32.0
rmApplied *-un-lft-identity32.0
Applied add-cube-cbrt32.1
Applied associate-*r*32.1
Applied times-frac31.3
Simplified31.3
rmApplied *-commutative31.3
Applied *-commutative31.3
Applied times-frac12.9
rmApplied clear-num13.2
Applied frac-times11.2
Simplified11.2
Final simplification8.0
herbie shell --seed 2020114
(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)))