\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}\;t \le -6.08335225107627983 \cdot 10^{-33}:\\
\;\;\;\;t \cdot \left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + \left(\sqrt[3]{\left(j \cdot 27\right) \cdot k} \cdot \sqrt[3]{\left(j \cdot 27\right) \cdot k}\right) \cdot \sqrt[3]{\left(j \cdot 27\right) \cdot k}\right)\right)\\
\mathbf{elif}\;t \le 2.11917748355629525 \cdot 10^{-168}:\\
\;\;\;\;t \cdot \left(0 - a \cdot 4\right) + \left(b \cdot c - \left(\left(x \cdot 4\right) \cdot i + j \cdot \left(27 \cdot k\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(t \cdot \left(\sqrt[3]{\left(\left(x \cdot 18\right) \cdot y\right) \cdot z - a \cdot 4} \cdot \sqrt[3]{\left(\left(x \cdot 18\right) \cdot y\right) \cdot z - a \cdot 4}\right)\right) \cdot \sqrt[3]{\left(\left(x \cdot 18\right) \cdot y\right) \cdot z - a \cdot 4} + \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 ((((((((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 VAR;
if ((t <= -6.08335225107628e-33)) {
VAR = ((t * ((((x * 18.0) * y) * z) - (a * 4.0))) + ((b * c) - (((x * 4.0) * i) + ((cbrt(((j * 27.0) * k)) * cbrt(((j * 27.0) * k))) * cbrt(((j * 27.0) * k))))));
} else {
double VAR_1;
if ((t <= 2.1191774835562953e-168)) {
VAR_1 = ((t * (0.0 - (a * 4.0))) + ((b * c) - (((x * 4.0) * i) + (j * (27.0 * k)))));
} else {
VAR_1 = (((t * (cbrt(((((x * 18.0) * y) * z) - (a * 4.0))) * cbrt(((((x * 18.0) * y) * z) - (a * 4.0))))) * cbrt(((((x * 18.0) * y) * z) - (a * 4.0)))) + ((b * c) - (((x * 4.0) * i) + ((j * 27.0) * k))));
}
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



Bits error versus i



Bits error versus j



Bits error versus k
Results
if t < -6.08335225107628e-33Initial program 2.1
Simplified2.1
rmApplied add-cube-cbrt2.3
if -6.08335225107628e-33 < t < 2.1191774835562953e-168Initial program 8.4
Simplified8.4
rmApplied associate-*l*8.4
Taylor expanded around 0 6.0
if 2.1191774835562953e-168 < t Initial program 3.7
Simplified3.7
rmApplied add-cube-cbrt4.1
Applied associate-*r*4.1
Final simplification4.5
herbie shell --seed 2020091
(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) y) z) t) (* (* a 4) t)) (* b c)) (* (* x 4) i)) (* (* j 27) k)))