\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}\;z \le 1.21671009256391355 \cdot 10^{-285}:\\
\;\;\;\;\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(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)\\
\mathbf{elif}\;z \le 6.02979551995023046 \cdot 10^{-157}:\\
\;\;\;\;\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(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}\;z \le 1.0718132790979475 \cdot 10^{-146}:\\
\;\;\;\;\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}\;z \le 5.0097592564704197 \cdot 10^{-52}:\\
\;\;\;\;\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}\;z \le 1.3684848252489064 \cdot 10^{211}:\\
\;\;\;\;\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(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)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(\left(\left(\left(x \cdot y - z \cdot t\right) \cdot \left(a \cdot b - c \cdot i\right) - 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(\left(t \cdot y2 - y \cdot y3\right) \cdot \left(\sqrt[3]{y4 \cdot c - y5 \cdot a} \cdot \sqrt[3]{y4 \cdot c - y5 \cdot a}\right)\right) \cdot \sqrt[3]{y4 \cdot c - y5 \cdot a}\right) + \left(k \cdot y2 - j \cdot y3\right) \cdot \left(y4 \cdot y1 - y5 \cdot y0\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 r145203 = x;
double r145204 = y;
double r145205 = r145203 * r145204;
double r145206 = z;
double r145207 = t;
double r145208 = r145206 * r145207;
double r145209 = r145205 - r145208;
double r145210 = a;
double r145211 = b;
double r145212 = r145210 * r145211;
double r145213 = c;
double r145214 = i;
double r145215 = r145213 * r145214;
double r145216 = r145212 - r145215;
double r145217 = r145209 * r145216;
double r145218 = j;
double r145219 = r145203 * r145218;
double r145220 = k;
double r145221 = r145206 * r145220;
double r145222 = r145219 - r145221;
double r145223 = y0;
double r145224 = r145223 * r145211;
double r145225 = y1;
double r145226 = r145225 * r145214;
double r145227 = r145224 - r145226;
double r145228 = r145222 * r145227;
double r145229 = r145217 - r145228;
double r145230 = y2;
double r145231 = r145203 * r145230;
double r145232 = y3;
double r145233 = r145206 * r145232;
double r145234 = r145231 - r145233;
double r145235 = r145223 * r145213;
double r145236 = r145225 * r145210;
double r145237 = r145235 - r145236;
double r145238 = r145234 * r145237;
double r145239 = r145229 + r145238;
double r145240 = r145207 * r145218;
double r145241 = r145204 * r145220;
double r145242 = r145240 - r145241;
double r145243 = y4;
double r145244 = r145243 * r145211;
double r145245 = y5;
double r145246 = r145245 * r145214;
double r145247 = r145244 - r145246;
double r145248 = r145242 * r145247;
double r145249 = r145239 + r145248;
double r145250 = r145207 * r145230;
double r145251 = r145204 * r145232;
double r145252 = r145250 - r145251;
double r145253 = r145243 * r145213;
double r145254 = r145245 * r145210;
double r145255 = r145253 - r145254;
double r145256 = r145252 * r145255;
double r145257 = r145249 - r145256;
double r145258 = r145220 * r145230;
double r145259 = r145218 * r145232;
double r145260 = r145258 - r145259;
double r145261 = r145243 * r145225;
double r145262 = r145245 * r145223;
double r145263 = r145261 - r145262;
double r145264 = r145260 * r145263;
double r145265 = r145257 + r145264;
return r145265;
}
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 r145266 = z;
double r145267 = 1.2167100925639136e-285;
bool r145268 = r145266 <= r145267;
double r145269 = x;
double r145270 = y;
double r145271 = r145269 * r145270;
double r145272 = t;
double r145273 = r145266 * r145272;
double r145274 = r145271 - r145273;
double r145275 = a;
double r145276 = b;
double r145277 = r145275 * r145276;
double r145278 = c;
double r145279 = i;
double r145280 = r145278 * r145279;
double r145281 = r145277 - r145280;
double r145282 = r145274 * r145281;
double r145283 = j;
double r145284 = r145269 * r145283;
double r145285 = k;
double r145286 = r145266 * r145285;
double r145287 = r145284 - r145286;
double r145288 = y0;
double r145289 = r145288 * r145276;
double r145290 = y1;
double r145291 = r145290 * r145279;
double r145292 = r145289 - r145291;
double r145293 = r145287 * r145292;
double r145294 = r145282 - r145293;
double r145295 = y2;
double r145296 = r145269 * r145295;
double r145297 = y3;
double r145298 = r145266 * r145297;
double r145299 = r145296 - r145298;
double r145300 = r145288 * r145278;
double r145301 = r145290 * r145275;
double r145302 = r145300 - r145301;
double r145303 = r145299 * r145302;
double r145304 = r145294 + r145303;
double r145305 = r145272 * r145283;
double r145306 = r145270 * r145285;
double r145307 = r145305 - r145306;
double r145308 = y4;
double r145309 = r145308 * r145276;
double r145310 = y5;
double r145311 = r145310 * r145279;
double r145312 = r145309 - r145311;
double r145313 = r145307 * r145312;
double r145314 = r145304 + r145313;
double r145315 = r145272 * r145295;
double r145316 = r145270 * r145297;
double r145317 = r145315 - r145316;
double r145318 = r145308 * r145278;
double r145319 = r145310 * r145275;
double r145320 = r145318 - r145319;
double r145321 = r145317 * r145320;
double r145322 = r145314 - r145321;
double r145323 = r145283 * r145310;
double r145324 = r145297 * r145323;
double r145325 = r145288 * r145324;
double r145326 = r145285 * r145310;
double r145327 = r145295 * r145326;
double r145328 = r145288 * r145327;
double r145329 = r145283 * r145308;
double r145330 = r145297 * r145329;
double r145331 = r145290 * r145330;
double r145332 = r145328 + r145331;
double r145333 = r145325 - r145332;
double r145334 = r145322 + r145333;
double r145335 = 6.02979551995023e-157;
bool r145336 = r145266 <= r145335;
double r145337 = r145290 * r145266;
double r145338 = r145297 * r145337;
double r145339 = r145275 * r145338;
double r145340 = r145297 * r145278;
double r145341 = r145266 * r145340;
double r145342 = r145288 * r145341;
double r145343 = r145295 * r145290;
double r145344 = r145269 * r145343;
double r145345 = r145275 * r145344;
double r145346 = r145342 + r145345;
double r145347 = r145339 - r145346;
double r145348 = r145294 + r145347;
double r145349 = r145348 + r145313;
double r145350 = r145349 - r145321;
double r145351 = r145285 * r145295;
double r145352 = r145283 * r145297;
double r145353 = r145351 - r145352;
double r145354 = r145308 * r145290;
double r145355 = r145310 * r145288;
double r145356 = r145354 - r145355;
double r145357 = r145353 * r145356;
double r145358 = r145350 + r145357;
double r145359 = 1.0718132790979475e-146;
bool r145360 = r145266 <= r145359;
double r145361 = r145266 * r145278;
double r145362 = r145279 * r145361;
double r145363 = r145272 * r145362;
double r145364 = r145270 * r145269;
double r145365 = r145278 * r145364;
double r145366 = r145279 * r145365;
double r145367 = r145266 * r145276;
double r145368 = r145272 * r145367;
double r145369 = r145275 * r145368;
double r145370 = r145366 + r145369;
double r145371 = r145363 - r145370;
double r145372 = r145371 - r145293;
double r145373 = r145372 + r145303;
double r145374 = r145373 + r145313;
double r145375 = r145374 - r145321;
double r145376 = r145375 + r145357;
double r145377 = 5.00975925647042e-52;
bool r145378 = r145266 <= r145377;
double r145379 = 0.0;
double r145380 = r145314 - r145379;
double r145381 = r145380 + r145357;
double r145382 = 1.3684848252489064e+211;
bool r145383 = r145266 <= r145382;
double r145384 = r145282 - r145379;
double r145385 = r145384 + r145303;
double r145386 = r145385 + r145313;
double r145387 = cbrt(r145320);
double r145388 = r145387 * r145387;
double r145389 = r145317 * r145388;
double r145390 = r145389 * r145387;
double r145391 = r145386 - r145390;
double r145392 = r145391 + r145357;
double r145393 = r145383 ? r145334 : r145392;
double r145394 = r145378 ? r145381 : r145393;
double r145395 = r145360 ? r145376 : r145394;
double r145396 = r145336 ? r145358 : r145395;
double r145397 = r145268 ? r145334 : r145396;
return r145397;
}



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 z < 1.2167100925639136e-285 or 5.00975925647042e-52 < z < 1.3684848252489064e+211Initial program 26.9
Taylor expanded around inf 29.8
if 1.2167100925639136e-285 < z < 6.02979551995023e-157Initial program 28.0
Taylor expanded around inf 30.7
if 6.02979551995023e-157 < z < 1.0718132790979475e-146Initial program 26.3
Taylor expanded around inf 26.2
if 1.0718132790979475e-146 < z < 5.00975925647042e-52Initial program 24.4
Taylor expanded around 0 30.1
if 1.3684848252489064e+211 < z Initial program 35.5
rmApplied add-cube-cbrt35.5
Applied associate-*r*35.5
Taylor expanded around 0 38.3
Final simplification30.2
herbie shell --seed 2020049
(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)))))