\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k
\begin{array}{l}
\mathbf{if}\;z \le -4.15846637414394352 \cdot 10^{131} \lor \neg \left(z \le 3.14022389051985691 \cdot 10^{43}\right):\\
\;\;\;\;\mathsf{fma}\left(t, \left(\left(x \cdot 18\right) \cdot y\right) \cdot z - a \cdot 4, b \cdot c - \mathsf{fma}\left(x, 4 \cdot i, j \cdot \left(27 \cdot k\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t, {\left(18 \cdot \left(x \cdot \left(z \cdot y\right)\right)\right)}^{1} - a \cdot 4, b \cdot c - \mathsf{fma}\left(x, 4 \cdot i, \left(j \cdot 27\right) \cdot k\right)\right)\\
\end{array}double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return ((((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k));
}
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double temp;
if (((z <= -4.1584663741439435e+131) || !(z <= 3.140223890519857e+43))) {
temp = fma(t, ((((x * 18.0) * y) * z) - (a * 4.0)), ((b * c) - fma(x, (4.0 * i), (j * (27.0 * k)))));
} else {
temp = fma(t, (pow((18.0 * (x * (z * y))), 1.0) - (a * 4.0)), ((b * c) - fma(x, (4.0 * i), ((j * 27.0) * k))));
}
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




Bits error versus i




Bits error versus j




Bits error versus k
Results
| Original | 5.4 |
|---|---|
| Target | 1.7 |
| Herbie | 3.7 |
if z < -4.1584663741439435e+131 or 3.140223890519857e+43 < z Initial program 7.9
Simplified7.9
rmApplied associate-*l*8.0
if -4.1584663741439435e+131 < z < 3.140223890519857e+43Initial program 4.4
Simplified4.4
rmApplied pow14.4
Applied pow14.4
Applied pow14.4
Applied pow14.4
Applied pow-prod-down4.4
Applied pow-prod-down4.4
Applied pow-prod-down4.4
Simplified1.9
Final simplification3.7
herbie shell --seed 2020057 +o rules:numerics
(FPCore (x y z t a b c i j k)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, E"
:precision binary64
:herbie-target
(if (< t -1.6210815397541398e-69) (- (- (* (* 18 t) (* (* x y) z)) (* (+ (* a t) (* i x)) 4)) (- (* (* k j) 27) (* c b))) (if (< t 165.68027943805222) (+ (- (* (* 18 y) (* x (* z t))) (* (+ (* a t) (* i x)) 4)) (- (* c b) (* 27 (* k j)))) (- (- (* (* 18 t) (* (* x y) z)) (* (+ (* a t) (* i x)) 4)) (- (* (* k j) 27) (* c b)))))
(- (- (+ (- (* (* (* (* x 18) y) z) t) (* (* a 4) t)) (* b c)) (* (* x 4) i)) (* (* j 27) k)))