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

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

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 r85581 = x;
        double r85582 = y;
        double r85583 = r85581 * r85582;
        double r85584 = z;
        double r85585 = t;
        double r85586 = r85584 * r85585;
        double r85587 = r85583 - r85586;
        double r85588 = a;
        double r85589 = b;
        double r85590 = r85588 * r85589;
        double r85591 = c;
        double r85592 = i;
        double r85593 = r85591 * r85592;
        double r85594 = r85590 - r85593;
        double r85595 = r85587 * r85594;
        double r85596 = j;
        double r85597 = r85581 * r85596;
        double r85598 = k;
        double r85599 = r85584 * r85598;
        double r85600 = r85597 - r85599;
        double r85601 = y0;
        double r85602 = r85601 * r85589;
        double r85603 = y1;
        double r85604 = r85603 * r85592;
        double r85605 = r85602 - r85604;
        double r85606 = r85600 * r85605;
        double r85607 = r85595 - r85606;
        double r85608 = y2;
        double r85609 = r85581 * r85608;
        double r85610 = y3;
        double r85611 = r85584 * r85610;
        double r85612 = r85609 - r85611;
        double r85613 = r85601 * r85591;
        double r85614 = r85603 * r85588;
        double r85615 = r85613 - r85614;
        double r85616 = r85612 * r85615;
        double r85617 = r85607 + r85616;
        double r85618 = r85585 * r85596;
        double r85619 = r85582 * r85598;
        double r85620 = r85618 - r85619;
        double r85621 = y4;
        double r85622 = r85621 * r85589;
        double r85623 = y5;
        double r85624 = r85623 * r85592;
        double r85625 = r85622 - r85624;
        double r85626 = r85620 * r85625;
        double r85627 = r85617 + r85626;
        double r85628 = r85585 * r85608;
        double r85629 = r85582 * r85610;
        double r85630 = r85628 - r85629;
        double r85631 = r85621 * r85591;
        double r85632 = r85623 * r85588;
        double r85633 = r85631 - r85632;
        double r85634 = r85630 * r85633;
        double r85635 = r85627 - r85634;
        double r85636 = r85598 * r85608;
        double r85637 = r85596 * r85610;
        double r85638 = r85636 - r85637;
        double r85639 = r85621 * r85603;
        double r85640 = r85623 * r85601;
        double r85641 = r85639 - r85640;
        double r85642 = r85638 * r85641;
        double r85643 = r85635 + r85642;
        return r85643;
}

Reproduce

herbie shell --seed 2020018 +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"
  :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)))))

Timeout in 10.0m

Reproduce