\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}\;y1 \le -7.867915292957589484241070210347083929061 \cdot 10^{50}:\\
\;\;\;\;\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}\;y1 \le -6.294234080972948076259974623937082195408 \cdot 10^{-158}:\\
\;\;\;\;\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(\left(x \cdot j - z \cdot k\right) \cdot \left(\sqrt[3]{y0 \cdot b - y1 \cdot i} \cdot \sqrt[3]{y0 \cdot b - y1 \cdot i}\right)\right) \cdot \sqrt[3]{y0 \cdot b - y1 \cdot i}\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}\;y1 \le -2.892383599226965167400921940423103999829 \cdot 10^{-251}:\\
\;\;\;\;\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}\;y1 \le 1.339934175676608121589714754571252408666 \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(a \cdot \left(y3 \cdot \left(y \cdot y5\right)\right) - \left(y \cdot \left(y3 \cdot \left(y4 \cdot c\right)\right) + y5 \cdot \left(a \cdot \left(y2 \cdot t\right)\right)\right)\right)\right) + \left(k \cdot y2 - j \cdot y3\right) \cdot \left(y4 \cdot y1 - y5 \cdot y0\right)\\
\mathbf{elif}\;y1 \le 1.720891527290092531582770702888956718479 \cdot 10^{-267}:\\
\;\;\;\;\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(\left(x \cdot j - z \cdot k\right) \cdot \left(\sqrt[3]{y0 \cdot b - y1 \cdot i} \cdot \sqrt[3]{y0 \cdot b - y1 \cdot i}\right)\right) \cdot \sqrt[3]{y0 \cdot b - y1 \cdot i}\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}\;y1 \le 5.949301058765557706195613572916436516488 \cdot 10^{53}:\\
\;\;\;\;\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(\sqrt[3]{\left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - y5 \cdot i\right)} \cdot \sqrt[3]{\left(t \cdot j - y \cdot k\right) \cdot \left(y4 \cdot b - y5 \cdot i\right)}\right) \cdot \sqrt[3]{\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{else}:\\
\;\;\;\;\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)\\
\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 r139762 = x;
double r139763 = y;
double r139764 = r139762 * r139763;
double r139765 = z;
double r139766 = t;
double r139767 = r139765 * r139766;
double r139768 = r139764 - r139767;
double r139769 = a;
double r139770 = b;
double r139771 = r139769 * r139770;
double r139772 = c;
double r139773 = i;
double r139774 = r139772 * r139773;
double r139775 = r139771 - r139774;
double r139776 = r139768 * r139775;
double r139777 = j;
double r139778 = r139762 * r139777;
double r139779 = k;
double r139780 = r139765 * r139779;
double r139781 = r139778 - r139780;
double r139782 = y0;
double r139783 = r139782 * r139770;
double r139784 = y1;
double r139785 = r139784 * r139773;
double r139786 = r139783 - r139785;
double r139787 = r139781 * r139786;
double r139788 = r139776 - r139787;
double r139789 = y2;
double r139790 = r139762 * r139789;
double r139791 = y3;
double r139792 = r139765 * r139791;
double r139793 = r139790 - r139792;
double r139794 = r139782 * r139772;
double r139795 = r139784 * r139769;
double r139796 = r139794 - r139795;
double r139797 = r139793 * r139796;
double r139798 = r139788 + r139797;
double r139799 = r139766 * r139777;
double r139800 = r139763 * r139779;
double r139801 = r139799 - r139800;
double r139802 = y4;
double r139803 = r139802 * r139770;
double r139804 = y5;
double r139805 = r139804 * r139773;
double r139806 = r139803 - r139805;
double r139807 = r139801 * r139806;
double r139808 = r139798 + r139807;
double r139809 = r139766 * r139789;
double r139810 = r139763 * r139791;
double r139811 = r139809 - r139810;
double r139812 = r139802 * r139772;
double r139813 = r139804 * r139769;
double r139814 = r139812 - r139813;
double r139815 = r139811 * r139814;
double r139816 = r139808 - r139815;
double r139817 = r139779 * r139789;
double r139818 = r139777 * r139791;
double r139819 = r139817 - r139818;
double r139820 = r139802 * r139784;
double r139821 = r139804 * r139782;
double r139822 = r139820 - r139821;
double r139823 = r139819 * r139822;
double r139824 = r139816 + r139823;
return r139824;
}
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 r139825 = y1;
double r139826 = -7.86791529295759e+50;
bool r139827 = r139825 <= r139826;
double r139828 = x;
double r139829 = y;
double r139830 = r139828 * r139829;
double r139831 = z;
double r139832 = t;
double r139833 = r139831 * r139832;
double r139834 = r139830 - r139833;
double r139835 = a;
double r139836 = b;
double r139837 = r139835 * r139836;
double r139838 = c;
double r139839 = i;
double r139840 = r139838 * r139839;
double r139841 = r139837 - r139840;
double r139842 = r139834 * r139841;
double r139843 = j;
double r139844 = r139828 * r139843;
double r139845 = k;
double r139846 = r139831 * r139845;
double r139847 = r139844 - r139846;
double r139848 = y0;
double r139849 = r139848 * r139836;
double r139850 = r139825 * r139839;
double r139851 = r139849 - r139850;
double r139852 = r139847 * r139851;
double r139853 = r139842 - r139852;
double r139854 = y2;
double r139855 = r139828 * r139854;
double r139856 = y3;
double r139857 = r139831 * r139856;
double r139858 = r139855 - r139857;
double r139859 = r139848 * r139838;
double r139860 = r139825 * r139835;
double r139861 = r139859 - r139860;
double r139862 = r139858 * r139861;
double r139863 = r139853 + r139862;
double r139864 = r139832 * r139843;
double r139865 = r139829 * r139845;
double r139866 = r139864 - r139865;
double r139867 = y4;
double r139868 = r139867 * r139836;
double r139869 = y5;
double r139870 = r139869 * r139839;
double r139871 = r139868 - r139870;
double r139872 = r139866 * r139871;
double r139873 = r139863 + r139872;
double r139874 = r139832 * r139854;
double r139875 = r139829 * r139856;
double r139876 = r139874 - r139875;
double r139877 = r139867 * r139838;
double r139878 = r139869 * r139835;
double r139879 = r139877 - r139878;
double r139880 = r139876 * r139879;
double r139881 = r139873 - r139880;
double r139882 = r139843 * r139869;
double r139883 = r139856 * r139882;
double r139884 = r139848 * r139883;
double r139885 = r139845 * r139869;
double r139886 = r139854 * r139885;
double r139887 = r139848 * r139886;
double r139888 = r139843 * r139867;
double r139889 = r139856 * r139888;
double r139890 = r139825 * r139889;
double r139891 = r139887 + r139890;
double r139892 = r139884 - r139891;
double r139893 = r139881 + r139892;
double r139894 = -6.294234080972948e-158;
bool r139895 = r139825 <= r139894;
double r139896 = r139831 * r139838;
double r139897 = r139839 * r139896;
double r139898 = r139832 * r139897;
double r139899 = r139829 * r139828;
double r139900 = r139838 * r139899;
double r139901 = r139839 * r139900;
double r139902 = r139831 * r139836;
double r139903 = r139832 * r139902;
double r139904 = r139835 * r139903;
double r139905 = r139901 + r139904;
double r139906 = r139898 - r139905;
double r139907 = cbrt(r139851);
double r139908 = r139907 * r139907;
double r139909 = r139847 * r139908;
double r139910 = r139909 * r139907;
double r139911 = r139906 - r139910;
double r139912 = r139911 + r139862;
double r139913 = r139912 + r139872;
double r139914 = r139913 - r139880;
double r139915 = r139845 * r139854;
double r139916 = r139843 * r139856;
double r139917 = r139915 - r139916;
double r139918 = r139867 * r139825;
double r139919 = r139869 * r139848;
double r139920 = r139918 - r139919;
double r139921 = r139917 * r139920;
double r139922 = r139914 + r139921;
double r139923 = -2.892383599226965e-251;
bool r139924 = r139825 <= r139923;
double r139925 = 1.3399341756766081e-285;
bool r139926 = r139825 <= r139925;
double r139927 = r139829 * r139869;
double r139928 = r139856 * r139927;
double r139929 = r139835 * r139928;
double r139930 = r139856 * r139877;
double r139931 = r139829 * r139930;
double r139932 = r139854 * r139832;
double r139933 = r139835 * r139932;
double r139934 = r139869 * r139933;
double r139935 = r139931 + r139934;
double r139936 = r139929 - r139935;
double r139937 = r139873 - r139936;
double r139938 = r139937 + r139921;
double r139939 = 1.7208915272900925e-267;
bool r139940 = r139825 <= r139939;
double r139941 = 5.949301058765558e+53;
bool r139942 = r139825 <= r139941;
double r139943 = cbrt(r139872);
double r139944 = r139943 * r139943;
double r139945 = r139944 * r139943;
double r139946 = r139863 + r139945;
double r139947 = r139946 - r139880;
double r139948 = r139947 + r139921;
double r139949 = r139942 ? r139948 : r139893;
double r139950 = r139940 ? r139922 : r139949;
double r139951 = r139926 ? r139938 : r139950;
double r139952 = r139924 ? r139893 : r139951;
double r139953 = r139895 ? r139922 : r139952;
double r139954 = r139827 ? r139893 : r139953;
return r139954;
}



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 y1 < -7.86791529295759e+50 or -6.294234080972948e-158 < y1 < -2.892383599226965e-251 or 5.949301058765558e+53 < y1 Initial program 29.0
Taylor expanded around inf 31.1
if -7.86791529295759e+50 < y1 < -6.294234080972948e-158 or 1.3399341756766081e-285 < y1 < 1.7208915272900925e-267Initial program 26.3
rmApplied add-cube-cbrt26.4
Applied associate-*r*26.4
Taylor expanded around inf 29.1
if -2.892383599226965e-251 < y1 < 1.3399341756766081e-285Initial program 27.2
Taylor expanded around inf 30.4
if 1.7208915272900925e-267 < y1 < 5.949301058765558e+53Initial program 24.7
rmApplied add-cube-cbrt24.8
Final simplification28.5
herbie shell --seed 2019298
(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)))))