Linear.Matrix:det44 from linear-1.19.1.3

Average Error: 25.0 → 28.6

Time: 1.7m

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}\;b \le -2.1348004490509246 \cdot 10^{-45}:\\ \;\;\;\;\left(y4 \cdot y1 - y5 \cdot y0\right) \cdot \left(y2 \cdot k - j \cdot y3\right) + \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(b \cdot a - i \cdot c\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right)\right) + \left(j \cdot t - y \cdot k\right) \cdot \left(b \cdot y4 - i \cdot y5\right)\right)\\ \mathbf{elif}\;b \le -4.1054470739167587 \cdot 10^{-112}:\\ \;\;\;\;\left(y4 \cdot y1 - y5 \cdot y0\right) \cdot \left(y2 \cdot k - j \cdot y3\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(b \cdot a - i \cdot c\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right)\right) + \left(i \cdot \left(\left(k \cdot y5\right) \cdot y\right) - \left(k \cdot \left(\left(b \cdot y4\right) \cdot y\right) + \left(\left(j \cdot y5\right) \cdot i\right) \cdot t\right)\right)\right) - \left(y4 \cdot c - y5 \cdot a\right) \cdot \left(t \cdot y2 - y \cdot y3\right)\right)\\ \mathbf{elif}\;b \le 1.5590799565832696 \cdot 10^{-117}:\\ \;\;\;\;\left(y4 \cdot y1 - y5 \cdot y0\right) \cdot \left(y2 \cdot k - j \cdot y3\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(b \cdot a - i \cdot c\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right)\right) + \left(j \cdot t - y \cdot k\right) \cdot \left(b \cdot y4 - i \cdot y5\right)\right) - \left(a \cdot \left(y3 \cdot \left(y \cdot y5\right)\right) - \left(y3 \cdot \left(y4 \cdot \left(y \cdot c\right)\right) + a \cdot \left(y2 \cdot \left(y5 \cdot t\right)\right)\right)\right)\right)\\ \mathbf{elif}\;b \le 2.572726652319361 \cdot 10^{-85}:\\ \;\;\;\;\left(y4 \cdot y1 - y5 \cdot y0\right) \cdot \left(y2 \cdot k - j \cdot y3\right) + \left(\left(\left(j \cdot t - y \cdot k\right) \cdot \left(b \cdot y4 - i \cdot y5\right) + \left(\left(x \cdot y - t \cdot z\right) \cdot \left(b \cdot a - i \cdot c\right) + \left(c \cdot y0 - y1 \cdot a\right) \cdot \left(x \cdot y2 - z \cdot y3\right)\right)\right) - \left(y4 \cdot c - y5 \cdot a\right) \cdot \left(t \cdot y2 - y \cdot y3\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\left(y4 \cdot y1 - y5 \cdot y0\right) \cdot \left(y2 \cdot k - j \cdot y3\right) + \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(b \cdot a - i \cdot c\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right)\right) + \left(j \cdot t - y \cdot k\right) \cdot \left(b \cdot y4 - i \cdot y5\right)\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 r2559004 = x;
        double r2559005 = y;
        double r2559006 = r2559004 * r2559005;
        double r2559007 = z;
        double r2559008 = t;
        double r2559009 = r2559007 * r2559008;
        double r2559010 = r2559006 - r2559009;
        double r2559011 = a;
        double r2559012 = b;
        double r2559013 = r2559011 * r2559012;
        double r2559014 = c;
        double r2559015 = i;
        double r2559016 = r2559014 * r2559015;
        double r2559017 = r2559013 - r2559016;
        double r2559018 = r2559010 * r2559017;
        double r2559019 = j;
        double r2559020 = r2559004 * r2559019;
        double r2559021 = k;
        double r2559022 = r2559007 * r2559021;
        double r2559023 = r2559020 - r2559022;
        double r2559024 = y0;
        double r2559025 = r2559024 * r2559012;
        double r2559026 = y1;
        double r2559027 = r2559026 * r2559015;
        double r2559028 = r2559025 - r2559027;
        double r2559029 = r2559023 * r2559028;
        double r2559030 = r2559018 - r2559029;
        double r2559031 = y2;
        double r2559032 = r2559004 * r2559031;
        double r2559033 = y3;
        double r2559034 = r2559007 * r2559033;
        double r2559035 = r2559032 - r2559034;
        double r2559036 = r2559024 * r2559014;
        double r2559037 = r2559026 * r2559011;
        double r2559038 = r2559036 - r2559037;
        double r2559039 = r2559035 * r2559038;
        double r2559040 = r2559030 + r2559039;
        double r2559041 = r2559008 * r2559019;
        double r2559042 = r2559005 * r2559021;
        double r2559043 = r2559041 - r2559042;
        double r2559044 = y4;
        double r2559045 = r2559044 * r2559012;
        double r2559046 = y5;
        double r2559047 = r2559046 * r2559015;
        double r2559048 = r2559045 - r2559047;
        double r2559049 = r2559043 * r2559048;
        double r2559050 = r2559040 + r2559049;
        double r2559051 = r2559008 * r2559031;
        double r2559052 = r2559005 * r2559033;
        double r2559053 = r2559051 - r2559052;
        double r2559054 = r2559044 * r2559014;
        double r2559055 = r2559046 * r2559011;
        double r2559056 = r2559054 - r2559055;
        double r2559057 = r2559053 * r2559056;
        double r2559058 = r2559050 - r2559057;
        double r2559059 = r2559021 * r2559031;
        double r2559060 = r2559019 * r2559033;
        double r2559061 = r2559059 - r2559060;
        double r2559062 = r2559044 * r2559026;
        double r2559063 = r2559046 * r2559024;
        double r2559064 = r2559062 - r2559063;
        double r2559065 = r2559061 * r2559064;
        double r2559066 = r2559058 + r2559065;
        return r2559066;
}

