Linear.Matrix:det44 from linear-1.19.1.3

Average Error: 25.2 → 27.0

Time: 2.0m

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 r27507637 = x;
        double r27507638 = y;
        double r27507639 = r27507637 * r27507638;
        double r27507640 = z;
        double r27507641 = t;
        double r27507642 = r27507640 * r27507641;
        double r27507643 = r27507639 - r27507642;
        double r27507644 = a;
        double r27507645 = b;
        double r27507646 = r27507644 * r27507645;
        double r27507647 = c;
        double r27507648 = i;
        double r27507649 = r27507647 * r27507648;
        double r27507650 = r27507646 - r27507649;
        double r27507651 = r27507643 * r27507650;
        double r27507652 = j;
        double r27507653 = r27507637 * r27507652;
        double r27507654 = k;
        double r27507655 = r27507640 * r27507654;
        double r27507656 = r27507653 - r27507655;
        double r27507657 = y0;
        double r27507658 = r27507657 * r27507645;
        double r27507659 = y1;
        double r27507660 = r27507659 * r27507648;
        double r27507661 = r27507658 - r27507660;
        double r27507662 = r27507656 * r27507661;
        double r27507663 = r27507651 - r27507662;
        double r27507664 = y2;
        double r27507665 = r27507637 * r27507664;
        double r27507666 = y3;
        double r27507667 = r27507640 * r27507666;
        double r27507668 = r27507665 - r27507667;
        double r27507669 = r27507657 * r27507647;
        double r27507670 = r27507659 * r27507644;
        double r27507671 = r27507669 - r27507670;
        double r27507672 = r27507668 * r27507671;
        double r27507673 = r27507663 + r27507672;
        double r27507674 = r27507641 * r27507652;
        double r27507675 = r27507638 * r27507654;
        double r27507676 = r27507674 - r27507675;
        double r27507677 = y4;
        double r27507678 = r27507677 * r27507645;
        double r27507679 = y5;
        double r27507680 = r27507679 * r27507648;
        double r27507681 = r27507678 - r27507680;
        double r27507682 = r27507676 * r27507681;
        double r27507683 = r27507673 + r27507682;
        double r27507684 = r27507641 * r27507664;
        double r27507685 = r27507638 * r27507666;
        double r27507686 = r27507684 - r27507685;
        double r27507687 = r27507677 * r27507647;
        double r27507688 = r27507679 * r27507644;
        double r27507689 = r27507687 - r27507688;
        double r27507690 = r27507686 * r27507689;
        double r27507691 = r27507683 - r27507690;
        double r27507692 = r27507654 * r27507664;
        double r27507693 = r27507652 * r27507666;
        double r27507694 = r27507692 - r27507693;
        double r27507695 = r27507677 * r27507659;
        double r27507696 = r27507679 * r27507657;
        double r27507697 = r27507695 - r27507696;
        double r27507698 = r27507694 * r27507697;
        double r27507699 = r27507691 + r27507698;
        return r27507699;
}

↓

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 r27507700 = y2;
        double r27507701 = -1.2961541745563326e+62;
        bool r27507702 = r27507700 <= r27507701;
        double r27507703 = y0;
        double r27507704 = c;
        double r27507705 = r27507703 * r27507704;
        double r27507706 = y1;
        double r27507707 = a;
        double r27507708 = r27507706 * r27507707;
        double r27507709 = r27507705 - r27507708;
        double r27507710 = x;
        double r27507711 = r27507710 * r27507700;
        double r27507712 = z;
        double r27507713 = y3;
        double r27507714 = r27507712 * r27507713;
        double r27507715 = r27507711 - r27507714;
        double r27507716 = r27507709 * r27507715;
        double r27507717 = y;
        double r27507718 = r27507710 * r27507717;
        double r27507719 = t;
        double r27507720 = r27507719 * r27507712;
        double r27507721 = r27507718 - r27507720;
        double r27507722 = b;
        double r27507723 = r27507707 * r27507722;
        double r27507724 = i;
        double r27507725 = r27507724 * r27507704;
        double r27507726 = r27507723 - r27507725;
        double r27507727 = r27507721 * r27507726;
        double r27507728 = r27507703 * r27507722;
        double r27507729 = r27507724 * r27507706;
        double r27507730 = r27507728 - r27507729;
        double r27507731 = j;
        double r27507732 = r27507710 * r27507731;
        double r27507733 = k;
        double r27507734 = r27507712 * r27507733;
        double r27507735 = r27507732 - r27507734;
        double r27507736 = r27507730 * r27507735;
        double r27507737 = r27507727 - r27507736;
        double r27507738 = r27507716 + r27507737;
        double r27507739 = r27507719 * r27507731;
        double r27507740 = r27507717 * r27507733;
        double r27507741 = r27507739 - r27507740;
        double r27507742 = y4;
        double r27507743 = r27507742 * r27507722;
        double r27507744 = y5;
        double r27507745 = r27507724 * r27507744;
        double r27507746 = r27507743 - r27507745;
        double r27507747 = r27507741 * r27507746;
        double r27507748 = r27507738 + r27507747;
        double r27507749 = r27507719 * r27507700;
        double r27507750 = r27507717 * r27507713;
        double r27507751 = r27507749 - r27507750;
        double r27507752 = r27507742 * r27507704;
        double r27507753 = r27507744 * r27507707;
        double r27507754 = r27507752 - r27507753;
        double r27507755 = r27507751 * r27507754;
        double r27507756 = r27507748 - r27507755;
        double r27507757 = 2.707765602532912e-42;
        bool r27507758 = r27507700 <= r27507757;
        double r27507759 = r27507712 * r27507706;
        double r27507760 = r27507713 * r27507759;
        double r27507761 = r27507700 * r27507706;
        double r27507762 = r27507710 * r27507761;
        double r27507763 = r27507760 - r27507762;
        double r27507764 = r27507763 * r27507707;
        double r27507765 = r27507704 * r27507712;
        double r27507766 = r27507713 * r27507765;
        double r27507767 = r27507766 * r27507703;
        double r27507768 = r27507764 - r27507767;
        double r27507769 = r27507768 + r27507737;
        double r27507770 = r27507769 + r27507747;
        double r27507771 = r27507770 - r27507755;
        double r27507772 = r27507706 * r27507742;
        double r27507773 = r27507703 * r27507744;
        double r27507774 = r27507772 - r27507773;
        double r27507775 = r27507700 * r27507733;
        double r27507776 = r27507713 * r27507731;
        double r27507777 = r27507775 - r27507776;
        double r27507778 = r27507774 * r27507777;
        double r27507779 = r27507771 + r27507778;
        double r27507780 = r27507724 * r27507765;
        double r27507781 = r27507719 * r27507780;
        double r27507782 = r27507717 * r27507704;
        double r27507783 = r27507710 * r27507782;
        double r27507784 = r27507783 * r27507724;
        double r27507785 = r27507722 * r27507712;
        double r27507786 = r27507719 * r27507785;
        double r27507787 = r27507707 * r27507786;
        double r27507788 = r27507784 + r27507787;
        double r27507789 = r27507781 - r27507788;
        double r27507790 = r27507789 - r27507736;
        double r27507791 = r27507716 + r27507790;
        double r27507792 = r27507791 + r27507747;
        double r27507793 = r27507792 - r27507755;
        double r27507794 = r27507793 + r27507778;
        double r27507795 = r27507758 ? r27507779 : r27507794;
        double r27507796 = r27507702 ? r27507756 : r27507795;
        return r27507796;
}

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