\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}\;x \le -1.2254010690602987 \cdot 10^{95} \lor \neg \left(x \le 1.9490316857655669 \cdot 10^{160}\right):\\
\;\;\;\;t \cdot \left(18 \cdot \left(x \cdot \left(z \cdot y\right)\right) - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + \left(\sqrt[3]{j} \cdot \left(\left(\sqrt[3]{\sqrt[3]{j}} \cdot \sqrt[3]{\sqrt[3]{j}}\right) \cdot \sqrt[3]{\sqrt[3]{j}}\right)\right) \cdot \left(\sqrt[3]{j} \cdot \left(27 \cdot k\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t \cdot \left(\left(x \cdot \left(18 \cdot y\right)\right) \cdot z - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \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 ((double) (((double) (((double) (((double) (((double) (((double) (((double) (((double) (x * 18.0)) * y)) * z)) * t)) - ((double) (((double) (a * 4.0)) * t)))) + ((double) (b * c)))) - ((double) (((double) (x * 4.0)) * i)))) - ((double) (((double) (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 VAR;
if (((x <= -1.2254010690602987e+95) || !(x <= 1.949031685765567e+160))) {
VAR = ((double) (((double) (t * ((double) (((double) (18.0 * ((double) (x * ((double) (z * y)))))) - ((double) (a * 4.0)))))) + ((double) (((double) (b * c)) - ((double) (((double) (((double) (x * 4.0)) * i)) + ((double) (((double) (((double) cbrt(j)) * ((double) (((double) (((double) cbrt(((double) cbrt(j)))) * ((double) cbrt(((double) cbrt(j)))))) * ((double) cbrt(((double) cbrt(j)))))))) * ((double) (((double) cbrt(j)) * ((double) (27.0 * k))))))))))));
} else {
VAR = ((double) (((double) (t * ((double) (((double) (((double) (x * ((double) (18.0 * y)))) * z)) - ((double) (a * 4.0)))))) + ((double) (((double) (b * c)) - ((double) (((double) (((double) (x * 4.0)) * i)) + ((double) (((double) (j * 27.0)) * k))))))));
}
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



Bits error versus i



Bits error versus j



Bits error versus k
Results
if x < -1.2254010690602987e+95 or 1.949031685765567e+160 < x Initial program 17.3
Simplified17.3
Taylor expanded around inf 9.3
rmApplied associate-*l*9.3
rmApplied add-cube-cbrt9.5
Applied associate-*l*9.5
rmApplied add-cube-cbrt9.5
if -1.2254010690602987e+95 < x < 1.949031685765567e+160Initial program 3.0
Simplified3.0
rmApplied associate-*l*3.1
Final simplification4.2
herbie shell --seed 2020140
(FPCore (x y z t a b c i j k)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1"
:precision binary64
(- (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) (* (* j 27.0) k)))