↓

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 r2559067 = b;
        double r2559068 = -2.1348004490509246e-45;
        bool r2559069 = r2559067 <= r2559068;
        double r2559070 = y4;
        double r2559071 = y1;
        double r2559072 = r2559070 * r2559071;
        double r2559073 = y5;
        double r2559074 = y0;
        double r2559075 = r2559073 * r2559074;
        double r2559076 = r2559072 - r2559075;
        double r2559077 = y2;
        double r2559078 = k;
        double r2559079 = r2559077 * r2559078;
        double r2559080 = j;
        double r2559081 = y3;
        double r2559082 = r2559080 * r2559081;
        double r2559083 = r2559079 - r2559082;
        double r2559084 = r2559076 * r2559083;
        double r2559085 = c;
        double r2559086 = r2559085 * r2559074;
        double r2559087 = a;
        double r2559088 = r2559071 * r2559087;
        double r2559089 = r2559086 - r2559088;
        double r2559090 = x;
        double r2559091 = r2559090 * r2559077;
        double r2559092 = z;
        double r2559093 = r2559092 * r2559081;
        double r2559094 = r2559091 - r2559093;
        double r2559095 = r2559089 * r2559094;
        double r2559096 = y;
        double r2559097 = r2559090 * r2559096;
        double r2559098 = t;
        double r2559099 = r2559098 * r2559092;
        double r2559100 = r2559097 - r2559099;
        double r2559101 = r2559067 * r2559087;
        double r2559102 = i;
        double r2559103 = r2559102 * r2559085;
        double r2559104 = r2559101 - r2559103;
        double r2559105 = r2559100 * r2559104;
        double r2559106 = r2559090 * r2559080;
        double r2559107 = r2559092 * r2559078;
        double r2559108 = r2559106 - r2559107;
        double r2559109 = r2559067 * r2559074;
        double r2559110 = r2559102 * r2559071;
        double r2559111 = r2559109 - r2559110;
        double r2559112 = r2559108 * r2559111;
        double r2559113 = r2559105 - r2559112;
        double r2559114 = r2559095 + r2559113;
        double r2559115 = r2559080 * r2559098;
        double r2559116 = r2559096 * r2559078;
        double r2559117 = r2559115 - r2559116;
        double r2559118 = r2559067 * r2559070;
        double r2559119 = r2559102 * r2559073;
        double r2559120 = r2559118 - r2559119;
        double r2559121 = r2559117 * r2559120;
        double r2559122 = r2559114 + r2559121;
        double r2559123 = r2559084 + r2559122;
        double r2559124 = -4.1054470739167587e-112;
        bool r2559125 = r2559067 <= r2559124;
        double r2559126 = r2559078 * r2559073;
        double r2559127 = r2559126 * r2559096;
        double r2559128 = r2559102 * r2559127;
        double r2559129 = r2559118 * r2559096;
        double r2559130 = r2559078 * r2559129;
        double r2559131 = r2559080 * r2559073;
        double r2559132 = r2559131 * r2559102;
        double r2559133 = r2559132 * r2559098;
        double r2559134 = r2559130 + r2559133;
        double r2559135 = r2559128 - r2559134;
        double r2559136 = r2559114 + r2559135;
        double r2559137 = r2559070 * r2559085;
        double r2559138 = r2559073 * r2559087;
        double r2559139 = r2559137 - r2559138;
        double r2559140 = r2559098 * r2559077;
        double r2559141 = r2559096 * r2559081;
        double r2559142 = r2559140 - r2559141;
        double r2559143 = r2559139 * r2559142;
        double r2559144 = r2559136 - r2559143;
        double r2559145 = r2559084 + r2559144;
        double r2559146 = 1.5590799565832696e-117;
        bool r2559147 = r2559067 <= r2559146;
        double r2559148 = r2559096 * r2559073;
        double r2559149 = r2559081 * r2559148;
        double r2559150 = r2559087 * r2559149;
        double r2559151 = r2559096 * r2559085;
        double r2559152 = r2559070 * r2559151;
        double r2559153 = r2559081 * r2559152;
        double r2559154 = r2559073 * r2559098;
        double r2559155 = r2559077 * r2559154;
        double r2559156 = r2559087 * r2559155;
        double r2559157 = r2559153 + r2559156;
        double r2559158 = r2559150 - r2559157;
        double r2559159 = r2559122 - r2559158;
        double r2559160 = r2559084 + r2559159;
        double r2559161 = 2.572726652319361e-85;
        bool r2559162 = r2559067 <= r2559161;
        double r2559163 = r2559105 + r2559095;
        double r2559164 = r2559121 + r2559163;
        double r2559165 = r2559164 - r2559143;
        double r2559166 = r2559084 + r2559165;
        double r2559167 = r2559162 ? r2559166 : r2559123;
        double r2559168 = r2559147 ? r2559160 : r2559167;
        double r2559169 = r2559125 ? r2559145 : r2559168;
        double r2559170 = r2559069 ? r2559123 : r2559169;
        return r2559170;
}

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 4 regimes

