\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}\;t \le -1.33910392604837064 \cdot 10^{-180}:\\
\;\;\;\;\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}\;t \le 6.31146352264548959 \cdot 10^{-228}:\\
\;\;\;\;\frac{\frac{\left(\left(x \cdot 9\right) \cdot y - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b}{z}}{c}\\
\mathbf{elif}\;t \le 31928303.1963208877:\\
\;\;\;\;\left(\frac{\sqrt[3]{b} \cdot \sqrt[3]{b}}{z} \cdot \frac{\sqrt[3]{b}}{c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \frac{a \cdot t}{c}\\
\mathbf{else}:\\
\;\;\;\;\left(\frac{b}{z \cdot c} + 9 \cdot \frac{x \cdot y}{z \cdot c}\right) - 4 \cdot \left(\frac{a}{c} \cdot t\right)\\
\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 ((t <= -1.3391039260483706e-180)) {
VAR = (((b / (z * c)) + (9.0 * ((x * y) / (z * c)))) - (4.0 * (a / (c / t))));
} else {
double VAR_1;
if ((t <= 6.31146352264549e-228)) {
VAR_1 = ((((((x * 9.0) * y) - (((z * 4.0) * t) * a)) + b) / z) / c);
} else {
double VAR_2;
if ((t <= 31928303.196320888)) {
VAR_2 = (((((cbrt(b) * cbrt(b)) / z) * (cbrt(b) / c)) + (9.0 * ((x * y) / (z * c)))) - (4.0 * ((a * t) / c)));
} else {
VAR_2 = (((b / (z * c)) + (9.0 * ((x * y) / (z * c)))) - (4.0 * ((a / c) * t)));
}
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.5 |
| Herbie | 10.3 |
if t < -1.3391039260483706e-180Initial program 23.2
Taylor expanded around 0 12.4
rmApplied associate-/l*10.6
if -1.3391039260483706e-180 < t < 6.31146352264549e-228Initial program 12.0
rmApplied associate-/r*10.0
if 6.31146352264549e-228 < t < 31928303.196320888Initial program 14.4
Taylor expanded around 0 8.3
rmApplied add-cube-cbrt8.6
Applied times-frac9.3
if 31928303.196320888 < t Initial program 29.5
Taylor expanded around 0 15.6
rmApplied associate-/l*11.3
rmApplied associate-/r/11.3
Final simplification10.3
herbie shell --seed 2020092
(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)))