\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}
t_1 := x \cdot y - z \cdot t\\
t_2 := t_1 \cdot \left(a \cdot b - c \cdot i\right)\\
t_3 := t \cdot y2 - y \cdot y3\\
t_4 := c \cdot \left(y4 \cdot t_3\right) - a \cdot \left(y5 \cdot t_3\right)\\
t_5 := x \cdot j - z \cdot k\\
t_6 := t_5 \cdot \left(b \cdot y0 - i \cdot y1\right)\\
t_7 := \left(a \cdot \left(y1 \cdot \left(z \cdot y3\right)\right) + c \cdot \left(y0 \cdot \left(x \cdot y2\right)\right)\right) - \left(c \cdot \left(y0 \cdot \left(z \cdot y3\right)\right) + y1 \cdot \left(a \cdot \left(x \cdot y2\right)\right)\right)\\
t_8 := t \cdot j - y \cdot k\\
t_9 := y4 \cdot \left(b \cdot t_8\right) - i \cdot \left(t_8 \cdot y5\right)\\
t_10 := k \cdot y2 - j \cdot y3\\
t_11 := t_2 - t_6\\
t_12 := t_11 + t_7\\
t_13 := \left(\left(t_12 + t_9\right) - t_4\right) + \left(y4 \cdot \left(y1 \cdot t_10\right) - y0 \cdot \left(y5 \cdot t_10\right)\right)\\
\mathbf{if}\;a \leq -1.7868532318154938 \cdot 10^{-151}:\\
\;\;\;\;t_13\\
\mathbf{else}:\\
\;\;\;\;\begin{array}{l}
t_14 := y1 \cdot y4 - y0 \cdot y5\\
t_15 := c \cdot y0 - a \cdot y1\\
t_16 := c \cdot y4 - a \cdot y5\\
t_17 := t_10 \cdot t_14\\
t_18 := t_8 \cdot \left(b \cdot y4 - i \cdot y5\right)\\
\mathbf{if}\;a \leq 4.7580947082678275 \cdot 10^{-307}:\\
\;\;\;\;\left(\left(\left(t_11 + \left(x \cdot y2 - z \cdot y3\right) \cdot t_15\right) + t_18\right) + y \cdot \left(y3 \cdot t_16\right)\right) + t_17\\
\mathbf{elif}\;a \leq 6.151744767147083 \cdot 10^{-235}:\\
\;\;\;\;t_13\\
\mathbf{elif}\;a \leq 3.904953834996993 \cdot 10^{-214}:\\
\;\;\;\;t_17 + \left(\left(t_18 + \left(t_11 - y3 \cdot \left(z \cdot t_15\right)\right)\right) - t_3 \cdot t_16\right)\\
\mathbf{elif}\;a \leq 6.503281304510843 \cdot 10^{-9}:\\
\;\;\;\;t_17 + \left(\left(t_9 + \left(t_7 + \left(t_2 + i \cdot \left(t_5 \cdot y1\right)\right)\right)\right) - t_4\right)\\
\mathbf{elif}\;a \leq 0.7947004650056086:\\
\;\;\;\;\left(\left(t_12 + t_18\right) - t_4\right) + \left(k \cdot y2\right) \cdot t_14\\
\mathbf{else}:\\
\;\;\;\;t_17 + \left(\left(t_9 + \left(t_7 + \left(\left(a \cdot \left(t_1 \cdot b\right) - c \cdot \left(t_1 \cdot i\right)\right) - t_6\right)\right)\right) - t_4\right)\\
\end{array}\\
\end{array}
(FPCore (x y z t a b c i j k y0 y1 y2 y3 y4 y5)
: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)))))(FPCore (x y z t a b c i j k y0 y1 y2 y3 y4 y5)
:precision binary64
(let* ((t_1 (- (* x y) (* z t)))
(t_2 (* t_1 (- (* a b) (* c i))))
(t_3 (- (* t y2) (* y y3)))
(t_4 (- (* c (* y4 t_3)) (* a (* y5 t_3))))
(t_5 (- (* x j) (* z k)))
(t_6 (* t_5 (- (* b y0) (* i y1))))
(t_7
(-
(+ (* a (* y1 (* z y3))) (* c (* y0 (* x y2))))
(+ (* c (* y0 (* z y3))) (* y1 (* a (* x y2))))))
(t_8 (- (* t j) (* y k)))
(t_9 (- (* y4 (* b t_8)) (* i (* t_8 y5))))
(t_10 (- (* k y2) (* j y3)))
(t_11 (- t_2 t_6))
(t_12 (+ t_11 t_7))
(t_13
(+ (- (+ t_12 t_9) t_4) (- (* y4 (* y1 t_10)) (* y0 (* y5 t_10))))))
(if (<= a -1.7868532318154938e-151)
t_13
(let* ((t_14 (- (* y1 y4) (* y0 y5)))
(t_15 (- (* c y0) (* a y1)))
(t_16 (- (* c y4) (* a y5)))
(t_17 (* t_10 t_14))
(t_18 (* t_8 (- (* b y4) (* i y5)))))
(if (<= a 4.7580947082678275e-307)
(+
(+
(+ (+ t_11 (* (- (* x y2) (* z y3)) t_15)) t_18)
(* y (* y3 t_16)))
t_17)
(if (<= a 6.151744767147083e-235)
t_13
(if (<= a 3.904953834996993e-214)
(+ t_17 (- (+ t_18 (- t_11 (* y3 (* z t_15)))) (* t_3 t_16)))
(if (<= a 6.503281304510843e-9)
(+ t_17 (- (+ t_9 (+ t_7 (+ t_2 (* i (* t_5 y1))))) t_4))
(if (<= a 0.7947004650056086)
(+ (- (+ t_12 t_18) t_4) (* (* k y2) t_14))
(+
t_17
(-
(+ t_9 (+ t_7 (- (- (* a (* t_1 b)) (* c (* t_1 i))) t_6)))
t_4)))))))))))double code(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) {
return (((((((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)));
}
double code(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 t_1 = (x * y) - (z * t);
double t_2 = t_1 * ((a * b) - (c * i));
double t_3 = (t * y2) - (y * y3);
double t_4 = (c * (y4 * t_3)) - (a * (y5 * t_3));
double t_5 = (x * j) - (z * k);
double t_6 = t_5 * ((b * y0) - (i * y1));
double t_7 = ((a * (y1 * (z * y3))) + (c * (y0 * (x * y2)))) - ((c * (y0 * (z * y3))) + (y1 * (a * (x * y2))));
double t_8 = (t * j) - (y * k);
double t_9 = (y4 * (b * t_8)) - (i * (t_8 * y5));
double t_10 = (k * y2) - (j * y3);
double t_11 = t_2 - t_6;
double t_12 = t_11 + t_7;
double t_13 = ((t_12 + t_9) - t_4) + ((y4 * (y1 * t_10)) - (y0 * (y5 * t_10)));
double tmp;
if (a <= -1.7868532318154938e-151) {
tmp = t_13;
} else {
double t_14 = (y1 * y4) - (y0 * y5);
double t_15 = (c * y0) - (a * y1);
double t_16 = (c * y4) - (a * y5);
double t_17 = t_10 * t_14;
double t_18 = t_8 * ((b * y4) - (i * y5));
double tmp_1;
if (a <= 4.7580947082678275e-307) {
tmp_1 = (((t_11 + (((x * y2) - (z * y3)) * t_15)) + t_18) + (y * (y3 * t_16))) + t_17;
} else if (a <= 6.151744767147083e-235) {
tmp_1 = t_13;
} else if (a <= 3.904953834996993e-214) {
tmp_1 = t_17 + ((t_18 + (t_11 - (y3 * (z * t_15)))) - (t_3 * t_16));
} else if (a <= 6.503281304510843e-9) {
tmp_1 = t_17 + ((t_9 + (t_7 + (t_2 + (i * (t_5 * y1))))) - t_4);
} else if (a <= 0.7947004650056086) {
tmp_1 = ((t_12 + t_18) - t_4) + ((k * y2) * t_14);
} else {
tmp_1 = t_17 + ((t_9 + (t_7 + (((a * (t_1 * b)) - (c * (t_1 * i))) - t_6))) - t_4);
}
tmp = tmp_1;
}
return tmp;
}




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
| Original | 26.8 |
|---|---|
| Target | 30.6 |
| Herbie | 27.1 |
if a < -1.78685323181549378e-151 or 4.7580947082678275e-307 < a < 6.1517447671470832e-235Initial program 27.0
Taylor expanded in x around inf 27.0
Applied sub-neg_binary6427.0
Applied distribute-rgt-in_binary6427.0
Simplified27.3
Simplified27.0
Applied sub-neg_binary6427.0
Applied distribute-rgt-in_binary6427.0
Simplified27.3
Simplified27.2
Applied sub-neg_binary6427.2
Applied distribute-rgt-in_binary6427.2
Simplified27.0
Simplified26.9
if -1.78685323181549378e-151 < a < 4.7580947082678275e-307Initial program 27.0
Taylor expanded in t around 0 28.4
Simplified28.4
if 6.1517447671470832e-235 < a < 3.9049538349969933e-214Initial program 28.0
Taylor expanded in x around 0 28.8
Simplified28.8
if 3.9049538349969933e-214 < a < 6.5032813045108433e-9Initial program 25.3
Taylor expanded in x around inf 26.4
Applied sub-neg_binary6426.4
Applied distribute-rgt-in_binary6426.4
Simplified26.7
Simplified28.3
Applied sub-neg_binary6428.3
Applied distribute-rgt-in_binary6428.3
Simplified28.5
Simplified28.5
Taylor expanded in i around inf 30.9
if 6.5032813045108433e-9 < a < 0.79470046500560865Initial program 24.7
Taylor expanded in x around inf 25.8
Applied sub-neg_binary6425.8
Applied distribute-rgt-in_binary6425.8
Simplified26.8
Simplified26.8
Taylor expanded in k around inf 28.2
if 0.79470046500560865 < a Initial program 28.1
Taylor expanded in x around inf 26.6
Applied sub-neg_binary6426.6
Applied distribute-rgt-in_binary6426.6
Simplified26.1
Simplified24.4
Applied sub-neg_binary6424.4
Applied distribute-rgt-in_binary6424.4
Simplified24.2
Simplified24.6
Applied sub-neg_binary6424.6
Applied distribute-rgt-in_binary6424.6
Simplified22.0
Simplified21.2
Final simplification27.1
herbie shell --seed 2022068
(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
: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.0 (- (* 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)))))