\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}\;x \le -2.64994232596744865 \cdot 10^{71}:\\
\;\;\;\;\mathsf{fma}\left(-4, \left(t \cdot a\right) \cdot \frac{1}{c}, \mathsf{fma}\left(9, \frac{x}{\frac{z \cdot c}{y}}, \frac{b}{z \cdot c}\right)\right)\\
\mathbf{elif}\;x \le -3.0898772136286032 \cdot 10^{-226}:\\
\;\;\;\;\mathsf{fma}\left(-4, t \cdot \frac{a}{c}, \frac{1}{z} \cdot \frac{\mathsf{fma}\left(9 \cdot x, y, b\right)}{c}\right)\\
\mathbf{elif}\;x \le 7.36682810843794424 \cdot 10^{-270}:\\
\;\;\;\;\mathsf{fma}\left(-4, \left(t \cdot a\right) \cdot \frac{1}{c}, \mathsf{fma}\left(9, \frac{x}{\frac{z \cdot c}{y}}, \frac{b}{z \cdot c}\right)\right)\\
\mathbf{elif}\;x \le 1.52002146635464452 \cdot 10^{84}:\\
\;\;\;\;\mathsf{fma}\left(-4, \frac{t}{\frac{c}{a}}, \frac{\mathsf{fma}\left(x, 9 \cdot y, b\right)}{z \cdot c}\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-4, \left(t \cdot a\right) \cdot \frac{1}{c}, \mathsf{fma}\left(9, \frac{x}{\frac{z \cdot c}{y}}, \frac{b}{z \cdot c}\right)\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 temp;
if ((x <= -2.6499423259674487e+71)) {
temp = fma(-4.0, ((t * a) * (1.0 / c)), fma(9.0, (x / ((z * c) / y)), (b / (z * c))));
} else {
double temp_1;
if ((x <= -3.0898772136286032e-226)) {
temp_1 = fma(-4.0, (t * (a / c)), ((1.0 / z) * (fma((9.0 * x), y, b) / c)));
} else {
double temp_2;
if ((x <= 7.366828108437944e-270)) {
temp_2 = fma(-4.0, ((t * a) * (1.0 / c)), fma(9.0, (x / ((z * c) / y)), (b / (z * c))));
} else {
double temp_3;
if ((x <= 1.5200214663546445e+84)) {
temp_3 = fma(-4.0, (t / (c / a)), (fma(x, (9.0 * y), b) / (z * c)));
} else {
temp_3 = fma(-4.0, ((t * a) * (1.0 / c)), fma(9.0, (x / ((z * c) / y)), (b / (z * c))));
}
temp_2 = temp_3;
}
temp_1 = temp_2;
}
temp = temp_1;
}
return temp;
}




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.8 |
|---|---|
| Target | 14.7 |
| Herbie | 10.1 |
if x < -2.6499423259674487e+71 or -3.0898772136286032e-226 < x < 7.366828108437944e-270 or 1.5200214663546445e+84 < x Initial program 24.6
Simplified16.6
rmApplied *-un-lft-identity16.6
Applied times-frac16.1
Simplified16.1
Taylor expanded around 0 16.1
Simplified16.1
rmApplied associate-/l*12.4
rmApplied div-inv12.4
Applied associate-*r*13.4
if -2.6499423259674487e+71 < x < -3.0898772136286032e-226Initial program 18.7
Simplified9.2
rmApplied *-un-lft-identity9.2
Applied times-frac8.0
Simplified8.0
rmApplied *-un-lft-identity8.0
Applied times-frac7.7
Simplified7.6
if 7.366828108437944e-270 < x < 1.5200214663546445e+84Initial program 17.1
Simplified8.1
rmApplied associate-/l*7.6
Final simplification10.1
herbie shell --seed 2020066 +o rules:numerics
(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)))