\left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \left(c \cdot t - i \cdot y\right)
\begin{array}{l}
\mathbf{if}\;x \le -4.9381594717483803 \cdot 10^{27}:\\
\;\;\;\;\left(x \cdot \left(y \cdot z - t \cdot a\right) - \left(\sqrt[3]{b} \cdot \sqrt[3]{b}\right) \cdot \left(\sqrt[3]{b} \cdot \left(c \cdot z - i \cdot a\right)\right)\right) + j \cdot \left(c \cdot t - i \cdot y\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\left(\sqrt[3]{\left(x \cdot y\right) \cdot z} \cdot \sqrt[3]{\left(x \cdot y\right) \cdot z}\right) \cdot \sqrt[3]{\left(x \cdot y\right) \cdot z} + \left(x \cdot \left(-t\right)\right) \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \left(c \cdot t - i \cdot y\right)\\
\end{array}double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
return (((x * ((y * z) - (t * a))) - (b * ((c * z) - (i * a)))) + (j * ((c * t) - (i * y))));
}
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double VAR;
if ((x <= -4.9381594717483803e+27)) {
VAR = (((x * ((y * z) - (t * a))) - ((cbrt(b) * cbrt(b)) * (cbrt(b) * ((c * z) - (i * a))))) + (j * ((c * t) - (i * y))));
} else {
VAR = (((((cbrt(((x * y) * z)) * cbrt(((x * y) * z))) * cbrt(((x * y) * z))) + ((x * -t) * a)) - (b * ((c * z) - (i * a)))) + (j * ((c * t) - (i * y))));
}
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
Results
if x < -4.9381594717483803e+27Initial program 7.8
rmApplied add-cube-cbrt8.0
Applied associate-*l*8.0
if -4.9381594717483803e+27 < x Initial program 13.3
rmApplied sub-neg13.3
Applied distribute-lft-in13.3
rmApplied distribute-lft-neg-in13.3
Applied associate-*r*12.5
rmApplied associate-*r*11.5
rmApplied add-cube-cbrt11.6
Final simplification11.0
herbie shell --seed 2020105
(FPCore (x y z t a b c i j)
:name "Linear.Matrix:det33 from linear-1.19.1.3"
:precision binary64
(+ (- (* x (- (* y z) (* t a))) (* b (- (* c z) (* i a)))) (* j (- (* c t) (* i y)))))