\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}\;y3 \le -6.6136655075287989 \cdot 10^{-142}:\\
\;\;\;\;\left(\left(\left(\left(\left(t \cdot \left(i \cdot \left(z \cdot c\right)\right) - \left(i \cdot \left(c \cdot \left(y \cdot x\right)\right) + a \cdot \left(t \cdot \left(z \cdot b\right)\right)\right)\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)\\
\mathbf{elif}\;y3 \le -4.81816913385106046 \cdot 10^{-224}:\\
\;\;\;\;\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(a \cdot \left(y3 \cdot \left(y1 \cdot z\right)\right) - \left(y0 \cdot \left(z \cdot \left(y3 \cdot c\right)\right) + a \cdot \left(x \cdot \left(y2 \cdot y1\right)\right)\right)\right)\right) + \left(k \cdot \left(i \cdot \left(y \cdot y5\right)\right) - \left(t \cdot \left(i \cdot \left(j \cdot y5\right)\right) + k \cdot \left(y4 \cdot \left(y \cdot b\right)\right)\right)\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)\\
\mathbf{elif}\;y3 \le -1.6540595938942003 \cdot 10^{-263}:\\
\;\;\;\;\left(\left(\left(\left(\left(t \cdot \left(i \cdot \left(z \cdot c\right)\right) - \left(i \cdot \left(c \cdot \left(y \cdot x\right)\right) + a \cdot \left(t \cdot \left(z \cdot b\right)\right)\right)\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)\\
\mathbf{elif}\;y3 \le 8.53527146760082345 \cdot 10^{-262}:\\
\;\;\;\;\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) + 0\\
\mathbf{elif}\;y3 \le 4.187276528621237 \cdot 10^{-242}:\\
\;\;\;\;\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) - 0\right) + \left(k \cdot y2 - j \cdot y3\right) \cdot \left(y4 \cdot y1 - y5 \cdot y0\right)\\
\mathbf{elif}\;y3 \le 4.0396091635139083 \cdot 10^{-62}:\\
\;\;\;\;\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) + 0\\
\mathbf{elif}\;y3 \le 1.17768737819116734 \cdot 10^{86}:\\
\;\;\;\;\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(a \cdot \left(y3 \cdot \left(y \cdot y5\right)\right) - \left(y \cdot \left(y3 \cdot \left(y4 \cdot c\right)\right) + y5 \cdot \left(a \cdot \left(y2 \cdot t\right)\right)\right)\right)\right) + \left(k \cdot y2 - j \cdot y3\right) \cdot \left(y4 \cdot y1 - y5 \cdot y0\right)\\
\mathbf{elif}\;y3 \le 1.13811263110995083 \cdot 10^{133}:\\
\;\;\;\;\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) + 0\right) + \left(k \cdot \left(i \cdot \left(y \cdot y5\right)\right) - \left(t \cdot \left(i \cdot \left(j \cdot y5\right)\right) + k \cdot \left(y4 \cdot \left(y \cdot b\right)\right)\right)\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)\\
\mathbf{else}:\\
\;\;\;\;\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(k \cdot \left(i \cdot \left(y \cdot y5\right)\right) - \left(t \cdot \left(i \cdot \left(j \cdot y5\right)\right) + k \cdot \left(y4 \cdot \left(y \cdot b\right)\right)\right)\right)\right) - \left(t \cdot y2 - y \cdot y3\right) \cdot \left(y4 \cdot c - y5 \cdot a\right)\right) + \left(y0 \cdot \left(y3 \cdot \left(j \cdot y5\right)\right) - \left(y0 \cdot \left(y2 \cdot \left(k \cdot y5\right)\right) + y1 \cdot \left(y3 \cdot \left(j \cdot y4\right)\right)\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 r148229 = x;
double r148230 = y;
double r148231 = r148229 * r148230;
double r148232 = z;
double r148233 = t;
double r148234 = r148232 * r148233;
double r148235 = r148231 - r148234;
double r148236 = a;
double r148237 = b;
double r148238 = r148236 * r148237;
double r148239 = c;
double r148240 = i;
double r148241 = r148239 * r148240;
double r148242 = r148238 - r148241;
double r148243 = r148235 * r148242;
double r148244 = j;
double r148245 = r148229 * r148244;
double r148246 = k;
double r148247 = r148232 * r148246;
double r148248 = r148245 - r148247;
double r148249 = y0;
double r148250 = r148249 * r148237;
double r148251 = y1;
double r148252 = r148251 * r148240;
double r148253 = r148250 - r148252;
double r148254 = r148248 * r148253;
double r148255 = r148243 - r148254;
double r148256 = y2;
double r148257 = r148229 * r148256;
double r148258 = y3;
double r148259 = r148232 * r148258;
double r148260 = r148257 - r148259;
double r148261 = r148249 * r148239;
double r148262 = r148251 * r148236;
double r148263 = r148261 - r148262;
double r148264 = r148260 * r148263;
double r148265 = r148255 + r148264;
double r148266 = r148233 * r148244;
double r148267 = r148230 * r148246;
double r148268 = r148266 - r148267;
double r148269 = y4;
double r148270 = r148269 * r148237;
double r148271 = y5;
double r148272 = r148271 * r148240;
double r148273 = r148270 - r148272;
double r148274 = r148268 * r148273;
double r148275 = r148265 + r148274;
double r148276 = r148233 * r148256;
double r148277 = r148230 * r148258;
double r148278 = r148276 - r148277;
double r148279 = r148269 * r148239;
double r148280 = r148271 * r148236;
double r148281 = r148279 - r148280;
double r148282 = r148278 * r148281;
double r148283 = r148275 - r148282;
double r148284 = r148246 * r148256;
double r148285 = r148244 * r148258;
double r148286 = r148284 - r148285;
double r148287 = r148269 * r148251;
double r148288 = r148271 * r148249;
double r148289 = r148287 - r148288;
double r148290 = r148286 * r148289;
double r148291 = r148283 + r148290;
return r148291;
}
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 r148292 = y3;
double r148293 = -6.613665507528799e-142;
bool r148294 = r148292 <= r148293;
double r148295 = t;
double r148296 = i;
double r148297 = z;
double r148298 = c;
double r148299 = r148297 * r148298;
double r148300 = r148296 * r148299;
double r148301 = r148295 * r148300;
double r148302 = y;
double r148303 = x;
double r148304 = r148302 * r148303;
double r148305 = r148298 * r148304;
double r148306 = r148296 * r148305;
double r148307 = a;
double r148308 = b;
double r148309 = r148297 * r148308;
double r148310 = r148295 * r148309;
double r148311 = r148307 * r148310;
double r148312 = r148306 + r148311;
double r148313 = r148301 - r148312;
double r148314 = j;
double r148315 = r148303 * r148314;
double r148316 = k;
double r148317 = r148297 * r148316;
double r148318 = r148315 - r148317;
double r148319 = y0;
double r148320 = r148319 * r148308;
double r148321 = y1;
double r148322 = r148321 * r148296;
double r148323 = r148320 - r148322;
double r148324 = r148318 * r148323;
double r148325 = r148313 - r148324;
double r148326 = y2;
double r148327 = r148303 * r148326;
double r148328 = r148297 * r148292;
double r148329 = r148327 - r148328;
double r148330 = r148319 * r148298;
double r148331 = r148321 * r148307;
double r148332 = r148330 - r148331;
double r148333 = r148329 * r148332;
double r148334 = r148325 + r148333;
double r148335 = r148295 * r148314;
double r148336 = r148302 * r148316;
double r148337 = r148335 - r148336;
double r148338 = y4;
double r148339 = r148338 * r148308;
double r148340 = y5;
double r148341 = r148340 * r148296;
double r148342 = r148339 - r148341;
double r148343 = r148337 * r148342;
double r148344 = r148334 + r148343;
double r148345 = r148295 * r148326;
double r148346 = r148302 * r148292;
double r148347 = r148345 - r148346;
double r148348 = r148338 * r148298;
double r148349 = r148340 * r148307;
double r148350 = r148348 - r148349;
double r148351 = r148347 * r148350;
double r148352 = r148344 - r148351;
double r148353 = r148316 * r148326;
double r148354 = r148314 * r148292;
double r148355 = r148353 - r148354;
double r148356 = r148338 * r148321;
double r148357 = r148340 * r148319;
double r148358 = r148356 - r148357;
double r148359 = r148355 * r148358;
double r148360 = r148352 + r148359;
double r148361 = -4.81816913385106e-224;
bool r148362 = r148292 <= r148361;
double r148363 = r148303 * r148302;
double r148364 = r148297 * r148295;
double r148365 = r148363 - r148364;
double r148366 = r148307 * r148308;
double r148367 = r148298 * r148296;
double r148368 = r148366 - r148367;
double r148369 = r148365 * r148368;
double r148370 = r148369 - r148324;
double r148371 = r148321 * r148297;
double r148372 = r148292 * r148371;
double r148373 = r148307 * r148372;
double r148374 = r148292 * r148298;
double r148375 = r148297 * r148374;
double r148376 = r148319 * r148375;
double r148377 = r148326 * r148321;
double r148378 = r148303 * r148377;
double r148379 = r148307 * r148378;
double r148380 = r148376 + r148379;
double r148381 = r148373 - r148380;
double r148382 = r148370 + r148381;
double r148383 = r148302 * r148340;
double r148384 = r148296 * r148383;
double r148385 = r148316 * r148384;
double r148386 = r148314 * r148340;
double r148387 = r148296 * r148386;
double r148388 = r148295 * r148387;
double r148389 = r148302 * r148308;
double r148390 = r148338 * r148389;
double r148391 = r148316 * r148390;
double r148392 = r148388 + r148391;
double r148393 = r148385 - r148392;
double r148394 = r148382 + r148393;
double r148395 = r148394 - r148351;
double r148396 = r148395 + r148359;
double r148397 = -1.6540595938942003e-263;
bool r148398 = r148292 <= r148397;
double r148399 = 8.535271467600823e-262;
bool r148400 = r148292 <= r148399;
double r148401 = r148370 + r148333;
double r148402 = r148401 + r148343;
double r148403 = r148402 - r148351;
double r148404 = 0.0;
double r148405 = r148403 + r148404;
double r148406 = 4.187276528621237e-242;
bool r148407 = r148292 <= r148406;
double r148408 = r148402 - r148404;
double r148409 = r148408 + r148359;
double r148410 = 4.039609163513908e-62;
bool r148411 = r148292 <= r148410;
double r148412 = 1.1776873781911673e+86;
bool r148413 = r148292 <= r148412;
double r148414 = r148292 * r148383;
double r148415 = r148307 * r148414;
double r148416 = r148292 * r148348;
double r148417 = r148302 * r148416;
double r148418 = r148326 * r148295;
double r148419 = r148307 * r148418;
double r148420 = r148340 * r148419;
double r148421 = r148417 + r148420;
double r148422 = r148415 - r148421;
double r148423 = r148402 - r148422;
double r148424 = r148423 + r148359;
double r148425 = 1.1381126311099508e+133;
bool r148426 = r148292 <= r148425;
double r148427 = r148370 + r148404;
double r148428 = r148427 + r148393;
double r148429 = r148428 - r148351;
double r148430 = r148429 + r148359;
double r148431 = r148401 + r148393;
double r148432 = r148431 - r148351;
double r148433 = r148292 * r148386;
double r148434 = r148319 * r148433;
double r148435 = r148316 * r148340;
double r148436 = r148326 * r148435;
double r148437 = r148319 * r148436;
double r148438 = r148314 * r148338;
double r148439 = r148292 * r148438;
double r148440 = r148321 * r148439;
double r148441 = r148437 + r148440;
double r148442 = r148434 - r148441;
double r148443 = r148432 + r148442;
double r148444 = r148426 ? r148430 : r148443;
double r148445 = r148413 ? r148424 : r148444;
double r148446 = r148411 ? r148405 : r148445;
double r148447 = r148407 ? r148409 : r148446;
double r148448 = r148400 ? r148405 : r148447;
double r148449 = r148398 ? r148360 : r148448;
double r148450 = r148362 ? r148396 : r148449;
double r148451 = r148294 ? r148360 : r148450;
return r148451;
}



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 y3 < -6.613665507528799e-142 or -4.81816913385106e-224 < y3 < -1.6540595938942003e-263Initial program 27.1
Taylor expanded around inf 28.8
if -6.613665507528799e-142 < y3 < -4.81816913385106e-224Initial program 25.3
Taylor expanded around inf 30.1
Taylor expanded around inf 32.4
if -1.6540595938942003e-263 < y3 < 8.535271467600823e-262 or 4.187276528621237e-242 < y3 < 4.039609163513908e-62Initial program 26.9
Taylor expanded around 0 29.8
if 8.535271467600823e-262 < y3 < 4.187276528621237e-242Initial program 27.4
Taylor expanded around 0 31.5
if 4.039609163513908e-62 < y3 < 1.1776873781911673e+86Initial program 24.0
Taylor expanded around inf 26.3
if 1.1776873781911673e+86 < y3 < 1.1381126311099508e+133Initial program 26.5
Taylor expanded around inf 28.4
Taylor expanded around 0 34.9
if 1.1381126311099508e+133 < y3 Initial program 32.8
Taylor expanded around inf 32.0
Taylor expanded around inf 32.6
Final simplification29.7
herbie shell --seed 2020100
(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"
:precision binary64
(+ (- (+ (+ (- (* (- (* 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)))))