\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}\;j \le -2.07354128237016473 \cdot 10^{46} \lor \neg \left(j \le 14736619545097708\right):\\
\;\;\;\;\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(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(t \cdot \left(j \cdot c\right) + -1 \cdot \left(i \cdot \left(y \cdot j\right)\right)\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 (((j <= -2.0735412823701647e+46) || !(j <= 14736619545097708.0))) {
VAR = (((x * ((y * z) - (t * a))) - ((cbrt(b) * cbrt(b)) * (cbrt(b) * ((c * z) - (i * a))))) + (j * ((c * t) - (i * y))));
} else {
VAR = (((x * ((y * z) - (t * a))) - (b * ((c * z) - (i * a)))) + ((t * (j * c)) + (-1.0 * (i * (y * j)))));
}
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 j < -2.0735412823701647e+46 or 14736619545097708.0 < j Initial program 7.6
rmApplied add-cube-cbrt7.8
Applied associate-*l*7.8
if -2.0735412823701647e+46 < j < 14736619545097708.0Initial program 15.1
rmApplied add-cube-cbrt15.3
Applied associate-*l*15.3
rmApplied sub-neg15.3
Applied distribute-lft-in15.3
Applied distribute-lft-in15.3
Simplified13.1
Simplified13.0
Taylor expanded around inf 10.5
Final simplification9.6
herbie shell --seed 2020102
(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)))))