\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}\;z \le -3.1961472743277424 \cdot 10^{-28}:\\
\;\;\;\;\left(\frac{\frac{b}{z}}{c} + \frac{9}{\frac{\frac{z \cdot c}{y}}{x}}\right) - 4 \cdot \frac{a \cdot t}{c}\\
\mathbf{elif}\;z \le 2.03165913720499016 \cdot 10^{-81}:\\
\;\;\;\;\frac{\frac{\left(\left(x \cdot 9\right) \cdot y - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b}{c}}{z}\\
\mathbf{elif}\;z \le 1.02952530981005718 \cdot 10^{127}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + \frac{9}{\frac{\frac{z \cdot c}{y}}{x}}\right) - 4 \cdot \frac{a}{\frac{c}{t}}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + \left(9 \cdot \frac{y}{c}\right) \cdot \frac{x}{z}\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 ((z <= -3.1961472743277424e-28)) {
VAR = ((((b / z) / c) + (9.0 / (((z * c) / y) / x))) - (4.0 * ((a * t) / c)));
} else {
double VAR_1;
if ((z <= 2.0316591372049902e-81)) {
VAR_1 = ((((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / c) / z);
} else {
double VAR_2;
if ((z <= 1.0295253098100572e+127)) {
VAR_2 = (((b / (z * c)) + (9.0 / (((z * c) / y) / x))) - (4.0 * (a / (c / t))));
} else {
VAR_2 = (((b / (z * c)) + ((9.0 * (y / c)) * (x / z))) - (4.0 * ((a * t) / c)));
}
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.7 |
|---|---|
| Target | 14.5 |
| Herbie | 8.3 |
if z < -3.1961472743277424e-28Initial program 29.0
Taylor expanded around 0 13.4
rmApplied associate-/l*11.3
Applied associate-*r/11.4
rmApplied associate-/l*11.3
rmApplied associate-/r*8.8
if -3.1961472743277424e-28 < z < 2.0316591372049902e-81Initial program 5.9
rmApplied *-commutative5.9
Applied associate-/r*5.9
if 2.0316591372049902e-81 < z < 1.0295253098100572e+127Initial program 14.3
Taylor expanded around 0 10.0
rmApplied associate-/l*8.9
Applied associate-*r/8.9
rmApplied associate-/l*8.9
rmApplied associate-/l*8.6
if 1.0295253098100572e+127 < z Initial program 37.2
Taylor expanded around 0 15.3
rmApplied *-commutative15.3
Applied *-commutative15.3
Applied times-frac11.3
Applied associate-*r*11.3
Final simplification8.3
herbie shell --seed 2020078
(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)))