double f(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double r13651173 = x;
double r13651174 = y;
double r13651175 = z;
double r13651176 = r13651174 * r13651175;
double r13651177 = t;
double r13651178 = a;
double r13651179 = r13651177 * r13651178;
double r13651180 = r13651176 - r13651179;
double r13651181 = r13651173 * r13651180;
double r13651182 = b;
double r13651183 = c;
double r13651184 = r13651183 * r13651175;
double r13651185 = i;
double r13651186 = r13651185 * r13651178;
double r13651187 = r13651184 - r13651186;
double r13651188 = r13651182 * r13651187;
double r13651189 = r13651181 - r13651188;
double r13651190 = j;
double r13651191 = r13651183 * r13651177;
double r13651192 = r13651185 * r13651174;
double r13651193 = r13651191 - r13651192;
double r13651194 = r13651190 * r13651193;
double r13651195 = r13651189 + r13651194;
return r13651195;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
double r13651196 = b;
double r13651197 = 4.281799399469648e-78;
bool r13651198 = r13651196 <= r13651197;
double r13651199 = j;
double r13651200 = t;
double r13651201 = c;
double r13651202 = r13651200 * r13651201;
double r13651203 = y;
double r13651204 = i;
double r13651205 = r13651203 * r13651204;
double r13651206 = r13651202 - r13651205;
double r13651207 = r13651199 * r13651206;
double r13651208 = z;
double r13651209 = r13651203 * r13651208;
double r13651210 = a;
double r13651211 = r13651200 * r13651210;
double r13651212 = r13651209 - r13651211;
double r13651213 = x;
double r13651214 = r13651212 * r13651213;
double r13651215 = cbrt(r13651196);
double r13651216 = r13651208 * r13651201;
double r13651217 = r13651215 * r13651216;
double r13651218 = r13651215 * r13651215;
double r13651219 = r13651217 * r13651218;
double r13651220 = r13651196 * r13651204;
double r13651221 = -r13651220;
double r13651222 = r13651210 * r13651221;
double r13651223 = r13651219 + r13651222;
double r13651224 = r13651214 - r13651223;
double r13651225 = r13651207 + r13651224;
double r13651226 = sqrt(r13651196);
double r13651227 = r13651210 * r13651204;
double r13651228 = r13651216 - r13651227;
double r13651229 = r13651228 * r13651226;
double r13651230 = r13651226 * r13651229;
double r13651231 = r13651214 - r13651230;
double r13651232 = r13651207 + r13651231;
double r13651233 = r13651198 ? r13651225 : r13651232;
return r13651233;
}
\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 4.281799399469648 \cdot 10^{-78}:\\
\;\;\;\;j \cdot \left(t \cdot c - y \cdot i\right) + \left(\left(y \cdot z - t \cdot a\right) \cdot x - \left(\left(\sqrt[3]{b} \cdot \left(z \cdot c\right)\right) \cdot \left(\sqrt[3]{b} \cdot \sqrt[3]{b}\right) + a \cdot \left(-b \cdot i\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;j \cdot \left(t \cdot c - y \cdot i\right) + \left(\left(y \cdot z - t \cdot a\right) \cdot x - \sqrt{b} \cdot \left(\left(z \cdot c - a \cdot i\right) \cdot \sqrt{b}\right)\right)\\
\end{array}


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
if b < 4.281799399469648e-78Initial program 12.6
rmApplied add-cube-cbrt12.9
Applied associate-*l*12.9
rmApplied sub-neg12.9
Applied distribute-lft-in12.9
Applied distribute-rgt-in12.9
Simplified12.1
if 4.281799399469648e-78 < b Initial program 8.1
rmApplied add-sqr-sqrt8.2
Applied associate-*l*8.3
Final simplification11.1
herbie shell --seed 2019102
(FPCore (x y z t a b c i j)
:name "Linear.Matrix:det33 from linear-1.19.1.3"
(+ (- (* x (- (* y z) (* t a))) (* b (- (* c z) (* i a)))) (* j (- (* c t) (* i y)))))