Linear.Matrix:det33 from linear-1.19.1.3

Average Error: 12.2 → 11.8

Time: 24.7s

Precision: 64

Math

\[\left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \left(c \cdot t - i \cdot y\right)\]

↓

\[\begin{array}{l} \mathbf{if}\;j \le -3.778273061350026247693292264821399017336 \cdot 10^{98}:\\ \;\;\;\;j \cdot \left(c \cdot t - i \cdot y\right) + \left(-b \cdot \left(c \cdot z - i \cdot a\right)\right)\\ \mathbf{elif}\;j \le -3.030521279964126763770820334228788464675 \cdot 10^{-144}:\\ \;\;\;\;\left(\left(\left(y \cdot z\right) \cdot x + \left(-a \cdot \left(x \cdot t\right)\right)\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(t \cdot \left(j \cdot c\right) + \left(-i \cdot \left(j \cdot y\right)\right)\right)\\ \mathbf{elif}\;j \le -2.084139175117578035405450843490766739903 \cdot 10^{-260}:\\ \;\;\;\;\left(x \cdot \left(y \cdot z - t \cdot a\right) - \left(z \cdot \left(b \cdot c\right) + \left(-i \cdot a\right) \cdot b\right)\right) + \left(\left(t \cdot j\right) \cdot c + \left(-i \cdot \left(j \cdot y\right)\right)\right)\\ \mathbf{elif}\;j \le 2.649066615389731199752472489106042554466 \cdot 10^{253}:\\ \;\;\;\;\left(x \cdot \left(y \cdot z - t \cdot a\right) - \left(b \cdot \left(\sqrt[3]{c \cdot z - i \cdot a} \cdot \sqrt[3]{c \cdot z - i \cdot a}\right)\right) \cdot \sqrt[3]{c \cdot z - i \cdot a}\right) + \left(\left(t \cdot j\right) \cdot c + \left(-i \cdot \left(j \cdot y\right)\right)\right)\\ \mathbf{else}:\\ \;\;\;\;j \cdot \left(c \cdot t - i \cdot y\right) + \left(-b \cdot \left(c \cdot z - i \cdot a\right)\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 r81572 = x;
        double r81573 = y;
        double r81574 = z;
        double r81575 = r81573 * r81574;
        double r81576 = t;
        double r81577 = a;
        double r81578 = r81576 * r81577;
        double r81579 = r81575 - r81578;
        double r81580 = r81572 * r81579;
        double r81581 = b;
        double r81582 = c;
        double r81583 = r81582 * r81574;
        double r81584 = i;
        double r81585 = r81584 * r81577;
        double r81586 = r81583 - r81585;
        double r81587 = r81581 * r81586;
        double r81588 = r81580 - r81587;
        double r81589 = j;
        double r81590 = r81582 * r81576;
        double r81591 = r81584 * r81573;
        double r81592 = r81590 - r81591;
        double r81593 = r81589 * r81592;
        double r81594 = r81588 + r81593;
        return r81594;
}

↓

double f(double x, double y, double z, double t, double a, double b, double c, double i, double j) {
        double r81595 = j;
        double r81596 = -3.778273061350026e+98;
        bool r81597 = r81595 <= r81596;
        double r81598 = c;
        double r81599 = t;
        double r81600 = r81598 * r81599;
        double r81601 = i;
        double r81602 = y;
        double r81603 = r81601 * r81602;
        double r81604 = r81600 - r81603;
        double r81605 = r81595 * r81604;
        double r81606 = b;
        double r81607 = z;
        double r81608 = r81598 * r81607;
        double r81609 = a;
        double r81610 = r81601 * r81609;
        double r81611 = r81608 - r81610;
        double r81612 = r81606 * r81611;
        double r81613 = -r81612;
        double r81614 = r81605 + r81613;
        double r81615 = -3.030521279964127e-144;
        bool r81616 = r81595 <= r81615;
        double r81617 = r81602 * r81607;
        double r81618 = x;
        double r81619 = r81617 * r81618;
        double r81620 = r81618 * r81599;
        double r81621 = r81609 * r81620;
        double r81622 = -r81621;
        double r81623 = r81619 + r81622;
        double r81624 = r81623 - r81612;
        double r81625 = r81595 * r81598;
        double r81626 = r81599 * r81625;
        double r81627 = r81595 * r81602;
        double r81628 = r81601 * r81627;
        double r81629 = -r81628;
        double r81630 = r81626 + r81629;
        double r81631 = r81624 + r81630;
        double r81632 = -2.084139175117578e-260;
        bool r81633 = r81595 <= r81632;
        double r81634 = r81599 * r81609;
        double r81635 = r81617 - r81634;
        double r81636 = r81618 * r81635;
        double r81637 = r81606 * r81598;
        double r81638 = r81607 * r81637;
        double r81639 = -r81610;
        double r81640 = r81639 * r81606;
        double r81641 = r81638 + r81640;
        double r81642 = r81636 - r81641;
        double r81643 = r81599 * r81595;
        double r81644 = r81643 * r81598;
        double r81645 = r81644 + r81629;
        double r81646 = r81642 + r81645;
        double r81647 = 2.6490666153897312e+253;
        bool r81648 = r81595 <= r81647;
        double r81649 = cbrt(r81611);
        double r81650 = r81649 * r81649;
        double r81651 = r81606 * r81650;
        double r81652 = r81651 * r81649;
        double r81653 = r81636 - r81652;
        double r81654 = r81653 + r81645;
        double r81655 = r81648 ? r81654 : r81614;
        double r81656 = r81633 ? r81646 : r81655;
        double r81657 = r81616 ? r81631 : r81656;
        double r81658 = r81597 ? r81614 : r81657;
        return r81658;
}

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

Derivation

1. Split input into 4 regimes

2. if j < -3.778273061350026e+98 or 2.6490666153897312e+253 < j

   1. Initial program 7.2

      \[\left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \left(c \cdot t - i \cdot y\right)\]

   2. Taylor expanded around 0 15.7

      \[\leadsto \left(\color{blue}{0} - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \left(c \cdot t - i \cdot y\right)\]

   if -3.778273061350026e+98 < j < -3.030521279964127e-144

   1. Initial program 11.1

      \[\left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \left(c \cdot t - i \cdot y\right)\]
   2. Using strategy rm

   3. Applied sub-neg 11.1

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \color{blue}{\left(c \cdot t + \left(-i \cdot y\right)\right)}\]

   4. Applied distribute-lft-in 11.1

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \color{blue}{\left(j \cdot \left(c \cdot t\right) + j \cdot \left(-i \cdot y\right)\right)}\]

   5. Simplified 9.9

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(\color{blue}{t \cdot \left(j \cdot c\right)} + j \cdot \left(-i \cdot y\right)\right)\]

   6. Simplified 9.0

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(t \cdot \left(j \cdot c\right) + \color{blue}{\left(-i \cdot \left(j \cdot y\right)\right)}\right)\]
   7. Using strategy rm

   8. Applied sub-neg 9.0

      \[\leadsto \left(x \cdot \color{blue}{\left(y \cdot z + \left(-t \cdot a\right)\right)} - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(t \cdot \left(j \cdot c\right) + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]

   9. Applied distribute-lft-in 9.0

      \[\leadsto \left(\color{blue}{\left(x \cdot \left(y \cdot z\right) + x \cdot \left(-t \cdot a\right)\right)} - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(t \cdot \left(j \cdot c\right) + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]

   10. Simplified 9.0

       \[\leadsto \left(\left(\color{blue}{\left(y \cdot z\right) \cdot x} + x \cdot \left(-t \cdot a\right)\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(t \cdot \left(j \cdot c\right) + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]

   11. Simplified 9.4

       \[\leadsto \left(\left(\left(y \cdot z\right) \cdot x + \color{blue}{\left(-a \cdot \left(x \cdot t\right)\right)}\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(t \cdot \left(j \cdot c\right) + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]

   if -3.030521279964127e-144 < j < -2.084139175117578e-260

   1. Initial program 17.5

      \[\left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \left(c \cdot t - i \cdot y\right)\]
   2. Using strategy rm

   3. Applied sub-neg 17.5

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \color{blue}{\left(c \cdot t + \left(-i \cdot y\right)\right)}\]

   4. Applied distribute-lft-in 17.5

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \color{blue}{\left(j \cdot \left(c \cdot t\right) + j \cdot \left(-i \cdot y\right)\right)}\]

   5. Simplified 14.8

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(\color{blue}{t \cdot \left(j \cdot c\right)} + j \cdot \left(-i \cdot y\right)\right)\]

   6. Simplified 11.3

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(t \cdot \left(j \cdot c\right) + \color{blue}{\left(-i \cdot \left(j \cdot y\right)\right)}\right)\]
   7. Using strategy rm

   8. Applied associate-*r* 11.3

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(\color{blue}{\left(t \cdot j\right) \cdot c} + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]
   9. Using strategy rm

   10. Applied sub-neg 11.3

       \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \color{blue}{\left(c \cdot z + \left(-i \cdot a\right)\right)}\right) + \left(\left(t \cdot j\right) \cdot c + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]

   11. Applied distribute-lft-in 11.3

       \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - \color{blue}{\left(b \cdot \left(c \cdot z\right) + b \cdot \left(-i \cdot a\right)\right)}\right) + \left(\left(t \cdot j\right) \cdot c + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]

   12. Simplified 11.6

       \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - \left(\color{blue}{z \cdot \left(b \cdot c\right)} + b \cdot \left(-i \cdot a\right)\right)\right) + \left(\left(t \cdot j\right) \cdot c + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]

   13. Simplified 11.6

       \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - \left(z \cdot \left(b \cdot c\right) + \color{blue}{\left(-i \cdot a\right) \cdot b}\right)\right) + \left(\left(t \cdot j\right) \cdot c + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]

   if -2.084139175117578e-260 < j < 2.6490666153897312e+253

   1. Initial program 12.9

      \[\left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \left(c \cdot t - i \cdot y\right)\]
   2. Using strategy rm

   3. Applied sub-neg 12.9

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + j \cdot \color{blue}{\left(c \cdot t + \left(-i \cdot y\right)\right)}\]

   4. Applied distribute-lft-in 12.9

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \color{blue}{\left(j \cdot \left(c \cdot t\right) + j \cdot \left(-i \cdot y\right)\right)}\]

   5. Simplified 12.6

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(\color{blue}{t \cdot \left(j \cdot c\right)} + j \cdot \left(-i \cdot y\right)\right)\]

   6. Simplified 11.7

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(t \cdot \left(j \cdot c\right) + \color{blue}{\left(-i \cdot \left(j \cdot y\right)\right)}\right)\]
   7. Using strategy rm

   8. Applied associate-*r* 11.5

      \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(\color{blue}{\left(t \cdot j\right) \cdot c} + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]
   9. Using strategy rm

   10. Applied add-cube-cbrt 11.9

       \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - b \cdot \color{blue}{\left(\left(\sqrt[3]{c \cdot z - i \cdot a} \cdot \sqrt[3]{c \cdot z - i \cdot a}\right) \cdot \sqrt[3]{c \cdot z - i \cdot a}\right)}\right) + \left(\left(t \cdot j\right) \cdot c + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]

   11. Applied associate-*r* 11.9

       \[\leadsto \left(x \cdot \left(y \cdot z - t \cdot a\right) - \color{blue}{\left(b \cdot \left(\sqrt[3]{c \cdot z - i \cdot a} \cdot \sqrt[3]{c \cdot z - i \cdot a}\right)\right) \cdot \sqrt[3]{c \cdot z - i \cdot a}}\right) + \left(\left(t \cdot j\right) \cdot c + \left(-i \cdot \left(j \cdot y\right)\right)\right)\]
