Linear.Matrix:det44 from linear-1.19.1.3

Average Error: 25.7 → 26.7

Time: 5.1m

Precision: 64

Math

\[\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 -5.80942904443867 \cdot 10^{-106}:\\ \;\;\;\;\left(y4 \cdot y1 - y0 \cdot y5\right) \cdot \left(k \cdot y2 - j \cdot y3\right) + \left(\left(\left(j \cdot t - k \cdot y\right) \cdot \left(b \cdot y4 - i \cdot y5\right) + \left(\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right) + \left(\mathsf{fma}\left(\left(x \cdot \left(c \cdot i\right)\right), \left(-y\right), \left(\left(t \cdot \left(c \cdot i - b \cdot a\right)\right) \cdot z\right)\right) - \left(b \cdot y0 - i \cdot y1\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right)\right) - \left(t \cdot y2 - y3 \cdot y\right) \cdot \left(c \cdot y4 - a \cdot y5\right)\right)\\ \mathbf{elif}\;y0 \le -2.015611789438149 \cdot 10^{-147}:\\ \;\;\;\;\left(\left(\left(j \cdot t - k \cdot y\right) \cdot \left(b \cdot y4 - i \cdot y5\right) + \left(\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right) + \left(\left(b \cdot a - c \cdot i\right) \cdot \left(y \cdot x - z \cdot t\right) - \left(b \cdot y0 - i \cdot y1\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right)\right) - \left(t \cdot y2 - y3 \cdot y\right) \cdot \left(c \cdot y4 - a \cdot y5\right)\right) + \left(\left(\left(y0 \cdot y5 - y4 \cdot y1\right) \cdot y3\right) \cdot j - \left(\left(y0 \cdot y5\right) \cdot y2\right) \cdot k\right)\\ \mathbf{elif}\;y0 \le 1.2019899616792464 \cdot 10^{-228}:\\ \;\;\;\;\left(y4 \cdot y1 - y0 \cdot y5\right) \cdot \left(k \cdot y2 - j \cdot y3\right) + \left(\left(\left(j \cdot t - k \cdot y\right) \cdot \left(b \cdot y4 - i \cdot y5\right) + \left(\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right) + \left(\mathsf{fma}\left(\left(x \cdot \left(c \cdot i\right)\right), \left(-y\right), \left(\left(t \cdot \left(c \cdot i - b \cdot a\right)\right) \cdot z\right)\right) - \left(b \cdot y0 - i \cdot y1\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right)\right) - \left(t \cdot y2 - y3 \cdot y\right) \cdot \left(c \cdot y4 - a \cdot y5\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\sqrt[3]{\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right)} \cdot \left(\sqrt[3]{\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right)} \cdot \sqrt[3]{\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right)}\right) + \left(\left(b \cdot a - c \cdot i\right) \cdot \left(y \cdot x - z \cdot t\right) - \left(b \cdot y0 - i \cdot y1\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right) + \left(j \cdot t - k \cdot y\right) \cdot \left(b \cdot y4 - i \cdot y5\right)\right) - \left(t \cdot y2 - y3 \cdot y\right) \cdot \left(c \cdot y4 - a \cdot y5\right)\right) + \left(y4 \cdot y1 - y0 \cdot y5\right) \cdot \left(k \cdot y2 - j \cdot y3\right)\\ \end{array}\]

C

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 r28650294 = x;
        double r28650295 = y;
        double r28650296 = r28650294 * r28650295;
        double r28650297 = z;
        double r28650298 = t;
        double r28650299 = r28650297 * r28650298;
        double r28650300 = r28650296 - r28650299;
        double r28650301 = a;
        double r28650302 = b;
        double r28650303 = r28650301 * r28650302;
        double r28650304 = c;
        double r28650305 = i;
        double r28650306 = r28650304 * r28650305;
        double r28650307 = r28650303 - r28650306;
        double r28650308 = r28650300 * r28650307;
        double r28650309 = j;
        double r28650310 = r28650294 * r28650309;
        double r28650311 = k;
        double r28650312 = r28650297 * r28650311;
        double r28650313 = r28650310 - r28650312;
        double r28650314 = y0;
        double r28650315 = r28650314 * r28650302;
        double r28650316 = y1;
        double r28650317 = r28650316 * r28650305;
        double r28650318 = r28650315 - r28650317;
        double r28650319 = r28650313 * r28650318;
        double r28650320 = r28650308 - r28650319;
        double r28650321 = y2;
        double r28650322 = r28650294 * r28650321;
        double r28650323 = y3;
        double r28650324 = r28650297 * r28650323;
        double r28650325 = r28650322 - r28650324;
        double r28650326 = r28650314 * r28650304;
        double r28650327 = r28650316 * r28650301;
        double r28650328 = r28650326 - r28650327;
        double r28650329 = r28650325 * r28650328;
        double r28650330 = r28650320 + r28650329;
        double r28650331 = r28650298 * r28650309;
        double r28650332 = r28650295 * r28650311;
        double r28650333 = r28650331 - r28650332;
        double r28650334 = y4;
        double r28650335 = r28650334 * r28650302;
        double r28650336 = y5;
        double r28650337 = r28650336 * r28650305;
        double r28650338 = r28650335 - r28650337;
        double r28650339 = r28650333 * r28650338;
        double r28650340 = r28650330 + r28650339;
        double r28650341 = r28650298 * r28650321;
        double r28650342 = r28650295 * r28650323;
        double r28650343 = r28650341 - r28650342;
        double r28650344 = r28650334 * r28650304;
        double r28650345 = r28650336 * r28650301;
        double r28650346 = r28650344 - r28650345;
        double r28650347 = r28650343 * r28650346;
        double r28650348 = r28650340 - r28650347;
        double r28650349 = r28650311 * r28650321;
        double r28650350 = r28650309 * r28650323;
        double r28650351 = r28650349 - r28650350;
        double r28650352 = r28650334 * r28650316;
        double r28650353 = r28650336 * r28650314;
        double r28650354 = r28650352 - r28650353;
        double r28650355 = r28650351 * r28650354;
        double r28650356 = r28650348 + r28650355;
        return r28650356;
}

