\left(\left(\left(\left(\left(x \cdot y - z \cdot t\right) \cdot \left(a \cdot b - c \cdot i\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(y0 \cdot b - y1 \cdot i\right)\right) + \left(x \cdot y2 - z \cdot y3\right) \cdot \left(y0 \cdot c - y1 \cdot a\right)\right) + \left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - y5 \cdot i\right)\right) - \left(t \cdot y2 - y \cdot y3\right) \cdot \left(y4 \cdot c - y5 \cdot a\right)\right) + \left(k \cdot y2 - j \cdot y3\right) \cdot \left(y4 \cdot y1 - y5 \cdot y0\right)
\begin{array}{l}
\mathbf{if}\;i \le -8.842172059440902 \cdot 10^{-229}:\\
\;\;\;\;\left(y1 \cdot y4 - y0 \cdot y5\right) \cdot \left(k \cdot y2 - j \cdot y3\right) + \left(\left(\left(\left(y0 \cdot c - y1 \cdot a\right) \cdot \left(x \cdot y2 - y3 \cdot z\right) + \left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - c \cdot i\right) - \sqrt[3]{\left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)} \cdot \left(\sqrt[3]{\left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)} \cdot \sqrt[3]{\sqrt[3]{b \cdot y0 - i \cdot y1} \cdot \left(\left(x \cdot j - z \cdot k\right) \cdot \left(\sqrt[3]{b \cdot y0 - i \cdot y1} \cdot \sqrt[3]{b \cdot y0 - i \cdot y1}\right)\right)}\right)\right)\right) + \left(y4 \cdot b - i \cdot y5\right) \cdot \left(j \cdot t - y \cdot k\right)\right) - \left(y2 \cdot t - y \cdot y3\right) \cdot \left(y4 \cdot c - a \cdot y5\right)\right)\\
\mathbf{elif}\;i \le 6.082724873086382 \cdot 10^{-126}:\\
\;\;\;\;\left(y1 \cdot y4 - y0 \cdot y5\right) \cdot \left(k \cdot y2 - j \cdot y3\right) + \left(\left(\left(y4 \cdot b - i \cdot y5\right) \cdot \left(j \cdot t - y \cdot k\right) + \left(\left(\left(\left(\left(c \cdot i - a \cdot b\right) \cdot t\right) \cdot z - c \cdot \left(\left(y \cdot i\right) \cdot x\right)\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right) + \left(y0 \cdot c - y1 \cdot a\right) \cdot \left(x \cdot y2 - y3 \cdot z\right)\right)\right) - \left(y2 \cdot t - y \cdot y3\right) \cdot \left(y4 \cdot c - a \cdot y5\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(y1 \cdot y4 - y0 \cdot y5\right) \cdot \left(k \cdot y2 - j \cdot y3\right) + \left(\left(\left(y4 \cdot b - i \cdot y5\right) \cdot \left(j \cdot t - y \cdot k\right) + \left(\left(\sqrt[3]{\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - c \cdot i\right)} \cdot \left(\sqrt[3]{\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - c \cdot i\right)} \cdot \sqrt[3]{\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - c \cdot i\right)}\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right) + \left(y0 \cdot c - y1 \cdot a\right) \cdot \left(x \cdot y2 - y3 \cdot z\right)\right)\right) - \left(y2 \cdot t - y \cdot y3\right) \cdot \left(y4 \cdot c - a \cdot y5\right)\right)\\
\end{array}double f(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k, double y0, double y1, double y2, double y3, double y4, double y5) {
double r19491254 = x;
double r19491255 = y;
double r19491256 = r19491254 * r19491255;
double r19491257 = z;
double r19491258 = t;
double r19491259 = r19491257 * r19491258;
double r19491260 = r19491256 - r19491259;
double r19491261 = a;
double r19491262 = b;
double r19491263 = r19491261 * r19491262;
double r19491264 = c;
double r19491265 = i;
double r19491266 = r19491264 * r19491265;
double r19491267 = r19491263 - r19491266;
double r19491268 = r19491260 * r19491267;
double r19491269 = j;
double r19491270 = r19491254 * r19491269;
double r19491271 = k;
double r19491272 = r19491257 * r19491271;
double r19491273 = r19491270 - r19491272;
double r19491274 = y0;
double r19491275 = r19491274 * r19491262;
double r19491276 = y1;
double r19491277 = r19491276 * r19491265;
double r19491278 = r19491275 - r19491277;
double r19491279 = r19491273 * r19491278;
double r19491280 = r19491268 - r19491279;
double r19491281 = y2;
double r19491282 = r19491254 * r19491281;
double r19491283 = y3;
double r19491284 = r19491257 * r19491283;
double r19491285 = r19491282 - r19491284;
double r19491286 = r19491274 * r19491264;
double r19491287 = r19491276 * r19491261;
double r19491288 = r19491286 - r19491287;
double r19491289 = r19491285 * r19491288;
double r19491290 = r19491280 + r19491289;
double r19491291 = r19491258 * r19491269;
double r19491292 = r19491255 * r19491271;
double r19491293 = r19491291 - r19491292;
double r19491294 = y4;
double r19491295 = r19491294 * r19491262;
double r19491296 = y5;
double r19491297 = r19491296 * r19491265;
double r19491298 = r19491295 - r19491297;
double r19491299 = r19491293 * r19491298;
double r19491300 = r19491290 + r19491299;
double r19491301 = r19491258 * r19491281;
double r19491302 = r19491255 * r19491283;
double r19491303 = r19491301 - r19491302;
double r19491304 = r19491294 * r19491264;
double r19491305 = r19491296 * r19491261;
double r19491306 = r19491304 - r19491305;
double r19491307 = r19491303 * r19491306;
double r19491308 = r19491300 - r19491307;
double r19491309 = r19491271 * r19491281;
double r19491310 = r19491269 * r19491283;
double r19491311 = r19491309 - r19491310;
double r19491312 = r19491294 * r19491276;
double r19491313 = r19491296 * r19491274;
double r19491314 = r19491312 - r19491313;
double r19491315 = r19491311 * r19491314;
double r19491316 = r19491308 + r19491315;
return r19491316;
}
double f(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k, double y0, double y1, double y2, double y3, double y4, double y5) {
double r19491317 = i;
double r19491318 = -8.842172059440902e-229;
bool r19491319 = r19491317 <= r19491318;
double r19491320 = y1;
double r19491321 = y4;
double r19491322 = r19491320 * r19491321;
double r19491323 = y0;
double r19491324 = y5;
double r19491325 = r19491323 * r19491324;
double r19491326 = r19491322 - r19491325;
double r19491327 = k;
double r19491328 = y2;
double r19491329 = r19491327 * r19491328;
double r19491330 = j;
double r19491331 = y3;
double r19491332 = r19491330 * r19491331;
double r19491333 = r19491329 - r19491332;
double r19491334 = r19491326 * r19491333;
double r19491335 = c;
double r19491336 = r19491323 * r19491335;
double r19491337 = a;
double r19491338 = r19491320 * r19491337;
double r19491339 = r19491336 - r19491338;
double r19491340 = x;
double r19491341 = r19491340 * r19491328;
double r19491342 = z;
double r19491343 = r19491331 * r19491342;
double r19491344 = r19491341 - r19491343;
double r19491345 = r19491339 * r19491344;
double r19491346 = y;
double r19491347 = r19491340 * r19491346;
double r19491348 = t;
double r19491349 = r19491348 * r19491342;
double r19491350 = r19491347 - r19491349;
double r19491351 = b;
double r19491352 = r19491337 * r19491351;
double r19491353 = r19491335 * r19491317;
double r19491354 = r19491352 - r19491353;
double r19491355 = r19491350 * r19491354;
double r19491356 = r19491340 * r19491330;
double r19491357 = r19491342 * r19491327;
double r19491358 = r19491356 - r19491357;
double r19491359 = r19491351 * r19491323;
double r19491360 = r19491317 * r19491320;
double r19491361 = r19491359 - r19491360;
double r19491362 = r19491358 * r19491361;
double r19491363 = cbrt(r19491362);
double r19491364 = cbrt(r19491361);
double r19491365 = r19491364 * r19491364;
double r19491366 = r19491358 * r19491365;
double r19491367 = r19491364 * r19491366;
double r19491368 = cbrt(r19491367);
double r19491369 = r19491363 * r19491368;
double r19491370 = r19491363 * r19491369;
double r19491371 = r19491355 - r19491370;
double r19491372 = r19491345 + r19491371;
double r19491373 = r19491321 * r19491351;
double r19491374 = r19491317 * r19491324;
double r19491375 = r19491373 - r19491374;
double r19491376 = r19491330 * r19491348;
double r19491377 = r19491346 * r19491327;
double r19491378 = r19491376 - r19491377;
double r19491379 = r19491375 * r19491378;
double r19491380 = r19491372 + r19491379;
double r19491381 = r19491328 * r19491348;
double r19491382 = r19491346 * r19491331;
double r19491383 = r19491381 - r19491382;
double r19491384 = r19491321 * r19491335;
double r19491385 = r19491337 * r19491324;
double r19491386 = r19491384 - r19491385;
double r19491387 = r19491383 * r19491386;
double r19491388 = r19491380 - r19491387;
double r19491389 = r19491334 + r19491388;
double r19491390 = 6.082724873086382e-126;
bool r19491391 = r19491317 <= r19491390;
double r19491392 = r19491353 - r19491352;
double r19491393 = r19491392 * r19491348;
double r19491394 = r19491393 * r19491342;
double r19491395 = r19491346 * r19491317;
double r19491396 = r19491395 * r19491340;
double r19491397 = r19491335 * r19491396;
double r19491398 = r19491394 - r19491397;
double r19491399 = r19491398 - r19491362;
double r19491400 = r19491399 + r19491345;
double r19491401 = r19491379 + r19491400;
double r19491402 = r19491401 - r19491387;
double r19491403 = r19491334 + r19491402;
double r19491404 = cbrt(r19491355);
double r19491405 = r19491404 * r19491404;
double r19491406 = r19491404 * r19491405;
double r19491407 = r19491406 - r19491362;
double r19491408 = r19491407 + r19491345;
double r19491409 = r19491379 + r19491408;
double r19491410 = r19491409 - r19491387;
double r19491411 = r19491334 + r19491410;
double r19491412 = r19491391 ? r19491403 : r19491411;
double r19491413 = r19491319 ? r19491389 : r19491412;
return r19491413;
}



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



Bits error versus y0



Bits error versus y1



Bits error versus y2



Bits error versus y3



Bits error versus y4



Bits error versus y5
Results
if i < -8.842172059440902e-229Initial program 26.0
rmApplied add-cube-cbrt26.2
rmApplied add-cube-cbrt26.1
Applied associate-*r*26.1
if -8.842172059440902e-229 < i < 6.082724873086382e-126Initial program 26.5
Taylor expanded around -inf 30.1
Simplified26.5
if 6.082724873086382e-126 < i Initial program 25.5
rmApplied add-cube-cbrt25.5
Final simplification26.1
herbie shell --seed 2019104
(FPCore (x y z t a b c i j k y0 y1 y2 y3 y4 y5)
:name "Linear.Matrix:det44 from linear-1.19.1.3"
(+ (- (+ (+ (- (* (- (* x y) (* z t)) (- (* a b) (* c i))) (* (- (* x j) (* z k)) (- (* y0 b) (* y1 i)))) (* (- (* x y2) (* z y3)) (- (* y0 c) (* y1 a)))) (* (- (* t j) (* y k)) (- (* y4 b) (* y5 i)))) (* (- (* t y2) (* y y3)) (- (* y4 c) (* y5 a)))) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0)))))