2. if b < -2.1348004490509246e-45 or 2.572726652319361e-85 < b

   1. Initial program 25.7

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

   if -2.1348004490509246e-45 < b < -4.1054470739167587e-112

   1. Initial program 20.7

      \[\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 23.3

      \[\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) + \color{blue}{\left(i \cdot \left(y \cdot \left(y5 \cdot k\right)\right) - \left(k \cdot \left(y \cdot \left(b \cdot y4\right)\right) + t \cdot \left(i \cdot \left(j \cdot y5\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)\]

   if -4.1054470739167587e-112 < b < 1.5590799565832696e-117

   1. Initial program 25.2

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

   if 1.5590799565832696e-117 < b < 2.572726652319361e-85

   1. Initial program 22.1

      \[\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 0 29.5

      \[\leadsto \left(\left(\left(\left(\left(x \cdot y - z \cdot t\right) \cdot \left(a \cdot b - c \cdot i\right) - \color{blue}{0}\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. Recombined 4 regimes into one program.

4. Final simplification 28.6

   \[\leadsto \begin{array}{l} \mathbf{if}\;b \le -2.1348004490509246 \cdot 10^{-45}:\\ \;\;\;\;\left(y4 \cdot y1 - y5 \cdot y0\right) \cdot \left(y2 \cdot k - j \cdot y3\right) + \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(b \cdot a - i \cdot c\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right)\right) + \left(j \cdot t - y \cdot k\right) \cdot \left(b \cdot y4 - i \cdot y5\right)\right)\\ \mathbf{elif}\;b \le -4.1054470739167587 \cdot 10^{-112}:\\ \;\;\;\;\left(y4 \cdot y1 - y5 \cdot y0\right) \cdot \left(y2 \cdot k - j \cdot y3\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(b \cdot a - i \cdot c\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right)\right) + \left(i \cdot \left(\left(k \cdot y5\right) \cdot y\right) - \left(k \cdot \left(\left(b \cdot y4\right) \cdot y\right) + \left(\left(j \cdot y5\right) \cdot i\right) \cdot t\right)\right)\right) - \left(y4 \cdot c - y5 \cdot a\right) \cdot \left(t \cdot y2 - y \cdot y3\right)\right)\\ \mathbf{elif}\;b \le 1.5590799565832696 \cdot 10^{-117}:\\ \;\;\;\;\left(y4 \cdot y1 - y5 \cdot y0\right) \cdot \left(y2 \cdot k - j \cdot y3\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(b \cdot a - i \cdot c\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right)\right) + \left(j \cdot t - y \cdot k\right) \cdot \left(b \cdot y4 - i \cdot y5\right)\right) - \left(a \cdot \left(y3 \cdot \left(y \cdot y5\right)\right) - \left(y3 \cdot \left(y4 \cdot \left(y \cdot c\right)\right) + a \cdot \left(y2 \cdot \left(y5 \cdot t\right)\right)\right)\right)\right)\\ \mathbf{elif}\;b \le 2.572726652319361 \cdot 10^{-85}:\\ \;\;\;\;\left(y4 \cdot y1 - y5 \cdot y0\right) \cdot \left(y2 \cdot k - j \cdot y3\right) + \left(\left(\left(j \cdot t - y \cdot k\right) \cdot \left(b \cdot y4 - i \cdot y5\right) + \left(\left(x \cdot y - t \cdot z\right) \cdot \left(b \cdot a - i \cdot c\right) + \left(c \cdot y0 - y1 \cdot a\right) \cdot \left(x \cdot y2 - z \cdot y3\right)\right)\right) - \left(y4 \cdot c - y5 \cdot a\right) \cdot \left(t \cdot y2 - y \cdot y3\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\left(y4 \cdot y1 - y5 \cdot y0\right) \cdot \left(y2 \cdot k - j \cdot y3\right) + \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(b \cdot a - i \cdot c\right) - \left(x \cdot j - z \cdot k\right) \cdot \left(b \cdot y0 - i \cdot y1\right)\right)\right) + \left(j \cdot t - y \cdot k\right) \cdot \left(b \cdot y4 - i \cdot y5\right)\right)\\ \end{array}\]

Reproduce

herbie shell --seed 2019154 
(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)))))