↓

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 r28650357 = y0;
        double r28650358 = -5.80942904443867e-106;
        bool r28650359 = r28650357 <= r28650358;
        double r28650360 = y4;
        double r28650361 = y1;
        double r28650362 = r28650360 * r28650361;
        double r28650363 = y5;
        double r28650364 = r28650357 * r28650363;
        double r28650365 = r28650362 - r28650364;
        double r28650366 = k;
        double r28650367 = y2;
        double r28650368 = r28650366 * r28650367;
        double r28650369 = j;
        double r28650370 = y3;
        double r28650371 = r28650369 * r28650370;
        double r28650372 = r28650368 - r28650371;
        double r28650373 = r28650365 * r28650372;
        double r28650374 = t;
        double r28650375 = r28650369 * r28650374;
        double r28650376 = y;
        double r28650377 = r28650366 * r28650376;
        double r28650378 = r28650375 - r28650377;
        double r28650379 = b;
        double r28650380 = r28650379 * r28650360;
        double r28650381 = i;
        double r28650382 = r28650381 * r28650363;
        double r28650383 = r28650380 - r28650382;
        double r28650384 = r28650378 * r28650383;
        double r28650385 = c;
        double r28650386 = r28650385 * r28650357;
        double r28650387 = a;
        double r28650388 = r28650387 * r28650361;
        double r28650389 = r28650386 - r28650388;
        double r28650390 = x;
        double r28650391 = r28650367 * r28650390;
        double r28650392 = z;
        double r28650393 = r28650392 * r28650370;
        double r28650394 = r28650391 - r28650393;
        double r28650395 = r28650389 * r28650394;
        double r28650396 = r28650385 * r28650381;
        double r28650397 = r28650390 * r28650396;
        double r28650398 = -r28650376;
        double r28650399 = r28650379 * r28650387;
        double r28650400 = r28650396 - r28650399;
        double r28650401 = r28650374 * r28650400;
        double r28650402 = r28650401 * r28650392;
        double r28650403 = fma(r28650397, r28650398, r28650402);
        double r28650404 = r28650379 * r28650357;
        double r28650405 = r28650381 * r28650361;
        double r28650406 = r28650404 - r28650405;
        double r28650407 = r28650390 * r28650369;
        double r28650408 = r28650392 * r28650366;
        double r28650409 = r28650407 - r28650408;
        double r28650410 = r28650406 * r28650409;
        double r28650411 = r28650403 - r28650410;
        double r28650412 = r28650395 + r28650411;
        double r28650413 = r28650384 + r28650412;
        double r28650414 = r28650374 * r28650367;
        double r28650415 = r28650370 * r28650376;
        double r28650416 = r28650414 - r28650415;
        double r28650417 = r28650385 * r28650360;
        double r28650418 = r28650387 * r28650363;
        double r28650419 = r28650417 - r28650418;
        double r28650420 = r28650416 * r28650419;
        double r28650421 = r28650413 - r28650420;
        double r28650422 = r28650373 + r28650421;
        double r28650423 = -2.015611789438149e-147;
        bool r28650424 = r28650357 <= r28650423;
        double r28650425 = r28650399 - r28650396;
        double r28650426 = r28650376 * r28650390;
        double r28650427 = r28650392 * r28650374;
        double r28650428 = r28650426 - r28650427;
        double r28650429 = r28650425 * r28650428;
        double r28650430 = r28650429 - r28650410;
        double r28650431 = r28650395 + r28650430;
        double r28650432 = r28650384 + r28650431;
        double r28650433 = r28650432 - r28650420;
        double r28650434 = r28650364 - r28650362;
        double r28650435 = r28650434 * r28650370;
        double r28650436 = r28650435 * r28650369;
        double r28650437 = r28650364 * r28650367;
        double r28650438 = r28650437 * r28650366;
        double r28650439 = r28650436 - r28650438;
        double r28650440 = r28650433 + r28650439;
        double r28650441 = 1.2019899616792464e-228;
        bool r28650442 = r28650357 <= r28650441;
        double r28650443 = cbrt(r28650395);
        double r28650444 = r28650443 * r28650443;
        double r28650445 = r28650443 * r28650444;
        double r28650446 = r28650445 + r28650430;
        double r28650447 = r28650446 + r28650384;
        double r28650448 = r28650447 - r28650420;
        double r28650449 = r28650448 + r28650373;
        double r28650450 = r28650442 ? r28650422 : r28650449;
        double r28650451 = r28650424 ? r28650440 : r28650450;
        double r28650452 = r28650359 ? r28650422 : r28650451;
        return r28650452;
}

