Linear.Matrix:det44 from linear-1.19.1.3

Average Error: 25.2 → 27.0

Time: 1.6m

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}\;y2 \le -1.2961541745563326 \cdot 10^{+62}:\\ \;\;\;\;\left(\left(\left(y0 \cdot c - 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(y0 \cdot b - i \cdot y1\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(t \cdot y2 - y \cdot y3\right) \cdot \left(y4 \cdot c - y5 \cdot a\right)\\ \mathbf{elif}\;y2 \le 2.707765602532912 \cdot 10^{-42}:\\ \;\;\;\;\left(\left(\left(\left(\left(y3 \cdot \left(z \cdot y1\right) - x \cdot \left(y2 \cdot y1\right)\right) \cdot a - \left(y3 \cdot \left(c \cdot z\right)\right) \cdot y0\right) + \left(\left(x \cdot y - t \cdot z\right) \cdot \left(a \cdot b - i \cdot c\right) - \left(y0 \cdot b - i \cdot y1\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(t \cdot y2 - y \cdot y3\right) \cdot \left(y4 \cdot c - y5 \cdot a\right)\right) + \left(y1 \cdot y4 - y0 \cdot y5\right) \cdot \left(y2 \cdot k - y3 \cdot j\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(\left(\left(y0 \cdot c - y1 \cdot a\right) \cdot \left(x \cdot y2 - z \cdot y3\right) + \left(\left(t \cdot \left(i \cdot \left(c \cdot z\right)\right) - \left(\left(x \cdot \left(y \cdot c\right)\right) \cdot i + a \cdot \left(t \cdot \left(b \cdot z\right)\right)\right)\right) - \left(y0 \cdot b - i \cdot y1\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(t \cdot y2 - y \cdot y3\right) \cdot \left(y4 \cdot c - y5 \cdot a\right)\right) + \left(y1 \cdot y4 - y0 \cdot y5\right) \cdot \left(y2 \cdot k - y3 \cdot j\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 r29723156 = x;
        double r29723157 = y;
        double r29723158 = r29723156 * r29723157;
        double r29723159 = z;
        double r29723160 = t;
        double r29723161 = r29723159 * r29723160;
        double r29723162 = r29723158 - r29723161;
        double r29723163 = a;
        double r29723164 = b;
        double r29723165 = r29723163 * r29723164;
        double r29723166 = c;
        double r29723167 = i;
        double r29723168 = r29723166 * r29723167;
        double r29723169 = r29723165 - r29723168;
        double r29723170 = r29723162 * r29723169;
        double r29723171 = j;
        double r29723172 = r29723156 * r29723171;
        double r29723173 = k;
        double r29723174 = r29723159 * r29723173;
        double r29723175 = r29723172 - r29723174;
        double r29723176 = y0;
        double r29723177 = r29723176 * r29723164;
        double r29723178 = y1;
        double r29723179 = r29723178 * r29723167;
        double r29723180 = r29723177 - r29723179;
        double r29723181 = r29723175 * r29723180;
        double r29723182 = r29723170 - r29723181;
        double r29723183 = y2;
        double r29723184 = r29723156 * r29723183;
        double r29723185 = y3;
        double r29723186 = r29723159 * r29723185;
        double r29723187 = r29723184 - r29723186;
        double r29723188 = r29723176 * r29723166;
        double r29723189 = r29723178 * r29723163;
        double r29723190 = r29723188 - r29723189;
        double r29723191 = r29723187 * r29723190;
        double r29723192 = r29723182 + r29723191;
        double r29723193 = r29723160 * r29723171;
        double r29723194 = r29723157 * r29723173;
        double r29723195 = r29723193 - r29723194;
        double r29723196 = y4;
        double r29723197 = r29723196 * r29723164;
        double r29723198 = y5;
        double r29723199 = r29723198 * r29723167;
        double r29723200 = r29723197 - r29723199;
        double r29723201 = r29723195 * r29723200;
        double r29723202 = r29723192 + r29723201;
        double r29723203 = r29723160 * r29723183;
        double r29723204 = r29723157 * r29723185;
        double r29723205 = r29723203 - r29723204;
        double r29723206 = r29723196 * r29723166;
        double r29723207 = r29723198 * r29723163;
        double r29723208 = r29723206 - r29723207;
        double r29723209 = r29723205 * r29723208;
        double r29723210 = r29723202 - r29723209;
        double r29723211 = r29723173 * r29723183;
        double r29723212 = r29723171 * r29723185;
        double r29723213 = r29723211 - r29723212;
        double r29723214 = r29723196 * r29723178;
        double r29723215 = r29723198 * r29723176;
        double r29723216 = r29723214 - r29723215;
        double r29723217 = r29723213 * r29723216;
        double r29723218 = r29723210 + r29723217;
        return r29723218;
}

↓

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 r29723219 = y2;
        double r29723220 = -1.2961541745563326e+62;
        bool r29723221 = r29723219 <= r29723220;
        double r29723222 = y0;
        double r29723223 = c;
        double r29723224 = r29723222 * r29723223;
        double r29723225 = y1;
        double r29723226 = a;
        double r29723227 = r29723225 * r29723226;
        double r29723228 = r29723224 - r29723227;
        double r29723229 = x;
        double r29723230 = r29723229 * r29723219;
        double r29723231 = z;
        double r29723232 = y3;
        double r29723233 = r29723231 * r29723232;
        double r29723234 = r29723230 - r29723233;
        double r29723235 = r29723228 * r29723234;
        double r29723236 = y;
        double r29723237 = r29723229 * r29723236;
        double r29723238 = t;
        double r29723239 = r29723238 * r29723231;
        double r29723240 = r29723237 - r29723239;
        double r29723241 = b;
        double r29723242 = r29723226 * r29723241;
        double r29723243 = i;
        double r29723244 = r29723243 * r29723223;
        double r29723245 = r29723242 - r29723244;
        double r29723246 = r29723240 * r29723245;
        double r29723247 = r29723222 * r29723241;
        double r29723248 = r29723243 * r29723225;
        double r29723249 = r29723247 - r29723248;
        double r29723250 = j;
        double r29723251 = r29723229 * r29723250;
        double r29723252 = k;
        double r29723253 = r29723231 * r29723252;
        double r29723254 = r29723251 - r29723253;
        double r29723255 = r29723249 * r29723254;
        double r29723256 = r29723246 - r29723255;
        double r29723257 = r29723235 + r29723256;
        double r29723258 = r29723238 * r29723250;
        double r29723259 = r29723236 * r29723252;
        double r29723260 = r29723258 - r29723259;
        double r29723261 = y4;
        double r29723262 = r29723261 * r29723241;
        double r29723263 = y5;
        double r29723264 = r29723243 * r29723263;
        double r29723265 = r29723262 - r29723264;
        double r29723266 = r29723260 * r29723265;
        double r29723267 = r29723257 + r29723266;
        double r29723268 = r29723238 * r29723219;
        double r29723269 = r29723236 * r29723232;
        double r29723270 = r29723268 - r29723269;
        double r29723271 = r29723261 * r29723223;
        double r29723272 = r29723263 * r29723226;
        double r29723273 = r29723271 - r29723272;
        double r29723274 = r29723270 * r29723273;
        double r29723275 = r29723267 - r29723274;
        double r29723276 = 2.707765602532912e-42;
        bool r29723277 = r29723219 <= r29723276;
        double r29723278 = r29723231 * r29723225;
        double r29723279 = r29723232 * r29723278;
        double r29723280 = r29723219 * r29723225;
        double r29723281 = r29723229 * r29723280;
        double r29723282 = r29723279 - r29723281;
        double r29723283 = r29723282 * r29723226;
        double r29723284 = r29723223 * r29723231;
        double r29723285 = r29723232 * r29723284;
        double r29723286 = r29723285 * r29723222;
        double r29723287 = r29723283 - r29723286;
        double r29723288 = r29723287 + r29723256;
        double r29723289 = r29723288 + r29723266;
        double r29723290 = r29723289 - r29723274;
        double r29723291 = r29723225 * r29723261;
        double r29723292 = r29723222 * r29723263;
        double r29723293 = r29723291 - r29723292;
        double r29723294 = r29723219 * r29723252;
        double r29723295 = r29723232 * r29723250;
        double r29723296 = r29723294 - r29723295;
        double r29723297 = r29723293 * r29723296;
        double r29723298 = r29723290 + r29723297;
        double r29723299 = r29723243 * r29723284;
        double r29723300 = r29723238 * r29723299;
        double r29723301 = r29723236 * r29723223;
        double r29723302 = r29723229 * r29723301;
        double r29723303 = r29723302 * r29723243;
        double r29723304 = r29723241 * r29723231;
        double r29723305 = r29723238 * r29723304;
        double r29723306 = r29723226 * r29723305;
        double r29723307 = r29723303 + r29723306;
        double r29723308 = r29723300 - r29723307;
        double r29723309 = r29723308 - r29723255;
        double r29723310 = r29723235 + r29723309;
        double r29723311 = r29723310 + r29723266;
        double r29723312 = r29723311 - r29723274;
        double r29723313 = r29723312 + r29723297;
        double r29723314 = r29723277 ? r29723298 : r29723313;
        double r29723315 = r29723221 ? r29723275 : r29723314;
        return r29723315;
}

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

Target

Original	25.2
Target	29.3
Herbie	27.0

\[\begin{array}{l} \mathbf{if}\;y4 \lt -7.206256231996481 \cdot 10^{+60}:\\ \;\;\;\;\left(\left(b \cdot a - i \cdot c\right) \cdot \left(y \cdot x - t \cdot z\right) - \left(\left(j \cdot x - k \cdot z\right) \cdot \left(y0 \cdot b - i \cdot y1\right) - \left(j \cdot t - k \cdot y\right) \cdot \left(y4 \cdot b - y5 \cdot i\right)\right)\right) - \left(\frac{y2 \cdot t - y3 \cdot y}{\frac{1}{y4 \cdot c - y5 \cdot a}} - \left(y2 \cdot k - y3 \cdot j\right) \cdot \left(y4 \cdot y1 - y5 \cdot y0\right)\right)\\ \mathbf{elif}\;y4 \lt -3.364603505246317 \cdot 10^{-66}:\\ \;\;\;\;\left(\left(\left(\left(t \cdot c\right) \cdot \left(i \cdot z\right) - \left(a \cdot t\right) \cdot \left(b \cdot z\right)\right) - \left(y \cdot c\right) \cdot \left(i \cdot x\right)\right) - \left(b \cdot y0 - i \cdot y1\right) \cdot \left(j \cdot x - k \cdot z\right)\right) + \left(\left(y0 \cdot c - a \cdot y1\right) \cdot \left(x \cdot y2 - z \cdot y3\right) - \left(\left(t \cdot y2 - y \cdot y3\right) \cdot \left(y4 \cdot c - a \cdot y5\right) - \left(y1 \cdot y4 - y5 \cdot y0\right) \cdot \left(k \cdot y2 - j \cdot y3\right)\right)\right)\\ \mathbf{elif}\;y4 \lt -1.2000065055686116 \cdot 10^{-105}:\\ \;\;\;\;\left(\left(\left(j \cdot t - k \cdot y\right) \cdot \left(y4 \cdot b - y5 \cdot i\right) - \left(y3 \cdot y\right) \cdot \left(y5 \cdot a - y4 \cdot c\right)\right) + \left(\left(y5 \cdot a\right) \cdot \left(t \cdot y2\right) + \left(k \cdot y2 - j \cdot y3\right) \cdot \left(y4 \cdot y1 - y5 \cdot y0\right)\right)\right) + \left(\left(x \cdot y2 - z \cdot y3\right) \cdot \left(c \cdot y0 - a \cdot y1\right) - \left(\left(b \cdot y0 - i \cdot y1\right) \cdot \left(j \cdot x - k \cdot z\right) - \left(y \cdot x - z \cdot t\right) \cdot \left(b \cdot a - i \cdot c\right)\right)\right)\\ \mathbf{elif}\;y4 \lt 6.718963124057495 \cdot 10^{-279}:\\ \;\;\;\;\left(\left(\left(\left(k \cdot y\right) \cdot \left(y5 \cdot i\right) - \left(y \cdot b\right) \cdot \left(y4 \cdot k\right)\right) - \left(y5 \cdot t\right) \cdot \left(i \cdot j\right)\right) - \left(\left(y2 \cdot t - y3 \cdot y\right) \cdot \left(y4 \cdot c - y5 \cdot a\right) - \left(y2 \cdot k - y3 \cdot j\right) \cdot \left(y4 \cdot y1 - y5 \cdot y0\right)\right)\right) + \left(\left(b \cdot a - i \cdot c\right) \cdot \left(y \cdot x - t \cdot z\right) - \left(\left(j \cdot x - k \cdot z\right) \cdot \left(y0 \cdot b - i \cdot y1\right) - \left(y2 \cdot x - y3 \cdot z\right) \cdot \left(c \cdot y0 - y1 \cdot a\right)\right)\right)\\ \mathbf{elif}\;y4 \lt 4.77962681403792 \cdot 10^{-222}:\\ \;\;\;\;\left(\left(\left(j \cdot t - k \cdot y\right) \cdot \left(y4 \cdot b - y5 \cdot i\right) - \left(y3 \cdot y\right) \cdot \left(y5 \cdot a - y4 \cdot c\right)\right) + \left(\left(y5 \cdot a\right) \cdot \left(t \cdot y2\right) + \left(k \cdot y2 - j \cdot y3\right) \cdot \left(y4 \cdot y1 - y5 \cdot y0\right)\right)\right) + \left(\left(x \cdot y2 - z \cdot y3\right) \cdot \left(c \cdot y0 - a \cdot y1\right) - \left(\left(b \cdot y0 - i \cdot y1\right) \cdot \left(j \cdot x - k \cdot z\right) - \left(y \cdot x - z \cdot t\right) \cdot \left(b \cdot a - i \cdot c\right)\right)\right)\\ \mathbf{elif}\;y4 \lt 2.2852241541266835 \cdot 10^{-175}:\\ \;\;\;\;\left(\left(\left(\left(k \cdot y\right) \cdot \left(y5 \cdot i\right) - \left(y \cdot b\right) \cdot \left(y4 \cdot k\right)\right) - \left(y5 \cdot t\right) \cdot \left(i \cdot j\right)\right) - \left(\left(y2 \cdot t - y3 \cdot y\right) \cdot \left(y4 \cdot c - y5 \cdot a\right) - \left(y2 \cdot k - y3 \cdot j\right) \cdot \left(y4 \cdot y1 - y5 \cdot y0\right)\right)\right) + \left(\left(b \cdot a - i \cdot c\right) \cdot \left(y \cdot x - t \cdot z\right) - \left(\left(j \cdot x - k \cdot z\right) \cdot \left(y0 \cdot b - i \cdot y1\right) - \left(y2 \cdot x - y3 \cdot z\right) \cdot \left(c \cdot y0 - y1 \cdot a\right)\right)\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(k \cdot \left(i \cdot \left(z \cdot y1\right)\right) - \left(j \cdot \left(i \cdot \left(x \cdot y1\right)\right) + y0 \cdot \left(k \cdot \left(z \cdot b\right)\right)\right)\right)\right) + \left(z \cdot \left(y3 \cdot \left(a \cdot y1\right)\right) - \left(y2 \cdot \left(x \cdot \left(a \cdot y1\right)\right) + y0 \cdot \left(z \cdot \left(c \cdot y3\right)\right)\right)\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)\\ \end{array}\]

Derivation

1. Split input into 3 regimes

2. if y2 < -1.2961541745563326e+62

   1. Initial program 27.8

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

      \[\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}{0}\]

   if -1.2961541745563326e+62 < y2 < 2.707765602532912e-42

   1. Initial program 24.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 25.9

      \[\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(a \cdot \left(y3 \cdot \left(y1 \cdot z\right)\right) - \left(a \cdot \left(x \cdot \left(y2 \cdot y1\right)\right) + c \cdot \left(z \cdot \left(y3 \cdot y0\right)\right)\right)\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 25.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) + \color{blue}{\left(a \cdot \left(y3 \cdot \left(y1 \cdot z\right) - x \cdot \left(y1 \cdot y2\right)\right) - \left(\left(z \cdot c\right) \cdot y3\right) \cdot y0\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 2.707765602532912e-42 < y2

   1. Initial program 26.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 27.7

      \[\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. Recombined 3 regimes into one program.

4. Final simplification 27.0

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

Reproduce

herbie shell --seed 2019162 
(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"

  :herbie-target
  (if (< y4 -7.206256231996481e+60) (- (- (* (- (* b a) (* i c)) (- (* y x) (* t z))) (- (* (- (* j x) (* k z)) (- (* y0 b) (* i y1))) (* (- (* j t) (* k y)) (- (* y4 b) (* y5 i))))) (- (/ (- (* y2 t) (* y3 y)) (/ 1 (- (* y4 c) (* y5 a)))) (* (- (* y2 k) (* y3 j)) (- (* y4 y1) (* y5 y0))))) (if (< y4 -3.364603505246317e-66) (+ (- (- (- (* (* t c) (* i z)) (* (* a t) (* b z))) (* (* y c) (* i x))) (* (- (* b y0) (* i y1)) (- (* j x) (* k z)))) (- (* (- (* y0 c) (* a y1)) (- (* x y2) (* z y3))) (- (* (- (* t y2) (* y y3)) (- (* y4 c) (* a y5))) (* (- (* y1 y4) (* y5 y0)) (- (* k y2) (* j y3)))))) (if (< y4 -1.2000065055686116e-105) (+ (+ (- (* (- (* j t) (* k y)) (- (* y4 b) (* y5 i))) (* (* y3 y) (- (* y5 a) (* y4 c)))) (+ (* (* y5 a) (* t y2)) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* x y2) (* z y3)) (- (* c y0) (* a y1))) (- (* (- (* b y0) (* i y1)) (- (* j x) (* k z))) (* (- (* y x) (* z t)) (- (* b a) (* i c)))))) (if (< y4 6.718963124057495e-279) (+ (- (- (- (* (* k y) (* y5 i)) (* (* y b) (* y4 k))) (* (* y5 t) (* i j))) (- (* (- (* y2 t) (* y3 y)) (- (* y4 c) (* y5 a))) (* (- (* y2 k) (* y3 j)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* b a) (* i c)) (- (* y x) (* t z))) (- (* (- (* j x) (* k z)) (- (* y0 b) (* i y1))) (* (- (* y2 x) (* y3 z)) (- (* c y0) (* y1 a)))))) (if (< y4 4.77962681403792e-222) (+ (+ (- (* (- (* j t) (* k y)) (- (* y4 b) (* y5 i))) (* (* y3 y) (- (* y5 a) (* y4 c)))) (+ (* (* y5 a) (* t y2)) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* x y2) (* z y3)) (- (* c y0) (* a y1))) (- (* (- (* b y0) (* i y1)) (- (* j x) (* k z))) (* (- (* y x) (* z t)) (- (* b a) (* i c)))))) (if (< y4 2.2852241541266835e-175) (+ (- (- (- (* (* k y) (* y5 i)) (* (* y b) (* y4 k))) (* (* y5 t) (* i j))) (- (* (- (* y2 t) (* y3 y)) (- (* y4 c) (* y5 a))) (* (- (* y2 k) (* y3 j)) (- (* y4 y1) (* y5 y0))))) (- (* (- (* b a) (* i c)) (- (* y x) (* t z))) (- (* (- (* j x) (* k z)) (- (* y0 b) (* i y1))) (* (- (* y2 x) (* y3 z)) (- (* c y0) (* y1 a)))))) (+ (- (+ (+ (- (* (- (* x y) (* z t)) (- (* a b) (* c i))) (- (* k (* i (* z y1))) (+ (* j (* i (* x y1))) (* y0 (* k (* z b)))))) (- (* z (* y3 (* a y1))) (+ (* y2 (* x (* a y1))) (* y0 (* z (* c y3)))))) (* (- (* t j) (* y k)) (- (* y4 b) (* y5 i)))) (* (- (* t y2) (* y y3)) (- (* y4 c) (* y5 a)))) (* (- (* k y2) (* j y3)) (- (* y4 y1) (* y5 y0))))))))))

  (+ (- (+ (+ (- (* (- (* 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)))))