\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}\;y0 \le -13614926747.548742:\\
\;\;\;\;\left(\left(\left(\left(c \cdot y0 - y1 \cdot a\right) \cdot \left(x \cdot y2 - z \cdot y3\right) + \left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - i \cdot c\right) - \left(b \cdot y0 - y1 \cdot i\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right) + \left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - i \cdot y5\right)\right) - \left(a \cdot \left(y3 \cdot \left(y \cdot y5\right)\right) - \left(\left(\left(y5 \cdot t\right) \cdot y2\right) \cdot a + y3 \cdot \left(\left(y \cdot c\right) \cdot y4\right)\right)\right)\right) + \left(k \cdot y2 - y3 \cdot j\right) \cdot \left(y4 \cdot y1 - y0 \cdot y5\right)\\
\mathbf{elif}\;y0 \le 1.7112105416372722 \cdot 10^{-256}:\\
\;\;\;\;\left(\left(\left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - i \cdot y5\right) + \left(\left(\left(\left(y3 \cdot y1\right) \cdot z - x \cdot \left(y2 \cdot y1\right)\right) \cdot a - \left(c \cdot z\right) \cdot \left(y3 \cdot y0\right)\right) + \left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - i \cdot c\right) - \left(b \cdot y0 - y1 \cdot i\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right)\right) - \left(\sqrt[3]{t \cdot y2 - y \cdot y3} \cdot \sqrt[3]{t \cdot y2 - y \cdot y3}\right) \cdot \left(\sqrt[3]{t \cdot y2 - y \cdot y3} \cdot \left(y4 \cdot c - y5 \cdot a\right)\right)\right) + \left(k \cdot y2 - y3 \cdot j\right) \cdot \left(y4 \cdot y1 - y0 \cdot y5\right)\\
\mathbf{elif}\;y0 \le 9.717693116775471 \cdot 10^{-203}:\\
\;\;\;\;\left(y3 \cdot \left(j \cdot \left(y0 \cdot y5\right)\right) - \left(\left(\left(y4 \cdot j\right) \cdot y3\right) \cdot y1 + k \cdot \left(\left(y0 \cdot y5\right) \cdot y2\right)\right)\right) + \left(\left(\left(\left(c \cdot y0 - y1 \cdot a\right) \cdot \left(x \cdot y2 - z \cdot y3\right) + \left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - i \cdot c\right) - \left(b \cdot y0 - y1 \cdot i\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right) + \left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - i \cdot y5\right)\right) - \left(\sqrt[3]{t \cdot y2 - y \cdot y3} \cdot \sqrt[3]{t \cdot y2 - y \cdot y3}\right) \cdot \left(\sqrt[3]{t \cdot y2 - y \cdot y3} \cdot \left(y4 \cdot c - y5 \cdot a\right)\right)\right)\\
\mathbf{elif}\;y0 \le 4.308997712685301 \cdot 10^{-71}:\\
\;\;\;\;\left(k \cdot y2 - y3 \cdot j\right) \cdot \left(y4 \cdot y1 - y0 \cdot y5\right) + \left(\left(\left(\sqrt[3]{\left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - i \cdot y5\right)} \cdot \sqrt[3]{\left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - i \cdot y5\right)}\right) \cdot \sqrt[3]{\left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - i \cdot y5\right)} + \left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - i \cdot c\right) - \left(b \cdot y0 - y1 \cdot i\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right) - \left(y4 \cdot c - y5 \cdot a\right) \cdot \left(t \cdot y2 - y \cdot y3\right)\right)\\
\mathbf{elif}\;y0 \le 1.3561933432615603 \cdot 10^{+56}:\\
\;\;\;\;\left(k \cdot y2 - y3 \cdot j\right) \cdot \left(y4 \cdot y1 - y0 \cdot y5\right) + \left(\left(\left(\left(\left(\left(\left(c \cdot z\right) \cdot i\right) \cdot t - \left(i \cdot \left(\left(y \cdot c\right) \cdot x\right) + a \cdot \left(t \cdot \left(b \cdot z\right)\right)\right)\right) - \left(b \cdot y0 - y1 \cdot i\right) \cdot \left(x \cdot j - z \cdot k\right)\right) + \left(c \cdot y0 - y1 \cdot a\right) \cdot \left(x \cdot y2 - z \cdot y3\right)\right) + \left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - i \cdot y5\right)\right) - \left(y4 \cdot c - y5 \cdot a\right) \cdot \left(t \cdot y2 - y \cdot y3\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\left(\left(c \cdot y0 - y1 \cdot a\right) \cdot \left(x \cdot y2 - z \cdot y3\right) + \left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - i \cdot c\right) - \left(b \cdot y0 - y1 \cdot i\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right) + \left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - i \cdot y5\right)\right) - \left(a \cdot \left(y3 \cdot \left(y \cdot y5\right)\right) - \left(\left(\left(y5 \cdot t\right) \cdot y2\right) \cdot a + y3 \cdot \left(\left(y \cdot c\right) \cdot y4\right)\right)\right)\right) + \left(k \cdot y2 - y3 \cdot j\right) \cdot \left(y4 \cdot y1 - y0 \cdot y5\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 r7061280 = x;
double r7061281 = y;
double r7061282 = r7061280 * r7061281;
double r7061283 = z;
double r7061284 = t;
double r7061285 = r7061283 * r7061284;
double r7061286 = r7061282 - r7061285;
double r7061287 = a;
double r7061288 = b;
double r7061289 = r7061287 * r7061288;
double r7061290 = c;
double r7061291 = i;
double r7061292 = r7061290 * r7061291;
double r7061293 = r7061289 - r7061292;
double r7061294 = r7061286 * r7061293;
double r7061295 = j;
double r7061296 = r7061280 * r7061295;
double r7061297 = k;
double r7061298 = r7061283 * r7061297;
double r7061299 = r7061296 - r7061298;
double r7061300 = y0;
double r7061301 = r7061300 * r7061288;
double r7061302 = y1;
double r7061303 = r7061302 * r7061291;
double r7061304 = r7061301 - r7061303;
double r7061305 = r7061299 * r7061304;
double r7061306 = r7061294 - r7061305;
double r7061307 = y2;
double r7061308 = r7061280 * r7061307;
double r7061309 = y3;
double r7061310 = r7061283 * r7061309;
double r7061311 = r7061308 - r7061310;
double r7061312 = r7061300 * r7061290;
double r7061313 = r7061302 * r7061287;
double r7061314 = r7061312 - r7061313;
double r7061315 = r7061311 * r7061314;
double r7061316 = r7061306 + r7061315;
double r7061317 = r7061284 * r7061295;
double r7061318 = r7061281 * r7061297;
double r7061319 = r7061317 - r7061318;
double r7061320 = y4;
double r7061321 = r7061320 * r7061288;
double r7061322 = y5;
double r7061323 = r7061322 * r7061291;
double r7061324 = r7061321 - r7061323;
double r7061325 = r7061319 * r7061324;
double r7061326 = r7061316 + r7061325;
double r7061327 = r7061284 * r7061307;
double r7061328 = r7061281 * r7061309;
double r7061329 = r7061327 - r7061328;
double r7061330 = r7061320 * r7061290;
double r7061331 = r7061322 * r7061287;
double r7061332 = r7061330 - r7061331;
double r7061333 = r7061329 * r7061332;
double r7061334 = r7061326 - r7061333;
double r7061335 = r7061297 * r7061307;
double r7061336 = r7061295 * r7061309;
double r7061337 = r7061335 - r7061336;
double r7061338 = r7061320 * r7061302;
double r7061339 = r7061322 * r7061300;
double r7061340 = r7061338 - r7061339;
double r7061341 = r7061337 * r7061340;
double r7061342 = r7061334 + r7061341;
return r7061342;
}
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 r7061343 = y0;
double r7061344 = -13614926747.548742;
bool r7061345 = r7061343 <= r7061344;
double r7061346 = c;
double r7061347 = r7061346 * r7061343;
double r7061348 = y1;
double r7061349 = a;
double r7061350 = r7061348 * r7061349;
double r7061351 = r7061347 - r7061350;
double r7061352 = x;
double r7061353 = y2;
double r7061354 = r7061352 * r7061353;
double r7061355 = z;
double r7061356 = y3;
double r7061357 = r7061355 * r7061356;
double r7061358 = r7061354 - r7061357;
double r7061359 = r7061351 * r7061358;
double r7061360 = y;
double r7061361 = r7061352 * r7061360;
double r7061362 = t;
double r7061363 = r7061362 * r7061355;
double r7061364 = r7061361 - r7061363;
double r7061365 = b;
double r7061366 = r7061349 * r7061365;
double r7061367 = i;
double r7061368 = r7061367 * r7061346;
double r7061369 = r7061366 - r7061368;
double r7061370 = r7061364 * r7061369;
double r7061371 = r7061365 * r7061343;
double r7061372 = r7061348 * r7061367;
double r7061373 = r7061371 - r7061372;
double r7061374 = j;
double r7061375 = r7061352 * r7061374;
double r7061376 = k;
double r7061377 = r7061355 * r7061376;
double r7061378 = r7061375 - r7061377;
double r7061379 = r7061373 * r7061378;
double r7061380 = r7061370 - r7061379;
double r7061381 = r7061359 + r7061380;
double r7061382 = r7061362 * r7061374;
double r7061383 = r7061360 * r7061376;
double r7061384 = r7061382 - r7061383;
double r7061385 = y4;
double r7061386 = r7061385 * r7061365;
double r7061387 = y5;
double r7061388 = r7061367 * r7061387;
double r7061389 = r7061386 - r7061388;
double r7061390 = r7061384 * r7061389;
double r7061391 = r7061381 + r7061390;
double r7061392 = r7061360 * r7061387;
double r7061393 = r7061356 * r7061392;
double r7061394 = r7061349 * r7061393;
double r7061395 = r7061387 * r7061362;
double r7061396 = r7061395 * r7061353;
double r7061397 = r7061396 * r7061349;
double r7061398 = r7061360 * r7061346;
double r7061399 = r7061398 * r7061385;
double r7061400 = r7061356 * r7061399;
double r7061401 = r7061397 + r7061400;
double r7061402 = r7061394 - r7061401;
double r7061403 = r7061391 - r7061402;
double r7061404 = r7061376 * r7061353;
double r7061405 = r7061356 * r7061374;
double r7061406 = r7061404 - r7061405;
double r7061407 = r7061385 * r7061348;
double r7061408 = r7061343 * r7061387;
double r7061409 = r7061407 - r7061408;
double r7061410 = r7061406 * r7061409;
double r7061411 = r7061403 + r7061410;
double r7061412 = 1.7112105416372722e-256;
bool r7061413 = r7061343 <= r7061412;
double r7061414 = r7061356 * r7061348;
double r7061415 = r7061414 * r7061355;
double r7061416 = r7061353 * r7061348;
double r7061417 = r7061352 * r7061416;
double r7061418 = r7061415 - r7061417;
double r7061419 = r7061418 * r7061349;
double r7061420 = r7061346 * r7061355;
double r7061421 = r7061356 * r7061343;
double r7061422 = r7061420 * r7061421;
double r7061423 = r7061419 - r7061422;
double r7061424 = r7061423 + r7061380;
double r7061425 = r7061390 + r7061424;
double r7061426 = r7061362 * r7061353;
double r7061427 = r7061360 * r7061356;
double r7061428 = r7061426 - r7061427;
double r7061429 = cbrt(r7061428);
double r7061430 = r7061429 * r7061429;
double r7061431 = r7061385 * r7061346;
double r7061432 = r7061387 * r7061349;
double r7061433 = r7061431 - r7061432;
double r7061434 = r7061429 * r7061433;
double r7061435 = r7061430 * r7061434;
double r7061436 = r7061425 - r7061435;
double r7061437 = r7061436 + r7061410;
double r7061438 = 9.717693116775471e-203;
bool r7061439 = r7061343 <= r7061438;
double r7061440 = r7061374 * r7061408;
double r7061441 = r7061356 * r7061440;
double r7061442 = r7061385 * r7061374;
double r7061443 = r7061442 * r7061356;
double r7061444 = r7061443 * r7061348;
double r7061445 = r7061408 * r7061353;
double r7061446 = r7061376 * r7061445;
double r7061447 = r7061444 + r7061446;
double r7061448 = r7061441 - r7061447;
double r7061449 = r7061391 - r7061435;
double r7061450 = r7061448 + r7061449;
double r7061451 = 4.308997712685301e-71;
bool r7061452 = r7061343 <= r7061451;
double r7061453 = cbrt(r7061390);
double r7061454 = r7061453 * r7061453;
double r7061455 = r7061454 * r7061453;
double r7061456 = r7061455 + r7061380;
double r7061457 = r7061433 * r7061428;
double r7061458 = r7061456 - r7061457;
double r7061459 = r7061410 + r7061458;
double r7061460 = 1.3561933432615603e+56;
bool r7061461 = r7061343 <= r7061460;
double r7061462 = r7061420 * r7061367;
double r7061463 = r7061462 * r7061362;
double r7061464 = r7061398 * r7061352;
double r7061465 = r7061367 * r7061464;
double r7061466 = r7061365 * r7061355;
double r7061467 = r7061362 * r7061466;
double r7061468 = r7061349 * r7061467;
double r7061469 = r7061465 + r7061468;
double r7061470 = r7061463 - r7061469;
double r7061471 = r7061470 - r7061379;
double r7061472 = r7061471 + r7061359;
double r7061473 = r7061472 + r7061390;
double r7061474 = r7061473 - r7061457;
double r7061475 = r7061410 + r7061474;
double r7061476 = r7061461 ? r7061475 : r7061411;
double r7061477 = r7061452 ? r7061459 : r7061476;
double r7061478 = r7061439 ? r7061450 : r7061477;
double r7061479 = r7061413 ? r7061437 : r7061478;
double r7061480 = r7061345 ? r7061411 : r7061479;
return r7061480;
}



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 y0 < -13614926747.548742 or 1.3561933432615603e+56 < y0 Initial program 28.0
Taylor expanded around inf 28.7
if -13614926747.548742 < y0 < 1.7112105416372722e-256Initial program 25.1
rmApplied add-cube-cbrt25.2
Applied associate-*l*25.2
Taylor expanded around inf 26.3
Simplified27.6
if 1.7112105416372722e-256 < y0 < 9.717693116775471e-203Initial program 27.0
rmApplied add-cube-cbrt27.1
Applied associate-*l*27.1
Taylor expanded around -inf 27.6
if 9.717693116775471e-203 < y0 < 4.308997712685301e-71Initial program 25.3
rmApplied add-cube-cbrt25.4
Taylor expanded around 0 28.3
if 4.308997712685301e-71 < y0 < 1.3561933432615603e+56Initial program 23.5
Taylor expanded around -inf 25.9
Final simplification27.8
herbie shell --seed 2019143
(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)))))