Error

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

Derivation

1. Split input into 3 regimes

2. if y0 < -5.80942904443867e-106 or -2.015611789438149e-147 < y0 < 1.2019899616792464e-228

   1. Initial program 25.4

      \[\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)\]

   2. Taylor expanded around inf 28.1

      \[\leadsto \left(\left(\left(\left(\color{blue}{\left(t \cdot \left(i \cdot \left(z \cdot c\right)\right) - \left(a \cdot \left(t \cdot \left(b \cdot z\right)\right) + i \cdot \left(x \cdot \left(c \cdot y\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)\]

   3. Simplified 27.2

      \[\leadsto \left(\left(\left(\left(\color{blue}{\mathsf{fma}\left(\left(x \cdot \left(c \cdot i\right)\right), \left(-y\right), \left(z \cdot \left(t \cdot \left(c \cdot i - a \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)\]

   if -5.80942904443867e-106 < y0 < -2.015611789438149e-147

   1. Initial program 26.3

      \[\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)\]

   2. Taylor expanded around -inf 27.5

      \[\leadsto \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) + \color{blue}{\left(y3 \cdot \left(j \cdot \left(y5 \cdot y0\right)\right) - \left(y1 \cdot \left(y3 \cdot \left(y4 \cdot j\right)\right) + k \cdot \left(y2 \cdot \left(y5 \cdot y0\right)\right)\right)\right)}\]

   3. Simplified 26.4

      \[\leadsto \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) + \color{blue}{\left(j \cdot \left(y3 \cdot \left(y5 \cdot y0 - y1 \cdot y4\right)\right) - k \cdot \left(\left(y5 \cdot y0\right) \cdot y2\right)\right)}\]

   if 1.2019899616792464e-228 < y0

   1. Initial program 25.9

      \[\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)\]
   2. Using strategy rm

   3. Applied add-cube-cbrt 26.0

      \[\leadsto \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) + \color{blue}{\left(\sqrt[3]{\left(x \cdot y2 - z \cdot y3\right) \cdot \left(y0 \cdot c - y1 \cdot a\right)} \cdot \sqrt[3]{\left(x \cdot y2 - z \cdot y3\right) \cdot \left(y0 \cdot c - y1 \cdot a\right)}\right) \cdot \sqrt[3]{\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)\]
3. Recombined 3 regimes into one program.

4. Final simplification 26.7

   \[\leadsto \begin{array}{l} \mathbf{if}\;y0 \le -5.80942904443867 \cdot 10^{-106}:\\ \;\;\;\;\left(y4 \cdot y1 - y0 \cdot y5\right) \cdot \left(k \cdot y2 - j \cdot y3\right) + \left(\left(\left(j \cdot t - k \cdot y\right) \cdot \left(b \cdot y4 - i \cdot y5\right) + \left(\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right) + \left(\mathsf{fma}\left(\left(x \cdot \left(c \cdot i\right)\right), \left(-y\right), \left(\left(t \cdot \left(c \cdot i - b \cdot a\right)\right) \cdot z\right)\right) - \left(b \cdot y0 - i \cdot y1\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right)\right) - \left(t \cdot y2 - y3 \cdot y\right) \cdot \left(c \cdot y4 - a \cdot y5\right)\right)\\ \mathbf{elif}\;y0 \le -2.015611789438149 \cdot 10^{-147}:\\ \;\;\;\;\left(\left(\left(j \cdot t - k \cdot y\right) \cdot \left(b \cdot y4 - i \cdot y5\right) + \left(\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right) + \left(\left(b \cdot a - c \cdot i\right) \cdot \left(y \cdot x - z \cdot t\right) - \left(b \cdot y0 - i \cdot y1\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right)\right) - \left(t \cdot y2 - y3 \cdot y\right) \cdot \left(c \cdot y4 - a \cdot y5\right)\right) + \left(\left(\left(y0 \cdot y5 - y4 \cdot y1\right) \cdot y3\right) \cdot j - \left(\left(y0 \cdot y5\right) \cdot y2\right) \cdot k\right)\\ \mathbf{elif}\;y0 \le 1.2019899616792464 \cdot 10^{-228}:\\ \;\;\;\;\left(y4 \cdot y1 - y0 \cdot y5\right) \cdot \left(k \cdot y2 - j \cdot y3\right) + \left(\left(\left(j \cdot t - k \cdot y\right) \cdot \left(b \cdot y4 - i \cdot y5\right) + \left(\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right) + \left(\mathsf{fma}\left(\left(x \cdot \left(c \cdot i\right)\right), \left(-y\right), \left(\left(t \cdot \left(c \cdot i - b \cdot a\right)\right) \cdot z\right)\right) - \left(b \cdot y0 - i \cdot y1\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right)\right) - \left(t \cdot y2 - y3 \cdot y\right) \cdot \left(c \cdot y4 - a \cdot y5\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\sqrt[3]{\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right)} \cdot \left(\sqrt[3]{\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right)} \cdot \sqrt[3]{\left(c \cdot y0 - a \cdot y1\right) \cdot \left(y2 \cdot x - z \cdot y3\right)}\right) + \left(\left(b \cdot a - c \cdot i\right) \cdot \left(y \cdot x - z \cdot t\right) - \left(b \cdot y0 - i \cdot y1\right) \cdot \left(x \cdot j - z \cdot k\right)\right)\right) + \left(j \cdot t - k \cdot y\right) \cdot \left(b \cdot y4 - i \cdot y5\right)\right) - \left(t \cdot y2 - y3 \cdot y\right) \cdot \left(c \cdot y4 - a \cdot y5\right)\right) + \left(y4 \cdot y1 - y0 \cdot y5\right) \cdot \left(k \cdot y2 - j \cdot y3\right)\\ \end{array}\]

Reproduce

herbie shell --seed 2019128 +o rules:numerics
(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)))))