3. Recombined 4 regimes into one program.

4. Final simplification 11.8

   \[\leadsto \begin{array}{l} \mathbf{if}\;j \le -3.778273061350026247693292264821399017336 \cdot 10^{98}:\\ \;\;\;\;j \cdot \left(c \cdot t - i \cdot y\right) + \left(-b \cdot \left(c \cdot z - i \cdot a\right)\right)\\ \mathbf{elif}\;j \le -3.030521279964126763770820334228788464675 \cdot 10^{-144}:\\ \;\;\;\;\left(\left(\left(y \cdot z\right) \cdot x + \left(-a \cdot \left(x \cdot t\right)\right)\right) - b \cdot \left(c \cdot z - i \cdot a\right)\right) + \left(t \cdot \left(j \cdot c\right) + \left(-i \cdot \left(j \cdot y\right)\right)\right)\\ \mathbf{elif}\;j \le -2.084139175117578035405450843490766739903 \cdot 10^{-260}:\\ \;\;\;\;\left(x \cdot \left(y \cdot z - t \cdot a\right) - \left(z \cdot \left(b \cdot c\right) + \left(-i \cdot a\right) \cdot b\right)\right) + \left(\left(t \cdot j\right) \cdot c + \left(-i \cdot \left(j \cdot y\right)\right)\right)\\ \mathbf{elif}\;j \le 2.649066615389731199752472489106042554466 \cdot 10^{253}:\\ \;\;\;\;\left(x \cdot \left(y \cdot z - t \cdot a\right) - \left(b \cdot \left(\sqrt[3]{c \cdot z - i \cdot a} \cdot \sqrt[3]{c \cdot z - i \cdot a}\right)\right) \cdot \sqrt[3]{c \cdot z - i \cdot a}\right) + \left(\left(t \cdot j\right) \cdot c + \left(-i \cdot \left(j \cdot y\right)\right)\right)\\ \mathbf{else}:\\ \;\;\;\;j \cdot \left(c \cdot t - i \cdot y\right) + \left(-b \cdot \left(c \cdot z - i \cdot a\right)\right)\\ \end{array}\]

Reproduce

herbie shell --seed 2019325 +o rules:numerics
(FPCore (x y z t a b c i j)
  :name "Linear.Matrix:det33 from linear-1.19.1.3"
  :precision binary64
  (+ (- (* x (- (* y z) (* t a))) (* b (- (* c z) (* i a)))) (* j (- (* c t) (* i y)))))