\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}\;b \le -3.31775490019201837 \cdot 10^{84} \lor \neg \left(b \le 4954531051790429180\right):\\
\;\;\;\;\left(\left(x \cdot \left(\sqrt[3]{y \cdot z - t \cdot a} \cdot \sqrt[3]{y \cdot z - t \cdot a}\right)\right) \cdot \sqrt[3]{y \cdot z - t \cdot a} - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \left(c \cdot t - i \cdot y\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x \cdot \left(y \cdot z - t \cdot a\right) - \left(\left(b \cdot c\right) \cdot z + \left(b \cdot i\right) \cdot \left(-a\right)\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 temp;
if (((b <= -3.3177549001920184e+84) || !(b <= 4.954531051790429e+18))) {
temp = ((((x * (cbrt(((y * z) - (t * a))) * cbrt(((y * z) - (t * a))))) * cbrt(((y * z) - (t * a)))) - (b * ((c * z) - (i * a)))) + (j * ((c * t) - (i * y))));
} else {
temp = (((x * ((y * z) - (t * a))) - (((b * c) * z) + ((b * i) * -a))) + (j * ((c * t) - (i * y))));
}
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
Results
if b < -3.3177549001920184e+84 or 4.954531051790429e+18 < b Initial program 7.3
rmApplied add-cube-cbrt7.5
Applied associate-*r*7.5
if -3.3177549001920184e+84 < b < 4.954531051790429e+18Initial program 14.7
rmApplied sub-neg14.7
Applied distribute-lft-in14.7
rmApplied distribute-rgt-neg-in14.7
Applied associate-*r*12.5
rmApplied associate-*r*10.5
Final simplification9.6
herbie shell --seed 2020060
(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)))))